UNCLASSIFIED 


AD  NUMBER 

AD878776 

NEW  LIMITATION  CHANGE 
TO 

Approved  for  public  release,  distribution 
unlimited 


FROM 

Distribution  authorized  to  U.S.  Gov't, 
agencies  and  their  contractors;  Critical 
Technology;  OCT  1970.  Other  requests  shall 
be  referred  to  Air  Force  Flight  Dynamics 
Laboratory,  Attn:  FBE,  Wright-Patterson 
AFB,  OH  45433. 

AUTHORITY 

AFFDL  ltr,  12  Nov  1971 


THIS  PAGE  IS  UNCLASSIFIED 


fiDRo. -  AD8787  76 


AFFDL-TR*70-10 
Volume  I,  Part  ! 


LOW  ALTITUDE  ATMOSPHERIC  TURBULEI 
LO-LOCAT  PHASE  III 


Volume  I,  Part  I,  Data  Analysis 


/.  W.  JONES,  R  H.  M1ELKE,  G.  W.  JONES,  et  al 
THE  BOEING  COMPANY 


TECHNICAL  REPORT  AFFDL-TR-70-10,  VOLUME  I,  PART  I 


NOVEMBER  1970 


D  D  C 

m  sa  i.-;i 


■yJl  IT2 


This  document  is  subject  to  special  export  controls  and  each  transmittal  to  foreign 
governments  or  foreign  nationals  may  be  made  only  with  prior  approval  of  the  Air 
Force  Flight  Dynamics  Laboratory  (FBE),  Wright-Patic.s&.t  Air  Fore'  Case, 
Ohio,  45433. 


AIR  FORCE  FLIGHT  DYNAMICS  LABORATORY 
AIR  FORCE  SYSTEMS  COMMAND 
WRIGHT-PATTERSON  AIR  FORCE  BASE,  OHIO 


Best 

Available 

Copy 


NOTICE 


When  Government  drawings,  specifications,  or  other  data  are  used  for  any 
purpose  other  than  in  connection  with  a  definitely  related  Government 
procurement  operation,  the  United  States  Government  thereby  incurs  no 
responsibility  nor  any  obligation  whatsoever;  and  the  fact  that  the 
Government  may  have  formulated,  furnished,  or  in  any  way  supplied  the 
said  drawings,  specifications,  or  other  data,  is  not  to  be  regarded  by 
implication  or  otherwise  as  in  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rights  or  permission  to 
manufacture,  use,  or  sell  any  patented  invention  that  may  in  any  way  be 
related  thereto. 

This  document  is  subject  to  special  export  controls  and  each  transmittal 
to  foreign  governments  or  foreign  nationals  may  be  made  only  with  prior 
approval  of  the  Air  Force  Flight  Dynamics  Laboratory  (FBE)  }  Wright- 
Patterson  Air  Force  Base,  Ohio,  45433. 

The  distribution  of  this  report  is  limited  because  it  contains  informa¬ 
tion  that  would  significantly  diminish  the  technological  lead  time  of 
the  United  States  and  friendly  foreign  nations  by  revealing  formulas, 
processes,  or  techniques  having  a  potential  strategic  or  economic  value 
not  generally  known  throughout  the  world. 


win  tat m  q 
off  so iw 

cH 


Copies  of  thiB  report  should  not  be  returned  unless  return  is  required 
by  security  considerations,  contractual  obligations,  or  notice  on  a 
specific  document. 


.  (  O<0\  , 


LOW  ALTITUDE  ATMOSPHERIC  TURBULENCE 
LO-LOCAT  PHASE  III 


Volume  I,  Part  I,  Data  Analysis 


/.  W.  JONES.  R.  H.  MIELKE,  G.  W.  JONES,  et  al 


This  document  is  subject  to  special  export  controls  and  each  transmittal  to  foreign 
governments  or  foreign  nationals  may  be  made  only  with  prior  approval  of  the 
Air  Force  Right  Dynamics  Laboratory  (FBE),  Wright- Patterson  Air  Force  Base, 
Ohio  45433. 


FOREWORD 


This  is  th«  final  report  ou  the  Low-Low  Altitude  CAT  research  work  con¬ 
ducted  under  Contract  Number  F33615-68-C-1468  (LO-LOCAT  Phase  III).  The 
report  was  prepared  by  The  Boeing  Company,  Wichita  Division,  for  the  Air 
Force  Flight  Dynamics  Laboratory,  Wrlght-Patterson  Air  Force  Base,  Ohio. 

The  LO-LOCAT  Phase  III  project  was  part  of  Advanced  Development  Program 
682E  (ALLCAT)  and  was  under  the  direction  of  ADP  682E  Program  Director, 

Mr.  E.  Brazier,  and  the  Technical  Coordinator,  Mr.  Neal  V.  Loving. 

Mr.  Jan  Garrison,  FBE,  was  the  Air  Force  Project  Engineer. 

The  research  effort  was  conducted  under  the  Boeing  supervision  of  Mr.  F.  K. 
Atnip,  Program  Manager.  Mr.  C.  F.  Peterson  was  the  Project  Pilot.  Mr. 

L.  B.  Marshall  was  in  charge  of  instrumentation,  Mr.  H.  H.  Depew  directed 
the  data  processing  effort,  Mr.  J.  D.  Gault  was  in  charge  of  data  analysis, 
and  Mr.  W.  B.  Moreland  (Boeing-Seattle)  directed  the  meteorological  fore¬ 
casts  and  analysis.  Airplane  maintenance  and  inspection  were  the  respon¬ 
sibility  of  Mr.  J.  Strain  and  Mr.  J.  Bonawltz,  respectively. 

Mr.  Gerald  A.  Comstock  piloted  the  observer  airplane  over  the  Peterson 
Field,  Colorado,  high  fountain  route  and  was  also  backup  pilot  for  this 
program.  While  acting  in  this  latter  capacity  during  a  data  gathering 
flight  over  the  Griff las  AFB,  New  York  route,  Mr.  Comstock  was  fatally 
injured  in  a  forced  landing  following  an  engine  flameout. 

Authors  of  this  report,  other  than  those  shown  on  the  cover,  were  Messrs 
D.  E.  Gunter  and  K.  R.  Monson. 

The  report  was  submitted  by  the  authors  17  August  1970.  It  was  reviewed 
by  Mr.  Jen  Garrison  and  Dr.  T.  Swaney  (Boeing)  who  made  many  constructuve 
comments  and  suggestions. 

This  technical  report  haa  been  reviewed  and  is  approved. 


GORDON  R.  KEGAARJr,  Major,  USAF 
Chief,  Design  Criteria  Branch 
Structures  Division 
Air  Force  Flight  dynamics  Laboratory 


ii 


ABSTRACT 


/ 

/ 

/ 


This  report  presents  procedures,  analysis  methods,  and  final  results  per¬ 
taining  to  the  LO-LOCAI  Phase  III  program.  Approximately  150  hours  of  low 
altitude  (100  -  1000  feet)  turbulence  and  associated  meteorological  data 
were  recorded  from  l£  August  1968  through  30  June  1969*  The  original  pro¬ 
gram  was  curtailed  by  approximately  six  weeks  due  to  the  crash  of  the  T-33 
research  airplane.  A  model  of  the  turbulence  environment  at  low-level  Is 
presented  In  terms  of  gust  velocity  primary  peaks,  level  crossings,  soli¬ 
tude  samples,  rms  values,  and  gust  maxima,  as  veil  as  derived  equivalent 
gusts,  turbulence  scale  lengths,  and  power  spectra.  Mathematical  expres¬ 
sions  for  turbulence  spectra,  scale  length  and  primary  peak  statistics  are 
shown.  Correlations  between  atmospheric  gust  velocities  and  meteorological 
and  geophysical  phenomena  are  evaluated.  It  was  found  that  gust  velocity 
magnitude  at  low  altitude  Is  most  affected  by  atmospheric  stability  and 
terrain.  Gust  velocity  rms  values  above  1.5  fps  may  be  approximated  by 
truncated  Gaussian  distributions.  For  wavelengths  less  than  15,000  feet, 
turbulence  spectra  are  best  represented  by  the  von  Karmen  mathematical 
expressions.  The  turbulence,  Baapled  for  4-1/2  minute  intervals  over  a 
distance  of  approximately  32  miles  at  absolute  altitudes  below  1,000  feet, 
was  found  to  be  basically  stationary,  Isotropic,  and  homogeneous.  A  high 
percentage  of  Phase  in  data  were  recorded  over  high  mountains  since  very 
little  high  mountain  data  were  recorded  under  contour  flight  conditions 
at  low  level  during  Phases  I  end  H.  Phase  in  data  are  compared,  with 
data  from  Phases  I  and  H  and  with  data  from  other  low  altitude  programs. 

This  report  consists  of  two  volumes  with  each  volume  divided  Into  two  parts. 
Parts  I  and  II  of  Volune  I  give  the  techniques  and  results  of  data  analysis. 
Part  I  of  Volume  n  provides  the  details  pertaining  to  data  acquisition,  In¬ 
strumentation,  calibrations  and  checks,  data  processing,  and  data  quality. 
Data  tabulations  and  plots  and  a  log  of  pertinent  Information  concerning  the 
program  are  al  so  presented  in  Volume  n  Part  I.  Part  n  of  Volume  II  con¬ 
tains  the  power  spectral  density  and  other  frequency  data  plots. 

(Distribution  of  this  abstract  Is  unlimited.) 
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Symbols  and  Abbreviations : 

Ei  normalized  standard  error  of  a  regression  coefficient 
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f 

f(r) 
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fc 

u 

f. 

GMT 
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H 


Hc 
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Hg 
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Total  crossings  of  the  level  a. 

B^se  of  natural  logarithms  (2.71828). 

F -ratio  or  F-test  of  significance  in  regression 
analysis. 

Cumulative  probability  dirtfjbution  function. 
Frequency  (cycles/second); function. 

Longitudinal  correlation  function. 
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Cutoff  frequency  of  a  lev-pass  filter  ( cycles/ second). 
Syqulst  or  folding  frequency  (cycles/second). 
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Greenwich  Mean  Time. 

Ground  speed  (feet/second). 

Earth's  gravitational  constant  at  sea  level  (32.17^ 
feet/second)2. 

lateral  or  vertical  correlation  function.  * 

Altitude  above  the  earth's  surface  (feet);  true  alti¬ 
tude  (feet). 

Calibrated  pressure  altitude  (feet). 

Filter  transfer  function. 

Symbol  for  the  element  mercury. 

Indicated  reading  of  tower  altimeter  (feet). 

Indicated  pressure  altitude  (feet). 

Constrained  filter  transfer  fimction. 
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(cycles/foot) . 

H0x,  Nunfcer  of  text)  level  crossings. 

Np  Total  ousfcer  of  peaks  per  obtained  using  extrap¬ 

olation  technique  (characteristic  frequency). 

lfp.  Total  nusfcer  of  peaks  obtained  using  extrapolation 

technique. 

K.  Crossings  per  alle  of  the  level  a. 

n  Acceleration  In  g  wits;  counter;  nusfcer  of  data 

points;  shape  parameter  In  spectra  mathematical 
expressions. 
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On 

PSD 

PST 

P. 
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Outside  air  teaperature  (degrees  Fahrenheit). 

Rtai  test  observation  of  the  mean. 

Run  test  observation  of  the  mean  square. 

Power  spectral  density. 

pacific  Standard  Time. 

Static  pressure  (inches  of  mercury). 

Indicated  static  pressure  (inches  of  mercury). 

Standard  atmospheric  pressure  at  eea  level  (29.9a 
Inches  of  mercury). 

Intercepts  of  the  exponential  curves  for  the  analyt¬ 
ical  expressions  of  the  peak  court  probability  dis¬ 
tributions. 

Niaber  of  Independent  variables  in  the  regression 
equation. 

Calibrated  Impact  pressure  (inches  of  mercury). 
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q  Indicated  Impact  pressure  (inches  of  sercury). 

R*  Radar  altitude  (feet). 

R  Richardson  sober. 
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Rotation  otrlx  used  to  transform  aeasurenente  froa 
the  airplane  reference  area  to  the  earth  reference 


R(i) 
RMy (i) 


R(r) 

R(r)/R(o) 


r 


Value  of  m  gust  velocity  at  a  given  wind  speed  as 
read  froa  least  square  line  (feet/second). 

Autocorrelation  faction. 

Crass  correlation  function, 

a% —  - - —  -  - 

noot  assn  square. 

Oust  velocity’  covariance  faction. 

dust  velocity  autocorrelation  faction  where  r  a  o«  • 

•  *  ,Ti' 

Distance  (feet). 

Value  of  r  at  which  the  autocorrelation  function 
crosses  sera., 


8  Airplane  plsnfxm  wing  area  (feet2);  horizontal 

distance  -  wiles. 

8SS  Sine  of  solar  elevation. 

f.  Aabient  air  temperature  (degrees  Rscklne). 

Trut,  Clse  width  of  filter  weight  faction. 

Y«  (Around  surface  temperature  (degrees  Fahrenheit). 

Tj  Osjmxmted  outside  air  teagarature  (degrees  Rsnkine). 

T,  faaereture  at  sea  level  wader  standard  conditions 

(53B.69  degrees  fenUn*}. 

T#  Air  tewparatore  at  the  surface  (degrees  centigrade). 
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Symbols  and  Abbreviations : 


Temperature  at 

Ta2 

Temperature  at 

Ta3 

Temperature  at 

TaV 

Temperature  at 
Fahrenheit) . 

250  feet  (degrees  Fahrenheit). 

750  feet^( degrees  Fahrenheit). 

1000  feet  i degrees  Fahrenheit), 
flight  altitude  plus  1000  feet  (degrees 


T*7 

Tp2 

Td3 

Tdv 


Tenqsereture  at  flight  jalti£ude  (degrees  Fahrenheit). 

I 

Dew  point  at  750  feet  ^.degrees  Fahrenheit). 


Dew  point  at  1000  feet  (degrees  Fahrenheit). 

Dew  point  at  flight  altitude  plus  1000  feet  (degrees 
Fahrenheit). 


TD5 

TOD 

t 

u 

V 


Ve 


v 


Dew  point  at  850ob  (degrees  centigrade) . 

Time  of  day. 

Time  (seconds);  standardized  variable. 

Derived  equivalent  gust  velocity  (feet/second). 
Longitudinal  gust  velocity  (feet/second);  positive-aft. 
Trim  airspeed  or  ground  speed  (feet/second). 

Equivalent  airspeed  (feet /second). 

True  airspeed  (feet/ second) . 

Represents  the  matrix  of  true  airspeed  components 
corrected  for  pitch  and  yaw. 

Lateral  gust  velocity  (feet/second),  positive  to  the 

right. 


W  Airplane  weight  (pounds);  wind  speed  (feet/second) . 

Vk  Angle  between  airplane  ground  track  and  wind  vector 

(degrees)  positive  -  wind  vector  from  the  left,  zero 
degrees  -  direct  tail  wind. 

V*K  Wind  direction  at  flight  altitude  plus  1000  feet 

(degrees  azimuth) . 
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WJjqcL  direction  (degrees  azimuth). 

Average  easterly  winds  (feet/ second). 

Average  northerly  winds  (feet/second). 

Wind  speed  at  the  surface  (knots). 

Wind  speed  at  700mb  (knots). 

Vertical  gust  velocity  (feet/second),  positive  - 
upward. 

Amplitude. 

Mid-band  value. 

1  th  value  of  x;  predictors;  coefficient  of  simple 
linear  correlation;  value  of  gust  velocity  rns  at 
a  given  wind  speed  (feet/second). 

Level  of  gust  velocity  (feet/second). 

kth  record  Included  In  the  ensemble  averaging 

scheme. 

Sampled  valve  of  tUae  series. 

Dependent  variable;  regression  function. 

t 

*  Angle  of  attack  (degrees),  positive  -  nose  above 
relatlve  wjnd;  levels  of  gust  velocity  used  in 
\  level  crossing  procedures;  confidence  limit. 

Angle  of  s Idea lip  (degree*),  positive  -  none  left 
of  relative  wind. 

« 

Standard  lapse  rate  (0.00356  degrees  Fahrenheit /foot) 
positive  -  temperature  decrease  with  increasing  alti¬ 
tude. 

Dry  adiabatic  lapse  rate  (0.0055  degrees  Fahrenheit/ 
foot). 

Ratio  of  specific  heats  for  air  (l.to). 

Coherency  function.  % 
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Symbols  and  Abbreviations: 

A 

Incremental  or  difference. 

Af 

Change  In  frequency  (cycles/second);  frequency  range 
spanned  by  filter  roll-off;  frequency  interval  be¬ 
tween  adjacent  PSD  data  points. 

AH 

Distance  above  or  below  reference  point  (feet); 
difference  In  terrain  elevation  (feet). 

Ah 

Static  pressure  error  (feet). 

AwZco 

Incremental  load  factor  (acceleration  units). 

AP. 

Static  pressure  position  error  (inches  of  mercury). 

AP./ftj 

Pressure  coefficient  (dimensionless). 

AP„ 

Angle  of  attack  differential  pressure  (inches  of 
mercury). 

AP* 

Angle  of  sideslip  differential  pressure  (inches  of 
mercury). 

AT., 

Vertical  temperature  gradient  (degrees  Fahrenheit/ 
foot). 

ath 

Horizontal  temperature  gradient  (degrees  Fahrenheit/ 
mile). 

AT, 

Atmospheric  stability  (degrees  Fahrenheit/lOOO  feet), 
positive  -  temperature  increase  with  increasing  alti¬ 
tude. 

At 

Time  interval  between  data  samples  (seconds). 

awd< 

Vertical  wind  direction  gradient  (degrees/foot). 

AW„ 

Horizontal  wind  gradient  (feet/second/mile). 

Awiz  v 

Horizontal  east-west  wind  gradient  (feet /second/ 

^  mile). 

AW, 

Vertical  wind  velocity  gradient  (feet/second/foot). 

XXXV 


SYMBOLS 


Symbols  and  Abbreviations; 

i  Drift  angle  (degrees),  positive  to  the  right. 

2  •s 

«  Viscous  dissipation  rate  (feet  /second3);  endsslvlty 

factor;  filter  empirical  error  bounds  (percent). 

if  Kolmogorov  microscale  (feet). 

i  Pitch  angle  (degrees),  positive  -  nose  up. 

$  Rate  of  pitch  (degrees/second),  positive  -  nose 

moving  up. 


X 


M 


m2,m3 

Mt) 


r 

0 

P 


0 


*1 


*P« 


Taylor  turbulence  microscale  length  (feet);  character¬ 
istic  wavelength  (feet);  terrain  wavelength  (feet). 

Air  viscosity  (pound  second/feet2);  rasan  value;  air¬ 
plane  mass  ratio. 

Constant  in  gust  velocity  rms  distribution  equation. 
Ensemble  average  time  function. 

Kinematic  viscosity  (feet2 /second);  degrees  of  freedom. 
3.lUl6.  ... 

Air  density  (slugs/foot3). 

Standard  air  density  at  sea  level  (0.002378  slugs/ 
foot3) • 

v*1 

Standard  deviation  of  a  statistical  sanq>le;  standard 
deviation  (fps)  from  gust  velocity  spectra  between 
0.0416  mod  10  cps;  Stefhn-Boltzman  constant. 

Standard  deviation  about  a  least  square  line. 

Standard  deviation  of  level  crossing  distribution. 

Standard  deviation  of  noise. 

Standard  deviation  of  primry  peaks  obtained  using 
extrapolated  value  of 

Standard  deviation  of  recorded  data. 


Standard  deviation  obtained  from  the  truncated  gust 
velocity  spectra  (feet./second). 
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5  1b  and  Abbreviations: 


Standard  deviation  of  the  gust  velocity  time  aeries 
(feet /second) . 

Standard  deviation  of  terrain  roughness  (feet). 

Standard  deviation  of  angle  of  attack  differential 
pressure  (inches  of  mercury) . 

Standard  deviation  of  sideslip  differential  pressure 
(inches  of  mercury). 

Standard  deviation  (dispersion  of  the  distribution) 
of  gust  velocity  rms  (feet/second). 

Standard  deviation  of  the  derivative  of  the  time 
function  divided  by  2»V  (cps  for  gust  velocity) . 

Constant  used  in  rm&  gust  velocity  distribution 
equation.  & 


*  (k^ 

+  (n) 

*’  (n) 

+xy  (n) 


Lag  time  for  the  weighting  operation  convolution 
(seconds). 

One-dimensional  gust  velocity  power  spectral  density. 
Power  spectral  estimates. 

Raw  estimate  of  power  spectral  density. 

Averaged  estimate  of  power  spectral  density. 

Estimate  of  cross  power  spectral  density. 


Power  spectra  90  per  cent  confidence  interval. 

Roll  angle  (degrees),  positive  -  right  wind  down. 

Rate  of  roll  (degrees/ second),  positive  -  right  wing 
moving  down. 

Chi  statistic. 

Yaw  angle  (degrees),  positive  -  nose  right. 

Rate  of  yaw  (degrees/second),  positive  -  airplane 
nose  moving  right. 


Frequency  (radians/second). 
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Subscripts  and  Superscript b : 


b 

0 

e 

flit. 

H 

1 

N 

n 

K 


grin 

P 

T 

t 

u 

▼ 
i r 
x 

y 

z 

o 

I 

n 

m 


Band  number. 

Baaed  on  the  Dryden  equation. 

Extrapolated  value. 

Filtered  value. 

Horizontal. 

Semple  mzaber. 

Noise. 

Counter. 

Baaed  on  the  von  Kansan  equation. 

Indicates  maximum  value. 

Indicates  nHn<— «  value. 

Peak  count;  baaed  an  the  Lumley-Penofsky  equation. 

Truncated;  obtained  from  tower. 

From  the  gust  velocity  time  series. 

From  the  longitudinal  gust  velocity  component. 

From  the  lateral  gust  velocity  coeyonent. 

From  the  vertical  gust  velocity  component, 
longitudinal. 

Lateral. 

Vertical. 

Initial  value. 

Pertaining  to  quadrant  1;  Indicates  LO-LOCAT  Phases  I  and  II. 
Pertaining  to  quadrant  2. 

Pertaining  to  quadrant  3;  Indicates  LD-ICCAT  Phase  m. 
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SECTION  I 
INTRODUCTION 


Knowledge  of  the  characteristics,  causes,  and  effects  of  atmospheric  tur¬ 
bulence  has  become  incaShsingly  important  during  the  last  few  years.  In¬ 
dications  are  that  the  need  for  this  knowledge  will  continue  to  Increase 
in  the  future.  This  Is  due  to  several  factors.  Many  economic  and  safety 
considerations  are  related  to  Increased  aircraft  life  through  decreased 
fatigue'  failure.  The  effects  of  atmospheric  turbulence  on  the  handling 
and  ride  qualities  of  airplanes  have  become  increasingly  significant  as 
the  demand  for  Improved  performance  has  Increased. 

Older  generation  aircraft,  due  to  their  relatively  slow  speeds,  operated 
'  well  Into  the  inertial  subrange  of  atmospheric  turbulence.  As  a  result, 
they  had  to  cope  mainly  with  the  low  power,  high  frequencies  of  turbulence. 
Requirements  for  high  performance  have  brought  about  the  use  of  highly  flex¬ 
ible  structures  in  aircraft  design.  Increased  speed  has  caused  airplanes  to 
be  responsive  to  the  higher  power  long  wavelength  turbulence  components  while 
Increased  flexibility  has  decreased  structural  response  frequencies.  As  a 
result,  aircraft  rigid  body  and  structural  response  frequencies  are  coinci¬ 
dent  with  frequencies  of  atmospheric  turbulence  which  are  capable  of  signif¬ 
icant  driving  forces.  It  appears  that  future  aircraft  may  be  responsive  to 
turbulence  frequencies  In  the  frequency  range  where  atmospheric  turbulence 
power  spectral  density,  in  theory,  remains  at  a  large  constant  value  with 
further  decreases  In  frequency.  The  detailed  definition  of  turbulence  char¬ 
acteristics  in  this  realm  has,  therefore,  become  increasingly  pertinent. 

Military  missions  Include  terrain  contour  flight  at  low  altitude  to  avoid 
radar  detection.  In  this  region  of  the  atmosphere,  near  ground  level,  there 
1b  a  frequent  occurrence  of  high  magnitude  turbulence.  Therefore,  detailed 
knowledge  of  the  atmospheric  environment  at  low  altitudes  has  become  partic¬ 
ularly  important.  The  Air  Force  has  long  recogi ized  this  fact  and  in  their 
continuing  efforts  to  design  better  aircraft,  contracted  The  Boeing  Company 
to  participate  in  an  extensive  investigation  of  low  altitude  turbulence.  The 
research  program  was  designated  LO-IDCAT  (Low-Low  Altitude  Critical  Atmospheric 
Turbulence).  ID-IOCAT  was  established  as  part  of  an  overall  program  known  as 
AiUCAT  (Reference  I.l).  ALLCAT  was  devised  to  define  a  model  of  atmospheric 
turbulence  over  a  wide  range  of  altitudes  as  follows: 


Program 

Altitude  (Feet) 

TOLCAT 

0  - 

250 

IA-IOCAT 

250  - 

1,000 

LOCAT 

1,000  - 

20,000 

MEDCAT 

20,000  - 

40,000 

HICAT 

140,000  - 

70,000 

HI-HICAT 

70,000  - 

200,000 

1 


Research  of  the  lower  atmosphere  had  "been  accomplished  in  the  past.  These 
programs  were  worthwhile  contributions,  but  the  data  were  limited  in  quan¬ 
tity  and/or  resolution  of  long  turbulence  wavelengths.  Advances  of  instru¬ 
ment  technology  and  calibration  techniques  following  these  programs  made 
the  measurement  of  lover  frequencies  possible.  Improved  computers  made 
the  accomplishment  of  a  turbulence  research  program  baaed  on  statistical 
analysis  of  large  quantities  of  data  feasible. 

The  purpose  of  the  LO-IOCAI  Program  was  to  determine  the  turbulence  envi¬ 
ronment  below  1,000  feet  above  the  ground.  This  was  done  by  utilising 
statistically  representative  samples  of  turbulence  data  obtained  over  a 
vide  range  of  meteorological,  topographical,  seasonal,  and  time  of  day 
conditions.  Basically,  the  following  general  items  were  studied: 

e  Probabilities  of  encountering  given  levels  of  turbulence  inten¬ 
sity. 

•  Frequency  analysis  of  atmospheric  turbulence  utilising  the  power 
spectral  method. 

o  Correlation  of  atmospheric  turbulence  with  meteorological  and 
geophysical  conditions,  and  associated  forecasting  techniques. 

The  ID-LOCAT  Program  consisted  of  three  phases.  Phases  I  and  U,  which _ 

are  discussed  in  Reference  1.2,  involved  the  use  of  ftour  C-131B  aircraft 
as  instrumentation  platforsa.  Bach  of  these  aircraft  was  stationed  at  a 
different  base  within  the  United  States,  and  was  flown  over  a  specific 
route  laid  out  near  the  base.  The  routes  were  located  in  Bew  York,  Kansas, 
Colorado,  and  California  and  were  established  to  provide  a  vide  range  in 
topographical  and  climatological  conditions.  Data  were  gathered  for  ap¬ 
proximately  15  months. 

The  C-131B  aircraft  were  flown  at  an  average  airspeed  of  approximately 
330  fps.  The  instrumentation  had  adequate  frequency  response  down  to 
approximately  0.04  cps;  therefore,  wavelengths  up  to  8,000  feet  were 
investigated. 

Phase  ill  of  the  ID-IDCAT  Program  was  accomplished  in  order  to  extend 
the  statistical  definition  of  the  turbulence  environment  and  investigate 
longer  wavelengths.  A  T-33A  airplane  was  utilised  and  was  flown  at  an 
average  speed  of  approximately  630  fps.  This  higher  speed,  with  the 
same  basic  instrumentation  used  during  Phases  I  and  U,  made  possible 
the  measurement  of  turbulence  wavelengths  up  to  15,500  feet. 

The  Phase  III  data  were  recorded  over  approximately  a  10 -month  period. 

Data  were  gathered  at  the  same  geographical  locations  as  during  Phases 
I  and  n. 

\  * 

A  specific  route  for  the  airplane  to  follow  was  laid  out  at  each  base. 

Each  route  consisted  of  eight  straight  legs  approximately  45  statute 
miles  long.  All  legs  were  traversed  in  the  same  direction  on  each  flight. 
Pour  and  one-half  minutes  of  data,  recorded  while  flying  over  ea^h  leg, 
constituted  a  turbulence  sample.  Prior  to  -flying  each  leg,  meteorological 
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data  were  obtained  at  100  and  1,000  feet  above  the  terrain.  These  data 
were  obtained  in  the  same  vertical  air  mass  at  the  beginning  of  each  leg, 
and  were  used  to  conpute  vertical  gradients  of  temperature,  pressure,  and 
wind. 

Plights  were  scheduled  nearly  every  week  day  not  required  for  routine  main¬ 
tenance  of  the  airplane  and  instrumentation .  Three  missions  were  scheduled 
each  flight  day;  one  at  dawn,  'one  at  mid-morning,  and  one  •  in'  the  afternoon. 
The  airplane  was  maintained  and  flown  by  Boeing  personnel.  Two  specific 
altitudes,  250  feet  and  750  feet  above  the  terrain,  were  used  for  the  gust 
data  gathering  portion  of  the  flights.  Only  one  altitude  was  flown  on  a 
given  day.  The  pilot  followed  the  terrain  contour,  as  closely  as  safety 
allowed,  using  a  radar  altimeter  to  maintain  a  constant  absolute  altitude. 

This  report  primarily  covers  the  results  of  the  research  conducted  under 
Phase  HI  of  the  program.  Resmts  from  additional  research  of  the  long 
wavelengths  of  turbulence  and  a  comparison  of  Phase  III  data  with  the  data 
from  Phases  I  and  n  and  other  turbulence  programs  are  also  included.  Vol¬ 
ume  I  is  divided  into  two  parts  and  gives  the  techniques  and  results  of 
data  analysis.  Volume  II  is  also  divided  into  two  parts.  Part  I  provides 
the  details  pertaining  to  data  acquisition,  instrumentation,  calibrations 
and  checks,  data  processing,  and  data  quality.  Data  tabulations  and  plots, 
and  a  log  of  pertinent  information  concerning  the  program  are  also  presented 
in  Volume  II  Part  I.  Part  H  of  Volume  II  contains  the  power  spectral  den¬ 
sity  and  other  frequency  data  plots. 

Pertinent  data  obtained  during  KJ-IOCAT  Phase  III  were  recorded  on  mag¬ 
netic  tape.  These  tapes  were  delivered  to  the  Air  Force  Flight  Dynamics 
Laboratory. 

1.  IMfflSUMmATIOW 

The  T-33A-  airplane  was  instrumented  to  measure  meteorological  and  gust 
velocity  data.  A  radar  altimeter,  Doppler  radar  system,  radiometer,  and 
outside  air  temperature  probe  comprised  the  basic  instrumentation  used 
to  obtain  meteorological  data.  Atmospheric  turbulence  was  sensed  by  in¬ 
strumentation  contained  in  a  Boeing-designed  gust  probe  mounted  on  the 
end  of  a  nose  boom.  Probe  data  were  supplemented  by  airplane  attitude 
and  roll  rate  information  obtained  at  the  boom  base.  Instrumentation 
contained  inside  the  gust  probe  sensed  the  following: 

o  Sideslip  and  attack  angles 

e  Vertical,  lateral,  and  longitudinal  accelerations 

•  Static  and  impact  pressures 

•  Pitch  and  yaw  rates 

Fine  and  coarse  measurements  were  provided,  in  some  instances,  to  improve 
resolution.  Also,’  the  temperature  of  each  transducer  was  maintained  at 
135°  F  by  thermostatically  controlled  heaters. 

Details  of  the  Instrumentation  are  presented  in  Appendix. I. 


All  data  vere  categorised  with  respect  to  selected  geophysical  conditions 
under  which  they  were  obtained.  This  foies  a  ccmnan  basis  for  ccsfiarlaon 
of  the  data,  points  out  specific  effects  of  scan  of  the  acre  Important 
category  panaeten,  and  provides  a  systematic  aethod  of  mnalysla.  At 
the  beginning  of  the  W-WCK2  Program,  six  gamer si  categories  were  chosen 
In  an  anticipated  order  of  significance.  These  Include  type  of  terrain, 
altitude  shoes  the  terrain  (absolute  altitude),  atmospheric  stability, 
time  of  day,  season  of  the  year,  and  geographic  location.  These  general 
anas  wen  divided  Into  the  components  Indicated  in  the  following  list: 

Type  of  Terrain 

1.  High  Maintain 

2.  Low  Mountain 

3*  Desert 

*.  Plains 

5*  Sump  (not  used) 

6.  Met ear 

Absolute  Altitude 

1.  290  Feet 

2.  790  Feet 


Atmospheric  Stability 


1.  Very  stable 

2.  Stable 

3.  fcutrel 
*.  Unstable 


t$ti  >  -2*F/1jOQO  ft. 

-2*F/1000  ft.  *  Zff.  >  -5*f/lOOO  ft. 
-5*f/idoo  ft.  *  >  -6* F/lflOO  ft. 

ZCC*  <  -6*F/lOOO  ft. 


TSm  of  Day 

1.  Dawn 

2.  IO.d-anm1.ng 

3.  Kid-afternoon 


SNOOT  TUBS  -  LOCAL  TDK  ±  30  NDUXS 


Summer  flgrlag-lWll  Winter 

ducrlse  “SErfsie  Sunrise 

0900  0930  1000 

1500  1*30  1*00 


Season  of  the  Tear 


1.  Spring  -  March  21  to  June  20 

2.  tamer  -  J»  21  to  September  20 

3.  Fail  -  Sept  saber  21  to  December  20 
*.  Winter  -  December  21  to  March  20 


Oaognphlc  Location 

1.  Eduards  Air  Pare*  Base,  California 

2.  Qrlfflss  Air  Pone  Baas,  lew  Tort 

3.  Peterson  Field,  Colorado 

*.  McConnell  Air  Force  Base,  Tenses 


* 


t! 


The  number  beside  each  category  component,  or  element,  was  used  as  a  code 
for  that  particular  element.  All  data  recorded  during  a  turbulence  sample 
were  then  categorized  by  the  appropriate  six-digit  nunker.  The  following 
example  Illustrates  this  method: 

Category  No.  422324 

4  -  Terrain  Feature,  plains 
2  -  Altitude,  750  Feet 

2  -  Atmospheric  Stability,  Stable 

3  -  Time  of  Day,  Afternoon 
2  -  Season,  Svaner 

4  -  Location,  McConnell  AFB 

Some  of  the  data  within  different  categories  were  pooled  when  the  data 
shoved  like  characteristics.  This  was  done  to  improve  statistical  anal¬ 
yses  by  Increasing  sample  sizes.  For  exaqple,  suppose  that  data  within 
category  422324  shows  the  same  characteristics  as  data  within  categories 
422124  and  422224,  i.e. ,  time  of  day  showed  no  effect  an  the  data.  Then 
all  the  data  were  put  together  within  the  category  422024,  where  the  0 
indicates  that  all  times  of  day  are  included. 

All  category  makers,  except  stability,  were  assigned  prior  to  data  pro¬ 
cessing.  The  maker  for  atmospheric  stability  was  assigned  after  the 
vertical  temperature  gradient  had  been  computed  for  each  turbulence  sam¬ 
ple.  i 

3.  DMA  ANALYSIS 


The  IO-LOCAT  Phase  HI  data  were  edited  In  a  maker  of  v&yB  depending  an 
the  statistic  being  analyzed.  The  various  methods  of  editing  and  analysis 
are  discussed  In  each  section  of  this  document  along  with  the  particular 
type  of  d  ta  Involved.  Time  histories  of  Individual  parameters  measured 
during  each  test  wen  evaluated  to  determine  if  instrumentation  malfucc- 
tior;  occurred.  Data  involving  malfunctions  were  eliminated  for  all  anal¬ 
yse* ,  The  time  series  standard  deviations  of  angle  of  Bideslip  and  angle 
of  attack  measurements  were  computed  for  each  sample.  When  these  values 
Indicated  a  low  signal  to  noise  ratio  for  particular  turbulence  samples, 
the  turbulence  samples  were  case  handled  to  conserve  computer  time.  They 
were  not  eliminated  from  statistical  calculations. 

A  schematic  shoving  the  primary  steps  Involved  in  the  data  processing, 
editing,  and  analysis  functions  is  presented  In  Figure  3-1- 

Oust  velocity  statistics  were  determined  based  on  prisary  peaks,  amplitude, 
and  level  crossing  counts.  Ibe  distributions  of  gust  velocity  peak  and 
amplitude  values  were  compared  to  normal  distributions.  The  level  cross - 
lag  data  were  coopered  to  Rice's  equation  (Reference  3* l) .  Values  of 
characteristic  frequency  calculated  from  peak  and  level  crossing  count 
data  were  compared  with  thoee  calculated  from  the  power  spectral  density. 
Frequency  of  exceedance  and  the  probability  of  exceeding  given  levels  of 
gust  velocity  were  determined. 
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Gust  velocity  rt  — ard  deviations,  determined  from  peak,  amplitude,  and 
level  crossings  count  data  were  compared.  The  cumulative  probability  of 
these  data  was  calculated. 

Power  spectral  densities  of  the  gust  velocity  components  were  calculated 
and  plotted  for  those  turbulence  samples  recorded  over  the  even  numbered 
legs  on  the  even  numbered  flights  and  for  those  recorded  over  the  odd  num¬ 
bered^  legs  on  the  odd  numbered  flights. 

The^pectra  were  normalized  by  the  gust  velocity  variance  and  also  by 
L  Mathematical  representation  of  the  spectra  was  investigated  by  com¬ 
paring  the  experimental  spectra  with  the  commonly  used  mathematical  expres¬ 
sions  as  suggested  by  von  Karman,  Dryden,  and  others.  The  integral  scale 
of  t Turbulence,  essential  to  the  representation  of  the  spectra  by  mathe¬ 
matical  expressions,  was  analyzed.  Turbulence  qualities  such  as  homoge¬ 
neity  and  isotropy  were  also  evaluated. 

Although  important  parameters  are  included  in  the  category  components, 
there  were  other  possible  parameters  which  might  affect  the  turbulence 
intensity.  The  turbulence ,  recorded  during  Phase  III  was  compared  to  some 
of  these  parameters.  They  included  wind  speed,  wind  direction,  stability 
ratie,  Richardson  number,  air  temperature,  ground  surface  temperature,  and 
other  meteorological  parameters  obtained  from  rayinsonde  equipment,  ground 
observers,  and  r/nopotic  charts. 

Analysis  of  the  data  by  category  required  a  pooling  of  data  showing  like 
characteristics.  Obviously,  many  combinations  of  the  category  components 
were  possible.  If-  all  possible  combinations  were  used,  only  a  small  amount 
of  data  would  be  in  each  combination.  Therefore,  the  data  were  pooled  for 
category  components  which  had  the  smallest  effects  in  order  to  determine 
reliable'  statistics  of  the  data. 

Peak  count,  gust  velocity  power  spectra,  turbulence  scale  length,  and  gust 
velocity  rms  data  obtained  during  Phase  III  were  compared  to  correspond¬ 
ing  data  obtained  during  Phases  I  and  II.  Corresponding  data  from  the 
various  phases  could  not  be  Indiscriminately  pooled  together  because  of 
differences  in  the  conditions  under  which  they  were  obtained  as  discussed 
in  Section  IX.. 


4. ,  DATA  GAPS  HI  THE  TURBULENCE  M3  EEL 


L0-L0CA1  Phases  I  and  II  data  can  be  used  to  define  a  turbulence  model  for 
all  terrain  types  except  high  mountain,  for  all  times  of  day  except  night¬ 
time,  fpr  various  weather  conditions  except  storms,  for  contour  flight  at 
altitudes  of  250  and  750  feet,  and  for  at  least  four  locations. 

Turbulence  data  were  rot  recorded  under  contour  flight  conditions  at  low 
level  over  the  high  mountain  terrain  at  the  Peterson  route  during  Phases 
I  and  II  as  shown  in  Figure  4.1.  Turbulence  'data  were  not  recorded  at 
nighttime  nor  during  periods  of  inclement  weather.  Although  the  data 
were  recorded  at  each  location  during  all  seasons,  the  climatology  for  “ 
that  particular  year  may  not  have  been  representative  of  the  average. 

It  was  the  purpose  of  Phase  III  to  extend  the  statistical  definition 
obtained  during  Phases  J.  and  II  and  to  define  longer  wavelengths. 

1  * 


r -■ 


Figure  4.1  Hata  Gaps 


Segments  of  data  significant  to  the  normal  atmospheric  syBtem,  for  which 
adequate  samples  for  a  statistical  model  were  not  obtained  are  defined  as 
data  gaps. 

The  best  method  to  develop  a  turbulence  model  based  on  LO-LOCAT  data  would 
be  to: 

,  •  Evaluate  the  variation  of  Phases  I  and  II  data  with  season. 

These  statistics  are  shown  in  Table  4.1. 

•  Assume  that  corresponding  seasonal  variations  would  have  oc¬ 
curred  during  Phase  HI  and  expand  the  data  (with  the  excep¬ 
tion  of  high  mountain  data)  accordingly. 

•  Evaluate  the  relationship  between  high  mountain  data  at 
Peterson  and  pertinent  climatology  which  existed  at  the  time 
it  was  gathered.  The  geophysical  situation  during  Phase  III 

,  is  discussed  in  Section  VI. 

•  Estimate  the  Peterson  high  mountain  data  in  the  gap  based  on 
the  climatology  which  normally  would  exist  during  the  time 
period  of  the  gap.  Records  of  normal  climatology  are  avail¬ 
able  from  the  (J.S.  Weather  Bureau. 

Once  the  Phase  IH  data  have  been  applied/  other  considerations  must  be 
made.  These  include: 

•  How  well  the  geophysical  situations  encountered  during  LO-LOCAT 
match  those  for  the  development  of  a  particular  model. 

•  The  compatibility  of  the  range  of  turbulence  wavelengths  mea¬ 
sured  with  that  of  an  airplane  being  designed. 

•  The  effects  of  not  gathering  data  under  storm  conditions  dur¬ 
ing  LO-LOCAT. 

Of  the  three  other  considerations  mentioned,  the  first  is  probably  the 
least  significant.  The  LO-LOCAT  Program  was  designed  to  cover  a  large 
number  of  geophysical  situations.  Chances  are,  the  geophysical  condi¬ 
tions  under  which  it  is  desired  to  develop  a  turbulence  model  may  be  well 
matched.  Not  a  great  amount  of  data  were  obtained  over  water  but  from 
all  indications,  low  altitude  atmospheric  turbulence  over  water  at  wave¬ 
lengths  less  than  15,000  feet  is  relatively  insignificant. 

The  difference  between  the  range  of  turbulence  wavelengths  measured  and 
the  range  to  which  an  airplane  being  designed  would  respond  is  a  highly 
significant  consideration.  It  is  necessary  to  extrapolate  the  data  avail¬ 
able  in  an  attempt  to  arrive  at  an  estimate  for  the  range  in  question. 
Effects  of  shifting  the  wavelength  window  thorough  which  the  turbulence 
is  viewed  is  discussed  throughout  this  report. 
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TABIE  4.1 


statistics  of  oust  velocity  intensity  (»t)  with 

SEASON  AMD  LOCATION  DURING  PHASES  I  AND  II 


Category 

Oust  Velocity  Component 

' 

yjfps) 

v(; 

Pps. 

1 

At* _ 

* 

*71 

100011 

3.14 

0.93 

3.38 

0.95 

3.61 

1,04 

100021 

3.21 

1.32 

3.59 

1.53 

3.13 

I.09 

100031 

2.93 

1.29 

3.17 

1.33 

3.00 

1.38 

100041 

2.78 

1.31 

3.05 

i.37 

2.74 

1.28 

200011 

2.78 

0.88 

3.09 

o.97 

3.20 

1.00 

200021 

2.74 

0.92 

2.77 

1.09 

3.22 

1.12 

200031 

2.52 

1.15 

2.83 

1.29 

2.71 

1.25 

200041 

2.69 

1.14 

3.07 

1.28 

2.72 

1.17 

300011 

2.72, 

0.99 

2.90 

1.02 

3.20 

1.22 

300021 

2.32 

1.01 

2.17 

1.00 

2.92 

1.28 

300031 

1.64 

0.86 

1.90 

0.92 

1.91 

1.11 

300041 

1.55 

0.70 

1.91 

0.79 

1.83 

1.00 

600011* 

0.93 

0.4l 

1.00 

0.30 

0.82 

0.31 

600021* 

0.56 

0.18 

0.72 

0.23 

0.65 

0.27 

600031 

0.89 

0752 

1.08 

0.56 

0.90 

0.56 

600041* 

0.6l 

0.25 

0.82 

0.26 

0.59 

0.26 

200012 

3.76 

1.46 

4.1C 

1.68 

3.76 

1.28 

200022 

3.03 

1.07 

3.22 

1.17 

3.21 

1.11 

200032 

3.00 

1.28 

3.31 

1.47 

2.99 

L.10 

200042 

3.21 

1.32 

3.59 

1.53 

3.13 

1.09 

400012 

3.67 

1.49 

4.07 

1.60 

3.43 

1.27 

400022 

2.78 

0.93 

2,82 

0.92 

2.76 

0.94 

400032 

2.81 

1.25 

3.02 

1.27 

2.6l 

1.04 

400042 

3.08 

1.31 

3.39 

1.44 

2.8l 

1.05 

400013 

2.46  1 

1.05 

2.95 

1.23 

2.71 

1.33 

400023 

2.14 

0.96 

2.48 

1.03 

2.40 

1.20 

400033 

1.89 

0.86 

2.21 

1.00 

1.84 

0.91 

400043 

2.43 

0.94 

2.77 

1.00 

2.40 

1.00 

400014 

2.80 

1.19 

2.98 

1.24 

2.73 

L.ll 

400024 

2.08 

0.84 

2.05 

1.05 

2.15 

0.91 

400034 

2.30 

1.02' 

2.49 

1.04 

2.22 

0.93 

400044 

2.79 

. 

1.08 

3.02 

1.15 

2.48 

0.97 

*  Less  than  30  sasqples  Involved 

NOTE:  ^  ■  mean  of  the  standard  deviations 

»  *  standard  deviation  of  the  gust  velocity  time  series 
standard  deviations 
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The  third  consideration  is  also  inqportant.  if  a  complete  model  is  being 
developed,  the  data  must  be  adjusted  to  account  for  the  fact  that  LQ-LOCAT 
data  were  not  measured  during  adverse  weather.  The  extent  of  the  adjust¬ 
ment  required  is  given  in  Tables  4.2  and  4.3.  These  data  were  taken  from 
Reference  4.1.  Turbulence  data  used  in  the  adjustment  must  come  from  pro¬ 
grams  other  than  LO-IDCAT.  Two  such  sources  are  References  4.2  and  4.3. 

A  multiplicity  of  types  of  data  were  obtained  and  analyzed  during  LO-UDCAT. 
These  data  may  be  applied  in  mmerous  vays  to  the  evaluation  of  atmospheric 
turbulence  at  low  altitude,  depending  on  the  application,  of  ap¬ 

plications  and  adjustments  that  may  be  made  are  shown  in  Sections  ix  and  X 
of  this  report.  In  these  sections,  the  LO-LOCAT  Phased  III  data  are  ana¬ 
lyzed  with  data  from  Phases  I  and  II  and  with  data  from  other  low  altitude 
turbulence  programs. 


TABUS  4.2 


HUMBER  OF  IFR  WEATHER  DAYS 


Location 

Spring 

Days 

Stainer 

Days 

Fall 
_  Days 

— 

Winter 

Days 

McConnell 

<9 

<9 

<9 

9-18 

45  (12.3) 

Edwards 

9-18 

9  -  36 

s 

i 

o\ 

9-27 

99  (27.1) 

Peterson 

<9 

<9 

<9 

<9 

36  (9.85) 

Griffias 

9  -  18 

'  <9 

9-27 

18  -  27 

8l  (22.2) 

TABLE  4.3 

NUMBER  OF  THUHLEBaTORM  DAYS 


Location 

Spring 

Days 

Summer 

Days 

Fall 

Winter 

Days 

Maximun 

Yearly 

Day®  (*) 

McConnell 

10 

20 

10 

1-5 

45  (12.3) 

Edwards 

1 

1 

1 

1 

(1.D 

Peterson 

10 

30 

5 

0  -  1 

k6  (12.6) 

Oriffiss 

5 

10 

1-5 

0-1 

21  (5.75) 

Oriffiss 
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SECTION  II 


SUMMARY 


It  was  found,  based  on  primary  peak,  amplitude,  and  level  crossing  counts 
and  rms  values  that  gust  velocity  magnitude  becomes  larger  as  the  terrain 
becomes  rougher,  as  altitude  above  the  terrain  is  decreased,  and  as  the 
atmosphere  becomes  less  stable.  Statistical  distributions  of  gust  veloc¬ 
ity  primary  peak,  amplitude^  and  level  crossing  data  were  found  to  be  in 
reasonable  agreement  with  each  other.  (See  Figures  11.10  through  11.12). 

The  distributions  of  the  primary  peaks  were  found  to  be  definable  by  an 
equation  using  the  component  Gaussian  process.  The  distribution  of  the 
gust  velocity  rms  values  was  also  defined  by  using  the  combination  of  two 
Gaussian  equations  (see  Figures  13*35  through  13. 37). 

Gust  velocities  in  excess  of  'JO  feet  per  second  were  encountered.  The 
maximum  value  for  vertical  gust  velocity  was  74.5  fps  and  occurred  over 
leg  3  of  the  Edwards  route.  The  maximum  gust  velocity  for  longitudinal 
and  lateral  components  were  71*0  and  76.9  fps,  respectively,  and  occurred 
over  leg  5  of  the  Peterson  route.  The  probability  distribution  of  the 
maximum  derived  equivalent  vertical  gust  velocity  agreed  well  with  that 
of  the  maximum  true  vertical  gust  velocity. 

The  nondlmenslonallzed  power  spectral  density  of  atmospheric  turbulence 
did  not  show  any  variation  which  might  be  associated  with  different  geo¬ 
physical  or  meteorological  phenomena.  Good  agreement  between  nondimen- 
Bionalized  von  Kaiman  and  measured  spectra  was  exhibited.  The  von  Kerman 
spectra,  however,  did  slightly  underestimate  the  .power  of  the  measured 
spectra  at  the  low  values  of  kL,  the  nondimensionalized  spatial  frequency. 
The  average  value  of  the  measured  scale  lengths  was  600  feet;  this  average, 
however,  is  more  representative  of  the  turbulence  found  around  the  high 
mountains  over  the  Peterson  route  in  Colorado  than  of  turbulence  in  general. 

Forecasting  of  low  altitude  turbulence  intensity  was  most  successful  when 
based  on  vertical  wind  gradient  and  lapse  rate.  It  was  found  that  the 
greatest  chance  of  encountering  significant  turbulence  exists  when  the 
Richardson  number  is  in  the  ±0.2  range.  As  the  Richardson  number  in¬ 
creases,  the  chance  of  high  intensity  turbulence  decreases  rapidly.  The 
stability  ratio  was  found  to  provide  a  good  indication  of  turbulence  in¬ 
tensity,  but  did  not,  as  suggested  by  other  investigators,  provide  any 
advantages  over  the  Richardson  number. 

Taylor  and  Kolmogorov  microscales  and  viscous  dissipation  rates  and  sta¬ 
bility  ratios  were  evaluated.  They  were  compared  to  values  obtained  by 
other  investigators. 

LO-IDCAT  Phases  I,  II,  and  III  data  were  compared.  Basic  information  con¬ 
cerning  the  three  phases  is  shown  in  Table  II. 1.  The  Phase  III  data  tended 
to  have  higher  turbulence  magnitudes.  Hiifl  increased  magnitude  is  attrib¬ 
uted  to  the  measurement  of  longer  wavelength  turbulence  and  measurement  of 
turbulence  over  high  mountains.  An  example  of  measuring  longer  wavelength 
turbulence  is  shown  in  Figure  17.1  where  the  vertical  gust  velocity  com¬ 
ponent  contains  a  wavelength  of  approximately  14,300  feet  with  a  magnitude 


exceeding  74.0  fps.  The  wavelength  appears  at  a  frequency  of  approximately 
0.04  cps  In  the  time  history.  This  same  long  wavelength  turbulence,  if  en~ 

•  countered  during  Phases  I  and  n,  would  have  appeared  at  a  frequency  of  0.02 
cps  due  to  the  lower  speed  at  which  the  C-131  airplanes  were  flown.  Phases 
I  and  II  gust  velocity  measurements  in  this  frequency  realm  would-  have  been 
attenuated  approximately  80  per  cent  by  the  data  filtering  techniques  used 
to  alleviate  the  effect  of  drift.  This  is  illustrative  of  the  fact  that 
there  is  turbulence,  some  of  which  is  significantly  high  in  intensity,  that 
contains  wavelengths  longer  than  7000  feet. 

a 

\ 

LO-LOCAT  data  were  compared  to  that  from  other  low  altitude  turbulence  pro¬ 
grams.  Because  of  the  variety  of  filtering  techniques  used,  only  those 
sauries  which  were  most  agreeable  to  Kolmogorov's  theory  were  used.  Based 
oh  the  method  of  selecting  the  spectra,  it  was  concluded  that  the  spectra 
of  turbulence  at  low  altitudes  from  various  independent  researchers  are 
consistent  in  that  their  shape  can  be  approximated  by  the  von  Harman  math¬ 
ematical  expressions  with  scale  lengths  generally  less  than  1000  feet. 

% 

During  the  course  of  this  program,  three  additional  turbulence  investiga¬ 
tions  were  accomplished.  These  were  the  BRER  Tower  Flyby,  Thunderstorm 
Turbulence  Investigation,  and  the  Wake  Turbulence  Investigation.  The  BREN 
Tower  Flyby  investigation  was  accomplished  to  compare  turbulence  data  re¬ 
corded  from  an  instrumented  airplane  and  an  instrumented  tower.  Although 
the  turbulence  levels  were  slightly  different,  the  normal.  1  zed  spectra  were 
in  good  agreement.  The  Thunderstorm  Turbulence  Investigation  was  accom¬ 
plished  to  investigate  turbulence  near  thunderstorms  and  convective  clouds. 
Unfortunately,  no  severe  thunderstorms  developed  in  the  area  during  the 
period  of  investigation.  It  was  found  that  turbulence  intensities  near 
convective  clouds  were  low  and  that  the  power  spectra  of  turbulence  near 
convective  clouds  exhibits  a  -5/3  logarithmic  slope  relationship  betveon 
power  density  and  spatial  frequency.  Air  temperature  changes  near  cumulus 
clouds  occurred,  in  some  instances,  when  turbulence  was  encountered.  The 
Wake  Turbulence  Investigation  consisted  of  flying  the  Phase  III  airplane 
behind  formations  of  C-l4l  airplanes  in  an  effort  to  qualitatively  evalu¬ 
ate  the  wake  turbulence  environment  generated  by  these  airplanes.  The 
maximum  gust  velocity  encountered  was  63  fps  and  the  estimated  size  of 
the  vortices  varied  from  19  to  4l  feet  in  diameter.  Although  high  gust 
velocities  were  encountered,  the  time  duration  of  the  gusts  were  so  short 
that  the  T-33  alrplatoe  did  not  respond  to  the  turbulence  and  the  pilot  re¬ 
ported  very  Little  turbulence. 

Since  the  major  findings  of  LO-LOCAT  Phase  III  were  the  significant  influ¬ 
ences  on  atmospheric  turbulence  statistics  of  the  longer  wavelengths  of 
turbulence  and  turbulence  over  high  mountains  in  Colorado,  additional  re¬ 
search  was  accomplished.  Gust  accelerations  and  wind  velocities  were  cal¬ 
culated  and  analyzed  to  determine  if  the  longer  wavelengths  could  be  better 
resolved.  The  most  pertinent  of  extreme  gust  data  obtained  during  the  High 
Intensity  Gust  Investigation  was  reprocessed  using  techniques  developed  dur¬ 
ing  the  LO-IOCAT  program  to  see  if  the  new  techniques  would  alter  the  results 
originally  obtained. 

Spectra  of  longitudinal  and  lateral  wind  components  showed  indications  of 
isotropy  for' the  two  components.  Integral  scale  lengths  were  found  to  be 
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,  on  the  order  of  6000  feet.  A  technique  w  developed  whereby  the  etatletica 
of  guet  velocities  can  be  determined  using  gust  accelerations*  Shis  tech¬ 
nique  proved  to  have  two  distinct  advantages.  The  gust  accelerations  re¬ 
quire  less  high-pass  filtering  due  to  their  lew  power  at  low  frequencies. 

The  technique  also  requires  considerably  less  computer  effort  since  filter¬ 
ing  asy  be  performed  in  the  frequency  rather  than  the  time  dosmin.  Analysis 
of  the  reprocessed  High  Intensity  Oust  Investigation  data  showed  that  the 
original  reaulte»were  effected  very  little  by  the  I0-LOGAT  processing  tech¬ 
niques. 

It  is  believed,  that  enough  statistical  Information  on  low  altitude  turbu¬ 
lence  is  available  from  various  programs  such  that  adequate  gust  criteria 
can  be  established  for  aircraft  which  are  sensitive  to  wavelengths  less 
then  15*000  feet.  In  the  application  of  these  data  to  gust  criteria  for 
the  mission  requirements  of  a  given  airplane  the  following  should  be  con¬ 
sidered: 

e  The  difference  in  speed  at  which  the  IO-IOCAT  data  were  obtained 
and  that  associated  with  the  mission.  Turbulence  data  recorded 
at  higher  speeds  will  usually  have  the  longer  wavelengths  better 
defined.  This  causes  the  gust  velocity  probability  distributions* 
average  scale  lengths*  etc.  to  vary  as  a  function  of  speed. 

e  The  geophysical  conditions  under  which  the  data  were  obtained. 
Turbulence  data  for  wavelengths  less  than  8,000  feet  were  ob¬ 
tained  during  rhaeea  I  and  H  for  all  saaaons  at  different 
locations*  for  all  terralsi  types  except  the  high  mountain*  for 
a  wide  range  of  meteorological  conditions  except  those  involv¬ 
ing  precipitation  and  storms*  and  fbr  a^l  times  of  day  exclud¬ 
ing  darkness,  phase  HI  data  were  not  recorded  during  all  sea- 
eons  at  each  of  the  route  locations. 


TftBUE  n.i 

basic  ifl-iflCAT  moaracB 


Phases  I  A  n 

Phase  HI 

Aircraft  used  to  obtain  data 

k  C-131's 

IT-33 

Flight  period 

9-15-66  -  12-20-67 

8-16-68  -  6-30-69 

Humber  of  turbulence  data 
flights 

1,244 

299 

lew  altitude  turbulence  samples 

8,871 

1*938 

BBS  tower  ewplee 

0 

8 

Thunderstorm  turbulence  staples 

0 

55 

ttoke  turbulence  esaplee 

0 

46 
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TABIE  H.l  (Cantd.) 
BASIC  LO-LOCAT  1HF0BMATZOH 


Phases  I  A  II 

Phase  in 

Turbulence  sample  average 
length 

21  ml. 

32  mi. 

* 

Turbulence  san$le  average 
time  Interval 

5.5  min. 

4.5  min. 

Average  turbulence  PSD  wave¬ 
length  band 

33  -  7,960  feet 

63  -  15,100  feet 

Aircraft  average .  ground  speed 

332  fps 

630  fps 

*t  bandwidth* 

0  -  0#  cps 

0  -  »  cps 

0 1  bandwidth 

0.33  -  10  cps 

.6&  -  10  cps 

9  bandwidth 

0.04  -  10  cps 

0.04  -  10  cps 

Hot  wire  data  bandwidth 

5  -  500  cps 

n/a 

u  time  series  samples  analyzed 
statistically 

6,553 

1,730 

v  time  series  sauries  analyzed 
statistically 

6,270 

1,709 

w  time  series  samples  analyzed 
statistically 

6,508 

1,716 

u  non-contour  sample;* 

1,113 

0 

v  non-contour  saiqplea 

1,091 

0 

w  non-contour  saaples 

1,116 

0 

V 

u  me  values  analyzed 
statistically 

7,670 

1,762 

v  rms  values  analyzed 
artistically 

6;62i 

1,740 

* 

v  rms  values  analyzed 
statistically 

7,630 

1,746 

Oust  velocity  samples  for  which 
spectra  were  calculated 

1,775 

904 

*  Influenced  below  0.046  and  above  10.0  cpa  by  high  end  low  pass  filters. 
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TABLE  II. 1  (Concluded) 
BASIC  LO-LOCAT  IPCBMITCT 


Phases  I  A  II 

Phase  ni 

u  spectra  eliminated  due  to 

296 

313 

irregularities** 

v  spectra  eliminated  due  to 

514 

380 

irregularities** 

v  spectra  ellainated  due  to 

306 

4o8 

irregularities** 

Homogeneous  u  spectra  analysed 

1,272 

591 

Homogeneous  v  spectra  analyzed 

l,06l 

524 

Homogeneous  v  spectra  analyzed 

1,259 

4g6 

u  spectra  obtained  during  non- 

207 

0 

contour  flight 

v 

J 

v  spectra  obtained  during  non- 

200 

0 

contour  flight 

v  spectra  obtained  during  non- 

210 

0 

contour  flight 

Oust  velocity  time  series  samples 

33,000 

27.000 

per  turbulence  sample 

Time  Interval  -  spectra 

0.01  sec. 

0.02  sec. 

calculations 

Oust  velocity  samples  in  spectra 

33,000 

13,500 

calculations 

Frequency  interval  -  spectra 

0.04l£  cps 

• 

0.04l6  cps 

calculations 

Hyquist  frequency 

50  cps 

25  cps 

ffimber  of  frequency  estimates  - 

1,200 

600 

PSD 

Miles  covered  recording  accept- 

135,900 

55,750 

able  turbulence  data 

Hours  of  acceptable  turbulence 

600 

130 

data 

Spectrum  degrees  of  freedom  - 

55 

45 

Individual  samples 

— . .  . —  ■ 

**  Low  slgnal/nolse,  high  airplane  notion  effect*,  low  turbulence  homogeneity 


etc. 
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SECTION  ni 


GUST  VELOCITY  TIME  FUNCTION 


5.  RUN  REST  FOR  SIATIONAlRITY 

One  of  the  assumptions  normily  nade  In  the  analysis  of  turbulence  data  Is 
that  It  Is  stationary.  Stationarlty  generally  means  that  statistical  pro- 
,  perties  computed  over  a  short  interval  of  time  do  not  vary  more  than  would 
be  expected  due  to  normal  "T1  variations. 

In  order  to  verify  that  the  stationarlty  assumption  was  valid,  it  was  nec¬ 
essary  to  perform  a  test  for  stationarlty  on  each  turbulence  sample.  Dur¬ 
ing  LO-LOCAT  Phase  m,  the  run  test  procedure  for  stationarlty,  as  described 
in  Reference  5. 1,  Ws  used.  Two  separate  run  tests,  one  on  the  .gust  velocity 
mean  values  and  one  on  the  mean-square  values,  were  performed  for  each  turbu¬ 
lence  sample.  The  use  of  mean-square  values  was  suggested  in  Reference  5. 2. 
Run  test  observations  were  obtained  using  the  following  equations: 


(Mean)  (5.1) 

*  v' 

b  -  1 

0'„  -  X*  (Mean-Square)  (5-2) 

i-a 

where: 

a  .  135u  -  134 
b  •  135n 

n  »  1,  3,  5.  ..  199 

A  total  of  lOO  observations  were  made  for  each  run  test.  Each  observation 
*  (0n  or  0'B).  das  compared  with  the  observation  mean  and  assigned  a  sign* 
plus  when  above  and  minus  when  below  the  mean.  The  number  of  runs  is  de¬ 
fined  as  the  mmtoer  of  sets  with  like  signs.  The  number  of  runs  and  whether 
or  not  the  maple  was  accepted  ss  stationary  at  the  0.02  level  of  signifi¬ 
cance  for  both  run  test  methods  was  then  determined.  This  was  done  by  com¬ 
paring  the  actual  number  of  runs  to  the  expected  number  of  r  vus  at  a  given 
level  of  significance. 

The  expected  n\mber  of  runs  at  a  level  of  significance  of  0.02  was  determined 
using  the  following  equations: 


* 


♦  1 


(5.3) 


»i  +  N2 


where: 


-.1 

2N1B2(2H1H2-H1-M2) 
(Hx  +  »2>2(1*i+M2-1) 

m 


(fU) 


Nx  =  Total  number  of  plus  values  (Ojj  or  0^) 

N2  m  Total  Number  of  minus  values  (0q  or  0^) 

Using  the  above  information,  the  expected  nuniber  of  runs  (E)  was  then  cal¬ 
culated: 

E  a  (j|  ±  2.326  s)  -0.5 

where  2.326  corresponds  to  the  0.02  level  of  significance. 

When  the  actual  number  of  runs  fell  between  the  two  values  of  E,  the  sta- 
tionarity  hypothesis  was  accepted  at  the  0.02  level  of  significance. 

Run  test  data  were  obtained  for  all  turbulence  samples  recorded  after  Flight 
Number  177.  Due  to  a  computer  programming  problem,  run  ?«.  ts  were  obtained 
for  only  some  of  the  flights  prior  to  Number  178.  Run  tests  were,  however, 
performed  on  a  total  of  002  turbulence  samples. 

The  results  for  both  run  tests  performed  during  L0-L0CAT  Phase  III  are 
shown  in  Table  5.1.  This  table  shows  that  over  00  per  cent  of  the  turbu¬ 
lence  samples  tested  by  a  run  test  of  the  means  were  accepted  as  stationary. 
The  mean-square  test  accepted  from  60  to  70  per  cent  of  the  samples  depend¬ 
ing  upon  the  component  of  gust  velocity  being  tested.  This  indicates  that 
the  mean-square  test  is  the  more  stringent.  The  average  of  these  tests  for 
all  components  indicat*  that  77  per  cent  of  the  samples  may  he  accepted  as 
stationary. 

When  a  table  such  as  the  one  in  Reference  $.1  1b  used  to  determine  the  ex¬ 
pected  number  of  runs  corresponding  to  a  given  level  of  significance,  the 
assumption  that  Nx  =*  in  Equations  5*3  and  5.4  is  made.  The  results  of 
making  this  assumption  are  illustrated  in  Table  5-1.  As  shown,  the  per¬ 
centage  of  turbulence  samples  accepted  as  stationary  is  reduced  by  using 
the  Nx  «  N?  assumption.  The  mean-square  test  is  affected  much  more  than 
the  test  of  the  means. 


A  number  indicating  the  degree  of  atationarity  for  each  turbulence  sample 
was  not  computed.  Hoover,  since  homogeneity  is  the  property  which  in 
spree  corrdinates  corresponds  to  the  stationarlty  property  in  time  coor¬ 
dinates,  these  two  properties  may  be  considered  similar.  Numerical  valves 
indicating  the  degree  of  homogeneity  were  computed  for  turbulence  samples 
for  which  power  spectra  were  computed  and  are  presented  in  Section  21  and 
in  Appendix  VTII. 
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TABLE  5.1 


PERCENT  OF  TUHBUIZNCE  SAMPLES  ACCEPTED  AS  SEATIONARY 
AT  TEE  0.02  LEVEL  OF  SIGNIFICANCE 


(Dest  of  Means 

Test  of  Mean  Squares  | 

u 

V 

w 

u 

V 

w 

Analytical 

Method 

85.1 

80.7 

88.3 

77.^ 

6o.l 

70.7 

Table 

Method 

77-0 

71.2 

79-0 

42.0 

28.5 

45.5 

6.  ENSEMBLE  AVERAGING 

A  random  process  may  be  classed  as  either  nonstationary,  stationary,  or 
stationary  and  ergo  die.  Atmospheric  turbulence,  over  long  periods  of  time, 
must  be  considered  to  be  nonstat ionary.  Stationarity  Implies  that  statis¬ 
tical  properties  computed  over  short  time  intervals  do  not  vary  more  than 
would  be  expected  due  to  normal  sampling  variations.  The  fact  that  turbu¬ 
lence  is  nonstationary  is  seen  by  comparing  mean  square  values  computed 
from  turbulence  samples  for  different  days.  Also,  in  mountainous  regions 
where  the  terrain  is  not  homogeneous  and  the  wind  currents  vary  considerably, 
mean  square  values  for  adjacent  4.5  minute  turbulence  samples  have  been 
found  to  vary  from  3  to  75  (fps)2  . 

Over  a  small  area  and  a  short  time  interval,  turbulence  can  usually  be  con¬ 
sidered  to  be  stationary.  Atmospheric  turbulence  cannot  be  considered  to 
fall  into  the  stationary  and  ergo die  class  since  ergo dicity  implies  that 
each  record  of  an  ensemble  is  statistically  equivalent  to  every  other  record 
and  that  ensemble  averages  over  a  large  number  of  records  at  fixed  times  may 
be  replaced  by  corresponding  time  average  on  a  single  representative  record 
of  the  ensemble. 

Reference  6.1  indicates  that  to  obtain  a  veil  defined  statistical  represen¬ 
tation  of  a  random  process  that  is  nonstationary  or  stationary  and  noty 
ergodic,  ensemble  averages  should  be  computed.  Since  atmospheric  turbu-^ 
lence  could  fall  into  either  class,  ensemble  averages  were  calculated  for 
several  categories  of  data.  This  averaging  of  a  number  of  samples  gave 
high  statistical  confidence,  as  discussed  in  Section  20. 

Ensemble  averaging  was  accomplished  by  taking  an  ensemble  (group)  of  records 
and  averaging  the  ensemble  at  corresponding  times.  Individual  points  of  the 
ensemble  average  time  function  j»,(t)  were  calculated  by: 

(6-1J 

k-1 
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Each  of  the  270-second  records  used  were  considered  to  be  self  stationary. 

t 

Ensemble  average  time  functions  were  confuted  for  each  of  the  eight  fol¬ 
lowing  category  combinations: 

High  mountain,  250  foot,  stable  (112000) 

High  mountain,  250  foot,  neutral  ( 113000) 

High  mountain,  250  foot,  unstable  (ll4000) 


High  mountain, 
High  mountain, 
High  mountain, 


750  foot,  stable  (122000) 
750  foot,  neutral  (123000) 
750  foot,  unstable  (121*000) 


Plains,  250  foot,  unstable  (414000) 
Plains,  750  foot,  unstable  (1*2 1*000) 


Power  spectra  and  coherency  functions  were  calculated  for  the  ensemble  aver¬ 
age  time  functions.  Isotropy  and  homogeneity  characteristics  were  evaluated 
and  the  spectra  were  compared  to  the  von  Karmen  mathematical  expressions. 

These  data  are  presented  in  Figures  6.1  through  6.8.  Isotropy,  homogeneity, 
and  coherency  checks  gave  essentially  the  same  results  as  those  from  the 
principle  analyses  (Section  v).^  The  observed  spectra  agreed  well  with  the 
von  Karman  spectra.  !£he  scale  lengths  from  the  ensemble  time  function,  with 
the  exception  of  two  points,  agreed  fairly  well  with  those  representing  the 
average  of  f.he  scale  lengths  for  the  individual  samples.  These  data  are  pre¬ 
sented  in  Figure  6.9.  The  ensemble  average  scale  lengths  varied  with  category 
in  the  same  way  that  irdividuf.1  sample  scale  lengths  varied.  That  is,  gener¬ 
ally  speaking,  scale  lengths  increase  with  increasing  terrain  roughness,  in¬ 
creasing  altitude  and  decreasing  stability. 

It  was  concluded  that  short,  relatively  self-stationary  segments  of  turbu¬ 
lence  data  can  be  accurately  evaluated  statistically  without  resorting  to 
ensemble  averaging  techniques. 
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Figure  6.1  IVmr  Spectra  Md  Associated  Data  of  Eos  cable 
Are  rage  Tlae  Foactioa  for  High  Mountain, 

250  Feet,  Stable  Category 


21 


_ _  _ A. ji  i  ijl  IM  J _ 


22 


»* 


Figure  6.3  Rover  Spectre  aid  Associated  Data  of  Quenble 
Average  Tine  Function  for  High  Mountain, 

230  Feet,  Unstable  Category 
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Figure  6.4  Poeer  Spectra  and  Associated  Data  of  Biicablc 
Average  Tie#  Function  for  High  Mountain, 

750  Feet,  Stable  Category 
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figure  6.3  Po»«r  Spectra  ad  Associated  Data  of  Enseablc 
Average  Tiae  Function  for  High  Mountain, 

730  Feet,  Neutral  Category 


Flgwr--  6.6  Spectra  m*  A**oclat*d  data  of  tumble 

A**<Kf*  Tim  hactlak  for  High  Mountain  t 
750  Fwt,  CttUhl*  Category 
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Figure  6,7  Power  Spectre  and  Associates  Data  ot  Diseeble 
A verace  Tlee  Function  lor  Plains,  CSC  Feet, 
Unstable  Category 
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Figure  6.8  Rower  Spectra  and  Associated  Data  of  &isaeble 
Average  Tine  Function  for  Plains ,  TjO  Feet , 
Unstable  Category 


LK  -  Feet  (average  of  individual  samples) 


Figure  6.9  Scale  Length  Coapariaon 


GAUSSIAN  DISTRIBUTION  CHECKS 


Gust  velocity  amplitude  values  computed  for  each  turbulence  sample  were 
checked  for  normality  using  the  chi-squared  goodness-of-flt  test.  (This 
test  compares  the  actual  distribution  with  a  normal  distribution. 

A  chi-squared  test  was  performed  on  each  four  and  ore-half  minute  turbu¬ 
lence  sample  using  28  bands  (Rj, )  with  widths  of  constant  probability.  The 
degrees  of  freedom  for  each  saiqple  were  (l^-l).  Two  hundred  (N)  gust  ve¬ 
locity  amplitude  valises  were  used  in  the  calculations  of  values  by  equa¬ 
tion  7,1: 


x2.£(v_m; 


(7.1) 


where: 


fb  =  observed  frequency 
Fb  -  theoretical  frequency 


Since: 


fri  *  fti  *  b  *b  • 


X2  - .lk  £  t£  -  200 


(7.2) 


A  cumulative  probability  distribution  of  the  chi-sqjare  vulues  for  all  L0- 
LOCAT  Phase  III  data  is  shown  in  Figure  7.1.  bey®}  =t  of  sienifi^ance  based 
on  the  number  of  degrees  of  freedom  for  these  data  are  shown.  Inspection 
of  this  figure  showB  that  approximately  87  per  curt  uf  the  vertical,  88  per 
cent  of  the  longitudinal,  and  89  per  cent  of  the  lateral  gust  velocity  sam¬ 
ples  were  accepted  as  Gaussian  at  the  0.02  level  of  significance.  Approxi¬ 
mately  10  per  cent  more  samples  were  accepted  during  Phases  I  and  U  (Ref. 
1,2).  15118  is  attributed  to  the  fact  that  turbulence  samp  baa  were  taken 
over  a  60  per  cent  longer  distance  ( approximately  32  udleo;  during  Phase  m. 
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TIVE  PROBABILITY 


figure  7.1 


CM -Square  Test  for  Statistical  Noreality 


8.  PRIMARY  PEAK  COURT 


In  the  peak  count  procedure  each  maximum  positive  or  negative  excursion  of 
a  gust  velocity  time  history,  between  adjacent  crossings  of  the  mean,  is 
defined  as  a  primary  peak.  A  graphical  illustration  of  this  procedure  is 
presented  in  Figure  8.1.  As  shown  in  this  figure  secondary  peaks  are  not 
considered,  and  only  the  primary  peaks  which  fall  within  each  2  tpa  wide 
band  are  counted. 

In  this  procedure,  high  frequency,  low  amplitude  data  are  counted  as  pri¬ 
mary  peaks  in  the  lowest  band.  Figure  8.2  shows  a  peak  count  distribution 
plotted  on  semi -logarithmic  scales  and  illustrates  how  it  was  extrapolated 
to  obtain  the  peaks  in  the  lowest  band.  The  shaded  points  indicate  the 
cumulative  number  of  peaks  in  the  lowest  band  as  they  were  originally 
counted.  The  unshaded  points  at  zero  gust  velocity  are  the  extrapolated 
values.  The  extrapolation  technique  had  the  effect  of  compensating  the 
error  induced  by  not  counting  the  secondary  peaks  and  thereby  adjusted 
the  cumulative  probability  of  the  peak  count  data  to  be  more  compatible 
with  that  of  the  level  crossings  data. 

The  extrapolation  was  performed  by  fitting  a  curve  through  the  data  points 
using  the  least  square  method.  The  curve  fitting  was  accomplished  by  com¬ 
puter  and  consisted  of  fitting  the  general  polynomial  (Equation  8.1)  to  the 
distribution. 

Yi  =  A0  +  A2  (8.1) 

The  coefficients  A0  and  A2  were  determined  such  that  the  following  quantity 
was  a  minimum: 

I  .2 

=  E  wi  fYi  <8-2> 

i-i  L 


where : 


Xt  =  (xt  +  20 )  in  fps,  where  is  the  gust  velocity 

Y±  =  log  (Fbi),  where  Fb  is  the  cumulative  number  of  peaks  per 
mile  occurring  in  a  given  band 


W  i  =  The  weight  applied  to  each  Y  j.  For  this  application  the 
following  values  were  used: 

wx  =  o  w6  =  6 

W2  =  20  W7  =  4 

w3  =  15  w8  =  2 

W{,  =  10  w9  through  =  1 

W5  =  8 
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e  =  error  or  difference  between  observed  and  analytical  vuluss 
of  Y 

N  -  Number  of  data  points  being  fitted 
A0,  A2  =  Constants 

The  following  equations  were  obtained  by  taking  the  partial  derivatives  of 
f(x)  with  respect  to  each  coefficient: 


I 

A0  E 

8 

+  A2  E  Wi  Xt 

8 

-  E Wi  Y*  - 0 

(8.3) 

1-1 

1-1 

i-i 

n 

8 

8 

Ao  E  wi 

X1  +  A2  E  W1 

xi  -  E  wi  xi  Yi  *  ° 

(8.4) 

1-1 

1-1 

i-i 

Simultaneous  solution  of  these  equations  yields  values  of  the  coefficients 
for  use  in  Equation  8.1.  The  coefficient  A0  is  the  extrapolated  log(r'p,e). 
Details  of  the  least  squares  curve  fit  are  given  in  References  8.1  through 
8.3. 


After  positive  and  negative  peaks  were  added  for  each  band  and  accumulated, 
cumulative  peaks  per  mile  were  calculated  by  dividing  the  cumulative  number 
of  peaks  in  each  band  by  the  miles  traveled  while  obtaining  the  turbulence 
sample.  Cumulative  nonliability  was  calculated  by  dividing  the  cumulative 
number  of  peaks  for  each  band  by  the  cumulative  number  of  peaks  for  the 
lowest  band  (for  peak  count,  this  is  an  extrapolated  value).  The  proba¬ 
bility  density  is  the  change  in  probability  from  one  band  to  the  next  higher 
one  divided  by  the  band  width. 

Figure  8.2  is  an  example  of  a  peak  count  cumulative  distribution  for  a  par¬ 
ticular  turbulence  sample.  The  curves  shown  on  this  plot  represent  those 
fitted  by  computer  to  the  data  points.  The  characteristic  frequency  of  the 
peak  count  data,  ,  is  defined  as  the  cumulative  peak  per  mile  value  de¬ 
termined  by  extrapolating  the  cumulative  distribution  curve  to  zero  gust 
velocity.  Figures  8.3  and  8.4  are  examples  of  peak  count  cumulative  prob¬ 
ability  distribution  and  probability  density  distribution  for  the  same  tur¬ 
bulence  sample  as  shown  in  Figure  8.2. 

The  turbulence  saddles  on  which  the  peak  count  was  performed  all  had  means 
which  had  been  set  to  zero.  A  standard  deviation  was  calculated  about  this 
mean  for  each  gust  velocity  Bangle  using  the  grouped  data  technique  dis¬ 
cussed  in  Reference  8.4.  Standard  deviation  values,  presented  in  Section  11, 
were  computed  from  the  following  equation: 


The  subscript  e  denotes  that  an  extrapolated  value  for  cunulative  peaks,  as 
discussed  previously,  was  used  in  this  calculation. 
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Those  turbulence  samples  set  aside  from  the  normal  data  processing  routine 
because  of  low  signal  to  noise  rat if j  'iiocvssed  in  Appendix  III,  Section 
III. 5)  were  included  in  the  gust  velocity  peak  count  distributions  so  that 
the  effects  of  low  intensity  turbulence  would  be  present  in  the  final  anal¬ 
ysis.  Inclusion  of  these  turbulence  samples  was  accomplished  in  the  fol¬ 
lowing  manner,  (l)  The  category  number  associated  with  each  of  these  tur- 
bulence  samples  was  determined  (Table  8.1  Is  a  list  of  the  number  of  samples 
involved  with  respect  to  category  number).  (2)  One  low  signal  to  noise 
sample  was  processed  for  each  different  category  number  to  obtain  peak  count 
data  representative  of  that  category.  (3)  The  representative  data,  multi¬ 
plied  by  a  factor  equal  to  the  number  of  samples  Involved  in  each  category 
(Table  8.1),  were  added  to  the  peak  count  distributions  initially  computed 
in  the  normal  data  processing  routine. 

Figure  8.5  is  a  plot  of  the  peak  distribution  for  the  overall  category. 

The  shaded  points  illustrate  the  effect  of  including  low  intensity  turbu¬ 
lence  samples.  Only  the  points  at  zero  gust  velocity  were  changed  a  no¬ 
ticeable  amount;  all  other  points  were  essentially  unchanged.  Also,  the 
Inclusion  of  the  low  intensity  turbulence  samples  slightly  lowered  the 
cumulative  probability  distribution  curves  of  those  data  in  the  categories 
involved.  Figure  8. 5-  also  compares  the  peak  count  distribution  curves  for 
the  gust  velocity  components  for  the  overall  category.  Lateral  gust  veloc¬ 
ity  has  the  highest  distribution  curve  and  longitudinal  the  lowest  with 
vertical  crossing  over  between  them.  The  relationship  between  components 
of  gust  velocity  for  all  categories  analyzed  Is  shown  in  Appendix  VH. 

Those  categories  which  contain  high  mountain  data  show  a  relationship  be¬ 
tween  components  similar  to  that  for  the  overall  category.  Categories  which 
do  not  contain  high  mountain  data  generally  show  a  closer  comparison  between 
lateral  and  vertical  distributions  with  longitudinal  always  slightly  lower. 
The  reason  lateral  is  higher  for  the  high  mountain  data  is  probably  due  to 
the  effects  of  flying  parallel  to  and  on  the  leeward  side  of  high  mountain 
ridges. 

The  distributions  of  the  three  components,  for  the  categories  other  than 
high  mountain,  agreed  much  more  closely  during  Phase  III  than  they  did  dur¬ 
ing  Phases  I  and  II.  High-pass  filtering  attenuates  the  long  wavelength 
amplitudes  of  the  longitudinal  campon.- >rt  more  than  the  vertical  or  lateral 
because  the  longitudinal  component  contains  higher  power  at  the  long  wave¬ 
lengths.  The  attenuation  of  these  long  wavelengths  was  less  severe  during 
Phase  III  than  during  Phases  I  and  n.  This  occurred  because  the  Phase  III 
data  were  recorded  using  a  higher  speed  airplane  which  made  possible  the 
measurement  of  longer  wavelengths  as  discussed  elsewhere  in  this  report. 
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TABUS  8.1 


HUJCER  OF  DOW  INTENSITY  TUHBUIZNCE  SAMPLES  IBOIEDED 
IN  EACH  CATEGORY  FOR  STATISTICAL  ANALYSIS  OF 
PEAK,  AMPLTTUIE,  AND  LEVEL  CROSSINGS  COUNT 


Category 

Number 

*  1 1  nl  ftl ;  laBPfci 

filial  WSMj, 

321131 

1 

611241 

1 

411124 

6 

611331 

5 

421113 

1 

611341 

1 

421124 

15 

621131 

6 

421134 

8 

62ll4l 

2 

421234 

3 

621231 

5 

422234 

1 

621241 

3 

424124 

1 

621331 

4 

6lll4l 

2 

621341 

1 

The  peak  count,  amplitude  count,  and  level  crossings  count  distributions 
obtained  for  each  turbulence  sample  were  compiled  on  a  master  computer 
tape  as  they  were  processed.  These  distributions  were  recorded  on  the 
master  tape  according  to  category  number.  For  DQ-DOCAT  Phase  HI,  data 
were  stored  in  26o  out  of  a  possible  1£20  different  categories  on  this 
tape.  Some  of  these  categories  contained  many  samples,  others  had  only 
one.  In  order  to  obtain  distributions  which  were  statistically  reliable, 
it  was  necessary  to  combine  categories  so  that  none  contained  fewer  than 
30  samples. 

The  master  tape  data  were  processed  through  a  computer  program  which  had 
the  capability  of  combining  categories  in  any  selected  combinations. 

Those  samples  not  categorized  with  respect  to  atmospheric  stability  be¬ 
came  of  Invalid  outside  air  temperature  data  could  not  be  handled  by  the 
computer  program  and  therefore,  were  not  Included  an  the  master'  tape,  nor 
included  in  the  statistical  analysis.  The  combined  category  distributions 
obtained  fron  this  program  were  then  plotted.  Table  8.2  lists  those  cate¬ 
gories  which  contained  30  or 'more  samples.  The  total  number  of  samples 
included  in  each  category,  as  well  as  the  number  of  low  intensity  turbu¬ 
lence  samples  included,  in  the  total  are  also  included  in  Table  8.2.  Plots 

and  listings  or  peak  count  distributions  and  probability  distributions  for 
these  categories  are  given  in  Appendices  VH  and  VHI,  respectively. 

Analytical,  expressions  of  the  peak  count  probability  distributions  were  de¬ 
termined  for  each  of  those  categories  involving  individual  category  compo¬ 
nents  as  well  as  the  *.n  data  category.  These  expressions  were  determined 
using  a  technique  similar  to  that  described  in  References  8.5  and  8.6.  This 
technique  involved  the  selection  by  least-squares  fit  techniques,  of  two  ex¬ 
ponential  curves  which,  when  added  together,  would  result  in  a  curve  repre¬ 
senting  the  distribution  in  question.  This  expression  takes  the  following 
form: 

F(x)  =  Px  exp(-  £_)  +  P  exp(-  i_)  (8.6) 

b  1  b2 
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where: 


F(x)  =  Cumulative  probability  distribution  function 
Pw  P2  =  Intercepts  of  the  exponential  curves 
l/b^  l/b2  =  Slopes  of  the  exponential  curves  on  semi-logarithmic  paper 
x  =  Gust  velocity  amplitude 

Table  8.3  shows  the  values  of  Pi,  P2,  bx,  and  b2  for  the  above  mentioned 
categories.  The  values  are  valid  out  to  the  maximum  gust  velocity  value 
measured  for  each  category  as  given  in  Appendix  VIII.  Experimental  data 
corresponding  to  the  fitted  curves  is  given  in  Appendix  m. 

Confidence  limits  were  computed  for  the  categorized  primary  peak  count 
cumulative  distribution  functions.  Since  each  distribution  analyzed  con¬ 
tained  many  turbulence  samples,  each  point  on  the  distribution  curve  is 
representative  of  the  mean  at  that  point  of  all  samples  contained  in  the 
category.  The  central  limit  theorem  states  that  for  large  samples,  mean 
values  are  distributed  normally  without  consideration  of  the  underlying 
distribution.  The  confidence  limits  were  therefore  calculated  using  the 
procedure  for  a  normal  population  with  unknown  standard  deviation,  as  de¬ 
scribed  in  Reference  8.4.  Peak  count  distributions  and  their  confidence 
limits  for  all  categories  analyzed  are  given  in  Appendix  I X.  Figures  8.6 
through  8.11  illustrate  the  relationship  of  the  99  and  95  per  cent  con¬ 
fidence  limits  to  the  distribution  of  all  of  Phase  III  data. 

An  analysis  of  the  peak  count  data  with  respect  to  those  categories  involv¬ 
ing  individual  category  components  as  well  as  the  all  data  category  is  con¬ 
tained  in  Section  11.  A  comparison  of  the  various  counting  techniques  is 
also  contained  in  that  section. 
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TABLE  8.2 


TOTAL  NUMBER  OF  TURBULENCE  SAMPLES  IN  EACH  COMBINED 
CATEGORY  DISTRIBUTION  AND  THE  NUMBER  OF  LOW 
INTENSITY  TURBULENCE  SAMPLES  INCLUDED  - 
PEAK,  AMPLITUDE,  AND  LEVEL  CROSSINGS  COUNT 


■■ 

Total 

Number  of  Samples 

Number  of 

Low  Turbulence 
Samples 

u 

V 

V 

OOOOOO 

1730 

1709 

1716 

66 

100000 

846 

84l 

846 

0 

200000 

297 

297 

297 

0 

300000 

115 

115 

115 

1 

400000 

396 

380 

382 

35 

600000 

76 

76 

76 

30 

010000 

899 

887 

891 

15 

020000 

831 

822 

825 

51 

001000 

561 

558 

560 

64 

002000 

479 

475 

477 

1 

003000 

393 

387 

390 

0 

oo4ooo 

297 

289 

289 

l 

000100 

1*65 

462 

i<65 

42 

000200 

644 

637 

639 

13 

000300 

621 

6lO 

612 

11 

000001 

599 

599 

599 

31 

000002 

135 

134 

135 

0 

000003 

743 

736 

742 

1 

000004 

253 

240 

240 

34 

niooc 

109 

108 

109 

0 

112000 

l4o 

l4o 

l4o 

0 

113000 

129 

128 

129 

0 

n4ooo 

64 

64 

64 

0 

121000 

1 12 

112 

112 

0 

122000 

147 

lte 

'-47 

0 

123000 

95 

94 

95 

0 

124000 

50 

49 

50 

0 

211000 

60 

60 

60 

0 

212000 

31 

31 

31 

0 

213000 

33 

33 

33 

0 

214000 

33 

33 

33 

0 

221000 

65 

65 

65 

0 

312000 

31 

31 

31 

0 

411000 

52 

51 

52 

6 

1-12000 

39 

39 

39 

0 

413000 

57 

54 

55 

0 

4l4ooo 

50 

44 

44 

0 

421000 

55 

54 

54 

27 

422000 

33 

30 

31 

1 

423000 

39 

38 

38 

0 

424000 

71 

70 

69 

1 

621000 

32 

32 

32 

21 
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Band 

No. 

Gust 

Velocity 

Plus  Peak 
Ident. 

Minus  Peak 
Ident. 

Cumulative 

Peaks 

Cumulative 
Peaks  per 
Mile* 

Probability 

Density 

■ 

0 

© 

7 

35 

1.0000 

.0714 

5 

2 

© 

© 

6 

30 

.8571 

.1428 

■ 

4 

mm 

4 

20 

.5714 

.0715 

■ 

6 

© 

®  © 

3 

15 

.4285 

.2143 

SOTS:  This  is  only  an  illustration  and  should  not  be  considered  as  typical  of 
turbulence  data. 

(2)  (6)  ^9)  «  Secondary  Peaks  -  Not  Counted 

*  For  this  exaaqple  ,  assisie  0.2  Btatute  miles  is  the  length  of  the  above  record. 


Figure  8.1  Primary  Peak  Count  Example 
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O  Longitudinal  55735 

□  Lateral  55076 

O  Vertical  55510 

CATEGORY  000000 


Shaded  points  are  those 
which  changed  due  to 
inclusion  of  low  turbulence 
saaples.  All  other  points 
were  essentially  unchanged. 
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f  lgur*  8.6  Londtttdlul  Oust  Velocity  FMk  Count  Distribution 
with  99i  Confidence  Units 
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TIVE  PEAKS  PER  STATUTE  MILE 


CUMULATIVE  PEAKS  PER  STATUTE  MILE 


TIVE  PEAKS  PER  STATUTE  MILE 


TIVE  PEAXS  PER  STATUTE  MILE 


AHPLITUIE  COURT 


The  amplitude  count  technique  of  analysis  vas  used  on  the  tine  series  of 
the  gust  velocity  coaponenta.  Values  for  each  component  were  computed  by 
the  gust  velocity  progrsn  (appendix  V)  at  ICO  saaples  per  second.  In  this 
procedure,  amplitude  hands  2  fpa  vide  were  placed  on  either  side  of  the 
zero  mean.  The  mmber  of  samples  contained  In  each  hand  vas  determined 
and  the  samples  in  the  corresponding  positive  and  negative  bands  were 
added  together.  Ibis  procedure  la  Illustrated  In  Figure  9.1.  Amplitude 
data  were  used  to  conpute  emulative  distributions,  emulative  probabil¬ 
ity  distributions  and  probability  density  distributions.  The  technique 
used  to  calculate  these  distributions  Is  also  ahovn  in  Figure  9.1.  The 
time  series  standard  deviation  (sj  for  each  gust  velocity  component  vas 
computed  from  the  following  equation: 


US'"-*’1]' 


(9.1) 


Examples  of  amplitude  count  emulative  distribution,  cumulative  proba¬ 
bility  distribution,  and  probability  density  distribution  for  an  indi¬ 
vidual  turbulence  saqple  are  ahovn  in  Figures  9.2,  9.3 ,  and  9.4,  respec¬ 
tively. 


Amplitude  data  were  combined  for  the  same  categories  shown  In  Table  8.2 
and  the  data  from  turbulence  samples  of  low  Intensity  were  accounted  for 
In  the  same  manner  as  described  in  Section  8. 


Amplitude  count  cumulative  probability  distribution  plots  for  the  Phase 
HI  data  axe  ahovn  In  Figure  9.9.  plots  for  the  k2  categories  of  Table 
8.2  are  shown  in  Appendix  VH  and  listings  of  these  data  are  presented 
in  Appendix  VIU. 


Comparisons  of  peak,  amplitude,  and  level  crossings  count  are  made  in 
Section  11. 
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HOTE:  This  i*  only  an  Ill  oat  rat  Ion  and  should  not  be  considered  aa  typical  of 
turbulence  data. 


Figure  9.1  ^Bftlltude  Count  Bxacyle 
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Figure  9.5  Amplitude  Count  Ciamlativa  Probability  Distribution* 


10.  LEVEL  CROSSINGS  COUNT 


Another  method  of  gust  velocity  amplitude  analysis  was  the  level  crossing 
technique.  Levels  of  gust  velocity  were  established  at  2  fps  intervals 
from  the  zero  mean.  Crossings  vith  a  positive  slope  were  counted.  Figure 
10.1  illustrates  the  level  crossings  technique.  This  figure  shove  hov 
crossings  of  corresponding  positive  and  negative  levels  were  combined. 
Crossings  of  the  zero  level  were  a  special  case  and  were  doubled  to  account 
for  crossings  with  both  positive  and  negative  slopes.  The  number  of  level 
crossings  per  mile,  probability  of  exceedance,  and  probability  density  dis¬ 
tributions  were  calculated  as  shown  in  Figure  10.1.  Detailed  discussions 
of  the  level  crossings  analysis  technique  can  be  found  in  References  10,1 
and  10.2. 

A  level  crossings  standard  deviation  (#  J  was  calculated  using  the  follow¬ 
ing  equation: 


(10. l) 


According  to  Rice  (Reference  3.1),  level  crossings  count  provides  a  de¬ 
scription  of  a  stationary  Gaussian  process.  The  level  crossings  distri¬ 
butions  obtained  for  each  turbulence  sample  were  compared  with  a  distri¬ 
bution  obtained  from  Equation  10.2  (Rice's  equation). 

N(x)  =  Nol  exp  (-xa2 /2#tZ)  (10.2) 

where: 

N(x)  a  Expected  number  of  crossings  per  mile  of  a  given  f> 
positive  and  negative  level  ^ 

Nql  =  Total  number  of  zero  crossings  with  positive  and 
negative  slopes  per  mile 

xa  =  Level  of  gust  velocity  in  ft. /sec. 

»t  a  Standard  deviation  of  time  series  (Equation  9.1) 

Figure  10.2  is  an  example  level  crossing  distribution  plot.  The  curves 
shown  on  these  plots  represent  a  distribution  obtained  from  Rice's  equa¬ 
tion.  Agreement  between  L0-IDCAT  data  and  Rice's  equation  was  found  to 
be  best  for  the  longitudinal  component.  Lice's  equation  underestimated 
nearly  every  LO-IOCAT  sample  measured.  For  a  further  discussion  of  Rice's 
equation,  see  Section  51. 

Figures  10.3  and  10.4  are  level  crossings  probability  distribution  and 
probability  density  distribution  plots,  respectively,  of  the  same  turbu¬ 
lence  sample  as  shown  in  Figure  10.2.  The  probability  distribution,  Fig¬ 
ure  10.3,  shows  a  higher  probability  of  occurrence  of  the  wmal]  rr  gust 
velocity  values  than  the  probability  distributions  obtained  from  either 
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the  peak  count  data  (Figure  8.3)  or  the  amplitude  count  data  (Figure  9*3)* 
This  is  also  shown  in  Figure  10.4  where  the  probability  density  is  lower 
at  one  than  it  is  at  nine  feet  squared  per  second  squared. 

Low  intensity  turbulence  samples  were  accounted  for  In  the  same  manner  as 
described  in  Section  8.  Also,  samples  were  categorized  using  the  same 
technique  as  discussed  in  Section  8. 

Comparison  of  the  peak,  amplitude,  and  level  crossings  count  are  given  in 
Section  11.  The  distributions  and  probability  distributions  of  the  level 
crossings  count  data  for  those  categories  shown  in  Table  8.2  are  plotted 
in  Appendix  VH  and  listed  in  Appendix  VIII. 


Oust 

Level 

(Ups) 

Number  of 
Crossings 

Level 

H2222IB 

EH35W3 

Rusher  of 
Crossings 
per  Milo* 

Probability 

of 

Exceedance 

Probability 

Density 

8 

0 

0 

8 

1.0000 

.1250 

6 

1 

2 

6 

30 

.7500 

.1250 

m 

2 

n 

4 

20 

.5000 

.0625 

2 

3 

6 

3 

15 

.3750 

.1875 

0 

4 

-2 

3 

-4 

-  2 

-6 

2 

-8 

0 

*  For  this  example,  assume  0.2  statute  miles  is  length  of  above  record. 

HOTS:  This  is  only  an  example  and  should  not  be  considered  as  typical  of 
turbulence  data. 


Figure  10.1  Level  Crossings  Court  Example 
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Figure  10-2  Oort  Velocity  Level  Crossing  distribution 
for  an  Individual  Turbulence  Sample 
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11.  COHPARISOH  or  COUMTIBO  CTCHHiqUKS 


The  peak,  aaplitude,  and  level  crossings  court  techniques  used  for  analyz¬ 
ing  gust  velocity  tine  aeries  data  were  compared,  lata  for  an  individual 
turbulence  easpla  and  the  all  data  category  are  shown  here  to  Illustrate 
the  relationship  between  counting  techniques,  nils  relationship  la  slnUar 
for  all  categories  so  they  will  not  be  shown.  Plots  and  listings  for  all 
the  categories  Investigated  are  available  In  Appendix  VH  and  VUI,  respec¬ 
tively. 

Occurrences  per  alia  of  peaks  and  level  croaaings  for  a  single  turbulence 
sample  are  coapaxed  In  figures  11.1  through  11.3*  These  plots  Show  that 
generally  a  greater  frequency  of  occurrence  of  level  croaaings  than  of 
peaks  above  a  given  gust  velocity  can  be  expected  for  all  three  coaponents 
of  gust  velocity,  figures  11.4  through  II. o  provide  a  comparison  of  peaks 
atwl  level  crossings  per  idle  for  the  all  data  category.  These  figures  show 
a  relationship  very  sis liar  to  that  shown  for  a  single  turtulence  aanple. 

The  emulative  mater  of  peaks  In  the  first  hand  should  agree  with  the  am¬ 
ber  of  level  crossings  counted  In  the  first  band.  This  is  because  the  peak 
count  technique  counts  tin  positive  and  negative  peak  excursions  following  a 
crossing  of  the  zero  assn  and  the  level  crossings  technique  doubles  the  tim¬ 
ber  of  crossings  of  the  zero  neon  with  a  positive  slope.  Therefore,  It  should 
be  noted  that  the  vnoxtnpalated  peak  count  emulative  occurrences  per  nlle 
equal  the  level  crossings  occurrences  per  adle  at  zero  gust  velocity.  This 
say  be  seen  by  eosparlng  the  uoextrapolated  peak  count  values  of  figure  8.2 
with  the  level  crossings  count  values  of  figure  10.2  at  the  zero  gust  ve¬ 
locity  level. 

Probabilities  of  encotaxterlng  a  gust  velocity  equal  to  or  greater  than  a 
given  v?lue  ftor  the  three  counting  techniques  are  shown  In  Figures  U«  7 
through  11.12.  figures  11.7  through  11.9  show  these  comparisons  for  a  sin¬ 
gle  turbulence  sample.  Figures  11.10  through  11.12  show  the  comparison  for 
tbs  all  data  category. 

The  relationship  between  the  counting  techniques  Is  very  similar  for  both 
the  single  — t1*  and  the  all  data  category.  Also,  there  is  little  varia¬ 
tion  firm  one  component  of  gust  velocity  to  another.  In  general,  the  ampli¬ 
tude  count  distributions  axe  lever  end  the  level  crossings  are  highest.  The 
peak  count  probability  curve  typcraxlaates  the  amplitude  count  curve  at  the 
lower  gust  velocity  values  sod  approaches  the  level  crossings  curve  at  the 
higher  gust  velocities. 

Probability  density  distributions,  obtained  froa  a  single  turbulence  sample 
vising  tbs  three  counting  techniques,  are  shown  in  figures  11.13  through  11.19 
for  the  three  goat  velocity  components,  figures  11.1b  through  U.lB  show 
the  asms  infoersetten  far  the  all  data  category.  Tbs  amplitude  count  proba¬ 
bility  density  distribution  curve  la  lower  than  tar  peak  oo«t  for  the 
■<»gt«  turbulence  saqple,  os  shorn  In  figures  11.13  through  11.15.  The 
level  rmsetrgs  d  nelty  carve  has  a  hi  ay  near  zero  gust  velocity  and  a 
slope  scaawha*  steeper  than  the  others.  This  is  representative  of  the 
level  crossings  probability  density  curve  for  an  Individual  sampl*  with 
a  aerial  gust  velocity  greater  than  mpproxlaetsly  10  fpe.  tar  a  sample 


with  a  maximal  gust  less  than  approximately  H)  fps,  the  density  curve  Is 
shaped  very  much  like  the  peak  and  amplitude  count  curves.  The  reason  for 
the  difference  In  curve  shape  Is  due  to  the  Increased  number  of  zero  cross- 
logs  for  those  turbulence  sasples  of  low  Intensity.  Figures  11.  l£  through 
11.18  show  that  generally  level  crossings  have  the  highest  probability 
density  curve  and  amplitude  court  the  lowest  with  peak  count  in  between. 

The  differences  between  the  three  curves  are  small  and  their  shape  is  gen¬ 
erally  the  Bam.  There  is  no  hup  in  the  level  crossings  curve  as  there 
was  for  seem  Individual  sample*.  This  is  due  to  the  combining  of  many 
turbulence  samples  of  both  high  and  low  intensity  into  this  category. 

Gust  velocity  standard  deviations  from  solitude  count,  peak  count,  and 
level  crossings  were  calculated  as  shown  by  Equations  8.5,  9«1>  and  10.1. 
Figures  11.19  through  11.21  show  comparisons  of  the  cumulative  distribu¬ 
tions  of  these  standard  deviations  for  the  three  gust  velocity  components. 
All  L0-L0CAT  Phase  IH  data  are  included. 

The  level  crossings  standard  deviations  show  a  higher  probability  of  oc¬ 
currence  than  for  peak  count  or  amplitude  count.  This  relationship  1s 
similar  for  all  three  components  of  gust  velocity,  except  the  peak  count 
distribution  is  higher  relatively  for  the  vertical  conponent. 

It  should  be  noted  that  the  characteristic  frequency,  S0,  as  referred  to 
in  this  report  is  not  the  actual  characteristic  frequency  of  the  turbulence 
but  is  equivalent  to  the  intercept  of  the  peak  count  cimu1  ative  distribu¬ 
tion  function.  Measurement  of  the  actual  characteristic  frequency  of  the 
turbulence  would  require  calculation  out  to  very  high  frequencies  (around 
30,000  cps). 

Values  of  the  characteristic  frequency  were  obtained  by  tnese  different 
methods.  In  the  first  method,  H0  (in  cycles  per  foot)  was  obtained  from 
PSD  data  using  Equation  11.1. 

Nq  /  *  cycles  per  foot  (ll.l) 

H0  was  converted  to  cycles  par  mile  and  multiplied  by  two  to  obtain  oc¬ 
currences  per  mile. 

In  the  second  method,  Np  (peaks  per  mile)  was  determined  by  extrapolation 
of  the  peak  court  distribution  curve  as  previously  described  in  Section  8. 

A  third  characteristic  frequency,  N0L  (crossings  per  mile),  was  calculated 
from  level  crossings  data  using  Equation  11.2. 

"ol  -  W-55  Js  (u.2) 

where: 

tm  m  Total  number  of  crossings  of  zero  level  with  both  positive 
and  negative  elope. 

V  ■  Average  speed  during  turbulence  sample. 
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Average  values  for  each  method  axe  given  in  Table  11.1. 

TABLE  11.1 

.V 

COMPARISON  OF  CHARACTERISTIC  FREQUENCIES  CALCULATED  FROM 

POWER  SPECTRA,  PEAK  COUNT,  AND  LEVEL  CROSSINGS  DATA 
# 

Characteristic  Frequenc 
Nq  (occurrences/mile) 

Np  (peaks/mile) 

Nol  (crossings/mlle) 


These  averages  represent  Phase  III  data.  Power  spectra  and  N0  values  were 
calculated  for  approximately  50  per  cent  of  the  data  processed.  In  addi¬ 
tion,  characteristic  frequency  values  were  not  calculated  for  Irregularly 
shaped  spectra  (l.e.  nonhcoogeneous  turbulence  or  turbulence  samples  having 
low  signal -to-noise  ratios).  Np  and  Nop  values  were  calculated  for  all  tur¬ 
bulence  samples  processed.  Extraneous  Np  and  N0l  values  defined  as  those 
which  deviated  more  than  four  standard  deviations  from  the  mean  were  removed 
before  averaging. 

The  characteristic  frequency  of  the  longitudinal  component  is  the  smallest 
and  vertical  is  the  largest.  This  Is  true  for  both  N0  and  Np.  Comparison 
of  corresponding  N0  and  Np  values  shows  reasonable  agreement  between  the 
two  methods  of  computation,  except  for  vertical  where  N0  is  significantly 
less  than  Np.  This  is  believed  to  be  due  to  the  removal  from  analysis  of 
a  higher  proportion  of  the  N0  than  Np  data.  These  N0  values  were  removed 
primarily  because  of  distortion  of  the  spectra  at  higher  frequencies. 

values  of  characteristic  frequency  obtained  using  level  crossings  data 
considerably  larger  than  those  obtained  from  peak  count.  This  is  be¬ 
cause  the  level  crossings  values  vorrespond  to  what  the  peak  count  values 
would  have  been  had  they  not  been  extrapolated.  The  reason  for  extrapolat¬ 
ing  to  obtain  the  characteristic  frequencies  from  peak  count  data  is  ex¬ 
plained  in  Section  8, 

It  should  be  noted  that  the  characteristic  frequencies  obtained  using  the 
peak  count  and  level  crossings  method  are  dependent  upon  the  frequency 
range  of  turbulence  being  investigated.  During  L0-L0CAT  Phases  I  and  II, 
the  C-13IB.  test  airplanes  were  flown  at  an  average  speed  of  approximately 
330  feet  per  second  while  gathering  gust  velocity  data.  TSrese  gust  velocity 
data  were  filtered  so  that  only  the  frequency  band  from  approximately  U.04 
to  10  cps  was  .retained.  This  frequency  band  corresponds  to  a  range  of  gust 
velocity  wavelengths  ranging  from  79^0  to  33  feet.  During  LO-LOCAT  Phase 
rZZ,  the  T-33A  test  airplane  was  flown  at  an  Average  speed  of  approximately 
630  feet  per  second.  The  same  band  of  frequencies  was  passed  during  data 
processing  as  for  Phases  I  and  II.  With  the  higher  airplane  speed  and  the 
Same  filtering  bands,  the  new  range  of  gust  velocity  wavelengths  obtained 
ranged  from  25, 100  to  63  feet.  The  segment  of  the  turbulence  spectrum 
investigated  during  LD-L/OCAT  Phase  III  wa#  therefore  shifted  toward  longer 
wavelengths  from  that  of  Phases  I  and  II . 
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The  average  values  of  Np  obtained  during  Phases  I  and  II  were  50  peaks 
per  mile  for  the  longitudinal  components  and  60  for  lateral  and  vertical 
components.  Comparison  of  these  values  with  those  obtained  during  Phase 
III  shows  the  difference  between  them  to  be  a  factor  of  approximately  2. 
This  difference  Is  due  to  the  high  frequency  of  occurence  per  mile  of 
the  33-63  foot  wavelength  gusts  which  were  measured  during  Phases  I  and  II. 
Wavelengths  less  than  60  feet  were  not  measured  during  Phase  IH. 
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.  CCSRELAT3DH  OF  PEAK  SISaSaBUTIGHS  WITH  CBSCFHXBICAL  GAZEQGHY 

1  111  •  1  . .  . . . i  . . — »< 

/ 

Peak,  amplitude,  and  lavel  crossings  court  distribution*  axfe  very  similar 
in  the  way  they  vary  with  different  geophysical  phenomena.  Therefore,  only 
the  peak  count  distributions  are  used"  In  this  dociamnt'to  illustrate  these 
variations.  Plots  and  listings  of  data  from  all  three  counting  techniques 
are  presented  In  Appendices  VH  and  VHI. 

Figures  12.1  through  12.3  dhow  the  terrain  effects  on  peak  count  distribu¬ 
tions  of  the  three  gust  velocity  components.  These  plots  show  that  the  num¬ 
ber  of  peaks  per  exceeding  given  values  are,  in  general,  larger  for 
high  mountain  data  and  decrease  for  low  mountains,  plains,  desert,  and  water 
in  that  order.  f 

Peak  probabilities  of  the  plains  data  are  higher  than  that  of  the  desert 
data,  even  though  the  average  desert  terrain  was  rougher  than  that  of  the 
plains.  Tbifl  was  probably  caused  by  the  fact  that  higher  wind  velocities' 
were  recorded  over  the  plains  resulting  in  a  higher  percentage  or  mechan¬ 
ically  generated  turbulence. 

/ 

A  distinct  change  was  seen  in  the  distribution  curves  for  terrain  effects 
between  the  Phase  m  Interim  Report  (Reference  12.1)  and  this  report. 

There  was  no  change  in  water  and  desert  distribution  curves  since  no  addi¬ 
tional  data  wer^  obtained  for  these  categories.  The  plains,  low  mountain, 
and  high  mountain  distribution  curves  were  raised  considerably.  Ibis  was 
due  primarily  to  the  effects  of’the  more  turbulent  air  encountered  over  the 
Peterson  route.  These  effects  were  not  included  in  Reference  12.1. 

Figures  12.4  through  12.6  presort  the  altitude  effects  on  peak  count  dis¬ 
tributions  .  These  figures  show  greater  frequencies  of  occurrence  of  gust 
velocity  peaks  at  250  feet  than  at  750  feet. 

* 

"iguros  12.7  through  12.9  show  the  effect  of  atmospheric  stability  on  peak 
covert  distributions.  In  general,  very  stable  air  has  the  lowest  distribu¬ 
tion,  whereas,  the  distributions  of  data  obtained  in  stable,  neutral,  and 
unstable  air  are  nearly  the  same.  This  was  not  the  case  for  the  Phase  III 
Interim  Report  data  where  the  difference  between  distributions  for  the  very 
stable  and  other  stabilities  was  not  so  distinct.  The  distribution  curves 
for  all  Tour  stabilities  were  raised  above  what  they  were  for  the  interim 
analysis.  But,  the  curves  for  stable,  neutral,  and  unstable  were  raised 
considerably  more  than  the  very  stable  curve.  This  w4a  probably  due  to  the 
fact  that  a  such  smaller  percentage  of  the  very  stable  data,  as  compared  to 
the  other  three  stabilities,  came  from  the  Petersen  route. 

Figures  12.10  throvgh  12.12  show  the  effect  of  time  of  day  on  peak  count 
distributions .  These  plots  show  that  the  distributions  do  not  vary  signif¬ 
icantly  with  time  of  day.  ihe  only  distinguishable  characteristic  Is  that 
the  distribution  fbr  dawn  starts  out  below  the  others  at  —ail  gust  veloc¬ 
ities,  but  crosses  over  at  about  30  fpa  and  generally  la  higher  from  there 
on.  This  is  probably  due  to  the  effects  of  the  high  percentage  of  the  low 
intensity  turbulence  twylss  that  ware  included  in  the  dawn  category.  Forty 
two  of  the  66  low  intensity  turbulence  eaaples  (Table  8.2)  were  obtained  at 
dawn.  •  Also,  the  effects  of  very  stable  samples  should  be  considered,  since 


68  per  cent  of  the  dawn  data  was  very  stable.  The  curve  for  dawn  data  was 
lower  at  the  smaller  gust  velocities,  however,  due  to  the  effects  of  the 
higher  intensity  turbulence  samples  obtained  at  Peterson,  the  curve  is 
higher  at  the  larger  gust  velocities. 

Figures  12.13  through  12.15  show  the  variation  of  peak  count  distributions 
with  geographical  location.  The  distribution  for  the  Peterson  route  was 
the  highest.  This  is  due  to  the  high  mountain  terrain  and  the  season  of 
the  year  during  which  this  route  was  flown.  The  Edwards  distribution  is 
the  lowest  at  the  smaller  gust  velocities  due  to  the  effects  of  the  water 
and  desert  legs,  but  increases  rapidly  with  increasing  velocity  due  to  the 
effects  of  the  mountain  legs.  The  distribution  for  the  Griff iss  route  is 
somewhat  higher  than  for  the  McConnell  route  due  to  the  effects  of  low 
mountain  legs  at  Griff iss. 

It  will  be  noted  that  the  lateral  component  for  the  Peterson  route  exhibits 
e  higher  distribution  than  that  of  the  vertical  component,  while,  for  the 
Edwards  route,  the  vertical  component  is  higher  than  the  lateral.  This  is 
due  to  location  of  those  legs  with  respect  to  the  more  rugged  terrain  char¬ 
acteristics.  The  legs  at  Peterson  paralleled  the  ridge  lines  while  at  Edwards 
the  legs  crossed  over  the  ridge  lines.  Thus,  higher  lateral  than  vertical 
gust  velocities  would  be  experienced  over  the  Peterson  route  and  higher  ver¬ 
tical  than  lateral  gust  velocities  would  be  experienced  over  the  Edwards 
route. 

The  variation  of  peak  count  distributions  with  season  was  not  determined 
during  Phase  III,  since  the  airplane  was  not  flown  over  each  route  for  a 
full  year. 


IVE  PEAKS  PER  ST 


Cumulative 
Miles  Represented 


G  High  Mountain  27469.0 
□  Low  Mountain  9384.0 
O  Desert  3569*0 
A  Plains  12399*0 
^  Water  2234.0 


Altitude  -  250  and  750  Feet 

Stability  -  All 

Time  of  Day  -  All 

Season  -  All 

Location  -  All 


XVE  PEAKS  PER 


91 

$ 

I 

: 

€ 


n,o 


FEAKS  PER 


1 


Cumulative 
Miles  Represented 


0 

0 


Dawn 

Mid-Morning 

Mid-Afternoon 


Terrain  -  All 

Altitude  -  250  and  750  Feet 
Stability  -  All 
Season  -  All 
Location  -  All 


15019.0 

20718.0 

20016.0 


20  to  60  80 

OUST  VELOCITY  AMPLITUDE  •  fPt 

Figure  12 .10  Tims  of  Day  Effects  on  Longitudinal  Gust  Velocity  Peeks 

98 


100 


0 


Cumulative 
Miles  Represented 


©  Edvards  Air  Force  Base,  California  18558-0 
□  Griffis s  Air  Force  Base,  Nev  York  4284.0 
O  Peterson  Field,  Colorado  2453-1.0 
A  McConnell  Air  Force  Base,  Kansas  7703-0 


Terrain  -  All 

Altitude  -  250  and  750  Feet 
Stability  -  All 
Tine  of  Day  -  All 
Season  -  All 


20  kQ  60  80 

GUST  VELOCITY  AMPLITUDE  -  FPS 

Figure  12.15 Location  Effects  ou  Vertical  Gust  Velocity  Peaks 

103 


SECTION  IV 


GUST  VEI  WVT  INTENSITY 


The  time  series  standard  deviation  (*t)  of  each  of  the  three  gust  velocity 
components  vas  calculated  by  the  following  equation: 


r  »  *i1/2 


=  rms  (with  x  =  0)  (9*l) 


The  values  for  x^  were  taken  from  the  time  histories  of  gust  velocities  which 
were  computed  at  the  rate  of  100  samples  per  second  giving  27,000  values 
for  each  four"  and  one-half  minute  turbulence  sample .  Therefore,  for  this  *' 
particular  calculation,  N  =  27,000.  The  mean  gust  velocity  value  is  indi¬ 
cated  by  x. 

Since  the  mean  value,  x,  is  made  essentially  equal  to  zero  by  high-pass 
filtering,  the  standard  deviation  is  actually  equal  to  the  root -mean-square 
(rms)  value  as  shown.  The  two  terms  are  used  interchangeably  in  this  report. 


13.  GUST  VEIPCITY  RMS  gI]ATISTICAL  ANALYSIS 


The  rms  values  were  grouped  into  amplitude  bands  of  0.5  fps  and  the  cumu¬ 
lative  probability  distribution  (CPD)  was  calculated.  These  CPD  values 
give  the  probabilities  of  a  gust  velocity  rms  value  being  equal  to  or 
greater  than  a  given  value.  The  cumulative  probabilities  determined  for 
all  data  recorded  during  the  Phase  III  program  are  shown  in  Figure  13.1. 

The  cumulative  probability  distribution  curve  for  lateral  rms  gust  velocity 
i6  higher  than  for  longitudinal  since  a  high  percentage  of  the  turbulence 
data  were  obtained  while  flying  parallel  to  mountain  ridge  lines.  The  curve 
for  vertical  rms  gust  velocity  was  lower  than  for  longitudinal  due  to  the 
compression  of  the  vertical  eddiv  r-  near  the  ground. 

For  categories  which  did  not  contain  any  high  mountain  data  the  cumulative 
probability  distribution  curves  for  longitudinal  and  lateral  rms  gust  veloc¬ 
ity  were  in  general  nearly  the  same. 

The  CPD  of  rms  values  for  various  geophysical  categories  was  calculated. 

In  order  to  assure  a  good  statistical  sampling,  a  minimum  of  30  data  sam¬ 
ples  per  category  vas  established  as  the  criteria  necessary  for  a  valid 
analysis. 

The  gust  velocity  rms  distributions  are  shown  as  functions  of  terrain  type, 
absolute  altitude,  atmospheric  stability,  time  of  day,  and  geographic  loca¬ 
tion  in  Figures  13.2  through  13.6.  As  can  be  seen,  the  results  are  similar 
for  all  three  gust  velocity  components. 
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The  highest  gust  velocity  rms  values  were  recorded  over  the  high  mountain 
legs  and  the  lowest  over  the  water  (Figure  13.2).  In  general,  the  gust 
velocity  rms  cumulative  probabilities  are  larger  at  a  given  rms  value  for 
data  recorded  over  the  rougher  terrain.  The  only  exception  occurs  for 
data  recorded  over  the  plains  and  desert.  Larger  gust  velocity  rms  values 
were  recorded  over  the  plains,  even  though  the  desert  terrain  was  some¬ 
what  rougher  (see  Section  32).  This  is  attributed  to  the  fact  that  the 
winds  were  significantly  greater  over  the  plains  (Section  33). 

The  data  were  recorded  at  the  two  nominal  altitudes  of  250  and  750  feet 
above  the  terrain.  Resultant  gust  velocity  rms  probabilities  are  shown 
in  Figure  13. 3,  and,  as  can  be  seen,  the  rms  values  are  greater  at  the 
lower  altitude. 

The  effect  of  atmospheric  stability  on  gust  velocity  rms  values  is  shown 
in  Figure  13. The  lowest  rms  values  were  recorded  when  the  atmospheric 
stability  was  classified  as  very  stable.  The  cumulative  probability  in¬ 
creases  with  decreasing  atmospheric  stability,  although  the  differences 
between  these  values  for  the  stability  classifications  of  stable,  neutral, 
and  unstable  are  very  small. 

As  the  time  of  day  during  which  the  data  were  recorded  varies  from  dawn 
to  mid-morning,  to  mid-afternoon,  the  cumulative  probability  of  encounter¬ 
ing  a  gust  rms  of  a  given  magnitude  increases.  This  is  shown  in  Figure 

13.5. 

Gust  velocity  rms  CPDs  are  shown  as  a  function  of  the  four  geographic  loca¬ 
tions.  The  greatest  gust  velocity  rms  values  are  associated  with  the 
Peterson  and  Griff iss  routes,  as  shown  in  Figure  13.6.  As  mentioned  pre¬ 
viously,  Peterson  was  the  only  location  where  data  were  obtained  over  more 
than  one  season.  These  winter  and  spring  data  are  shown  in  Figure  13.7* 

The  cumulative  probability  is  greater  for  a  given  gust  velocity  rms  value 
during  the  winter  season. 

Based  on  tne  findings  of  Phases  I  ana  II  and  on  the  data  in  Figures  13.2 
through  13.6,  it  appears  that  combinations  of  terrain,  altitude,  and  sta¬ 
bility  are  primary  factors  in  the  determination  of  gust  velocity  rms  cumu¬ 
lative  probability  distributions.  The  distributions  for  the  combinations 
of  the  three  category  components  are  shown  in  Figures  13.8  through  13.27. 

In  each  of  these  figures,  only  one  of  these  category  components  is  varied, 
while  the  other  two  are  held  constant.  Only  data  obtained  over  high  moun¬ 
tains  and  plains  are  analyzed  in  this  meaner.  Combinations  of  these  three 
category  components  involving  data  obtained  over  the  other  types  of  terrain 
did  not  generally  contain  a  sufficient  number  of  turbulence  samples.  Like¬ 
wise,  the  data  obtained  over  the  plains,  at  750  feet,  during  the  condition  of 
very  Stable  atmospheric  stability  did  not  contain  the  necessary  30  samples 
and  is  not  shown  in  Figures  13.8  through  13.27.  However,  these  figures  show 
the  effects  of  variations  in  the  geophysical  parameters  for  rough  and  smooth 
types  of  terrain. 

The  data  in  Figures  13.8  through  13.27  indicate  approximately  the  same  re¬ 
sults  as  were  indicated  In  Figures  13.2  through  13.6.  Gust  velocity  rms 
values  are  greater  over  the  rougher  terrain,  at  the  lower  altitude,  and  for 
atmospheric  stability  conditions  other  than  very  stable. 
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The  data  obtained  over  the  lew  ncuntairs,  desert,  and  rater  were  not  suffi¬ 
cient  to  permit  an  extensive  investigation  of  geophysical  effects.  However,, 
these  data  were  analyzed  as  a  function  of  altitude  and  type  of  terrain  (Fig- 
vires  13.28  through  13.33).  The  high  mountain  and  plains  data  are  also  shown 
for  comparative  purposes.  In  general,  it  ras  found  that  for  the  various  com¬ 
binations  of  terrain,  altitude,  and  atmospheric  stability,  gust  velocity  rms 
cumulative  probability  becomes  greater,  for  a  given  rms  value,  as  terrain 
becomes  rougher,  altitude  decreases,  and  atmospheric  stability  becomes  more 
unstable. 


It  has  been  found  that  the  effects  of  noise  are  significant  at  the  smaller 
gust  velocity  rms  magnitudes.  According  to  theory  for  normal  distributions: 


« 


2 

R 


•;  *  *; 


(13.1) 


where: 


o’  t  =  standard  deviation  of  the  turbulence 

=  standard  deviation  of  the  noise 

=  standard  deviation  of  the  recorded  data 

If  is  assumed  constant  for  a  given  category,  this  value  becomes  a  greater 
and  greater  percentage  of  oR  as  #t  becomes  smaller  and  smaller.  Therefore, 
at  lower  turbulence  intensities  the  distribution  is  distorted  by  the  high 
percentage  of  noise  content  in  the  data  where  the  noise  has  a  greater  effect 
on  the  distribution  than  does  the  turbulence. 

During  the  LO-LOCAT  Phase  III  Program,  gust  velocities  were  not  computed 
for  turbulence  sauqcles  when  the  level  of  turbulence  and  hence,  the  r.ignal- 
to-noise  ratio,  was  known  to  be  small.  The  criteria  used  to  determine  the 
level  of  turbulence  were  the  values  of  the  angles  of  attack  and  sideslip 
differentia.;,  pressures  as  sensed  by  the  gust  probe  pressure  ports.  Gust 
velocities  rare  calculated  only  for  turbulence  samples  where  the  standard 
deviation  of  both  the  angle  of  attack  differential  pressure  (<ra)  and  angle 
of  sideslip  differential  pressure  (9  J  were  greater  than  0.07  and  0.05 
inches  of  Hg,  respectively.  This  Is  discussed  further  in  Appendix  III, 
Section  III. 5. 

Because  of  the  low  signai-to-noise  ratio  of  some  turbulence  samples,  the 
cumulative  probability  distributions  have  not  been  plotted  below  a  gust 
velocity  rms  value  of  1.5  fps,  since  this  is  the  lowest  band  limit  which 
is  greater  than  the  minimum  valid  gust  velocity  rms  value  of  1.4  fps.  The 
nuntoer  of  turbulence  samples  having  rms  values  below  1.5  feet  per  second 
were  used  to  compute  the  probabilities  for  the  data  above  3,5  fps ,  The 
number  of  turbulence  samples  not  processed  due  to  the  and  *  $  criteria 
were  assumed  to  have  gust  velocity  rms  values  below  1. 5  fpe  and  were  in¬ 
cluded  in  the  cumulative  probability  distribution  calculations  by  being 
added  into  the  total  number  of  turbulence  samples  within  the  category 
being  investigated.  For  example,  *  ....  ■  the  following  distribution 
which  excludes  those  samples  not  ne  9  a  and  9  s  criteri*  but  in¬ 
cludes  those  having  a  less  than  and  meeting  the  »a  and  9^  criUri 
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Suppose  that,  because  of  the  #„  and  criteria,  50  samples  were  not  pro 
cessed.  Then  the  following  tabulation  illustrates  the  method  used  to  ac 
count  for  these  50  sauries  when  calculating  the  gust  velocity  rms  proba¬ 
bility  distributions . 
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It  should  be  pointed  out  that  this  procedure  was  only  necessary  for  cate¬ 
gories  which  included  data  obtained  over  plain3,  desert,  and  water  legs. 

It  was  over  these  types  of  terrain  where  the  data  having  low  signal-to- 
noise  ratios  was  gathered.  Categories  consisting  exclusively  of  data  ob¬ 
tained  over  high  or  low  mountains  were  not  involved.  The  number  of  low 
curb 'ilence  samples  within  the  categories  pertaining  to  the  peak  count, 
amplitude  count,  and  level  crossing  data  analysis  are  shown  In  Tables 
8.1  and  8.2.  The  number  of  samples  included  In  those  tables  is  not  entirely 
applicable  to  the  gust  velocity  ras  value  analysis  because  of  the  following 
reasons . 

•  Peak  count,  amplitude  count,  and  level  crossings  data  were 
put  on  a  master  tape,  as  explained  in  Section  III.  The  pro¬ 
gram  by  which  this  was  accomplished  was  designed  to  put  t.ue 
data  on  tape  in  the  order  of  the  six-digit  category.  The 
atmospheric  stability  category  component  was  not  defined  for 
some  of  these  low  turbulence  samples  because  of  instrumenta¬ 
tion  malfunctions.  Therefore,  these  particular  samples  could 
not  be  used  in  the  peak  count,  amplitude  count,  and  level 
crossings  data.  The  gust  velocity  rms  values  were  not  con¬ 
strained  by  this  six -digit  category  number.  Thus,  even  though 
stability  was  not  defined  for  some  low  turbulence  samples, 
tnese  samples  were  st 511  included  in  category  combinations 
including  components  other  than  stability. 

•  When  the  gust  velocity  nn3  values  were  being  individually  re¬ 
viewed  it  was  noticed  that  a  number  of  values  obtained  over 
the  water  leg  at  Edwards  were  loss  than  1.0  fps.  These  data, 
questionable  because  of  their  low  signal -to-noise  ratios,  were 
treated  as  low  intensity  turbulence  samples  in  the  gust  veloc¬ 
ity  nns  analysis. 

Table  13.1  shows  the  number  of  low  intensity  turbulence  samples  and  the 
corresponding  total  number  of  samples  included  in  tne  category  combina¬ 
tions  utilized  in  the  gust  velocity  ms  analysis.  A  tabulation  of  all 
valid  values  Is  presented  in  Appendix  VI.  The  number  of  low  intensity 
turbulence  samples  consists  of  those  not  processed  due  to  ra  and  cri¬ 
teria  and  those  low  intensity  samples  recorded  over  the  water  legs  as 
mentioned  above. 

According  to  the  theory  of  isotropic  turbulence,  #tl  =  #  -  #tw.  Refer¬ 

ence  13.1  suggests  that  this  relationship  does  not  hold  true  for  turbulence 
up  to  about  1000  feet  above  the  terrain.  This  reference  proposes  that  the 
relationship  between  the  rma  values  of  the  three  components  is  a  function 
of  atmospheric  stability  and  absolute  altitude. 

For  very  stable,  stable,  and  neutral  atmospheric  conditions: 

0  0t 

— ^  «  — -  =  1.2  -  0.00017  a  (13.2) 
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TABLE  13.1 


TOTAL  NUMBER  OF  SAMPLES  AND  NUMBER  OF  LOW  INTENSITY 
'TURBULENCE  SAMPLE, S  INCLUDED  IN  GUST  VELOCITY  RMS  ANALYSIS 


Category 

Total 

Number  of  Sam] 

Dies 

Number  of 

Low  Intensity 
Turbulence  Sandies 

u 

V 

w 

OOOOOO 

1762 

1740 

1746 

113 

100000 

653 

851 

853 

0 

200000 

297 

297 

297 

0 

300000 

115 

115 

115 

1 

400000 

415 

395 

401 

50 

600000 

83 

83 

31 

6-2 

010000 

909 

894 

900 

37 

020000 

853 

846 

846 

76 

001000 

570 

564 

568 

35 

002000 

4S0 

477 

478 

7 

003000 

394 

389 

391 

0 

004000 

296 

28o 

288 

0 

000100 

491 

485 

490 

67 

000200 

648 
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For  unstable  atmospheric  conditions 


~  -  -r1  =  1.3  -  0.00058  h 

*tr  *tv 


(13.3) 


Values  from  Equations  13.2  and  13.3  are  compared  to  LO-LOCAT  Phase  III 
data  in  Table  13.2. 


TABID  13.2 


LO-LOCAT  GUST  VELOCITY  RMS  RATIOS  COMPARED  TO  T3DSE 
_ BBCQjtjEKEED  EY  REFERENCE  13.1 _ 


Absolute 

Values 

Predicted 

LO-LOCAT 

Phase  III  Data 

Altitude 
(H)  ~  Ft. 

Atmospheric 

Stability 

by  Equations 

13.2  and  13.3 

•t« 

•tu 

•tv 

•tv 

250 

Very  Stable, 
Stable,  and 
Neutral 

l.l£ 

1.07 

1.13 

250 

Unstable 

1.15 

1.11 

1.14  ' 

750 

Very  Stable, 
Stable,  and 
Neutral 

1.07  . 

1.05 

1.15 

750 

Unstable 

0.86 

.. 

0.9 9 

1.06 

_ 

The  data  show  approximately  the  same  magnitudes  as  predicted  by  the  equa¬ 
tions.  At  the  750  foot  altitude,  smaller  ratios  occur  for  unstable  than 
for  very  stable,  neutral,  and  unstable  conditions. 

Reference  13.1  also  presents  the  vertical  gust  velocity  rms  cumulative 
probability  distributions  from  several  different  turbulence  research  pro¬ 
grams,  as  veil  as  a  distribution  ^recommended  for  data  at  an  absolute  alti¬ 
tude  of  500  feet.  The  data  from  the  LO-LOCAT  Phase  III  Program  are  compared 
to  those  distributions  in  Figure  13.34. 

The  gust  velocity  rms  cumulative  probability*-  distributions  for  all  three 
components  can  be  closely  represented  by  portions  of  two  normal  distribu¬ 
tions,  i.e.: 


no 


where: 


9  =  dispersion  of  the  distributions 
M  =  mean  of  the  distributions 

and:  subscript  2  refers  to  that  portion  of  the  probability  distribution 
between  rms  values  from  1.5  to  a,  and  subscript  3  refers  to  that  portion 
of  the  distribution  between  rms  values  from  a  to  b. 

The  constants  for  Equation  13.4  are  shown  in  Table  13.3  for  each  component 
The  agreement  between  the  distribution  of  the  recorded  data  and  the  distri 
buttons  represented  by  the  equation  is  shown  in  Figures  13.35,  13*36,  and 
13.37* 


TABLE  13.3 

GUST  VELOCITY  KMS  DISTRIBUTION  CONSTANTS 


Constant 

Gust  Velocity  Conponent  | 

u 

V 

w 

a 

4.79 

4.75 

5.30 

3.35 

3.46 

3.23 

a2 

1.70 

1.65 

1.70 

m3 

1.79 

0.98 

0.65 

ff3 

3.59 

4.79 

3.79 

b 

14.00 

17.00 

14.00 
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Figure  13.1  Gust  Velocity  I34S  Cumulative  Probability 
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Figure  13.2  Gust  Velocity  FWS  Cumulative  Probability  Associated 
with  Type  of  Terrain 
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Figure  13,5  Gust  Velocity  HKS  Cumulative  Probability  Associated 
with  Time  of  Day 
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Figure  13.  l4  Terrain  Effects  on  Gust  Velocity  RMS  Values  -  250  Feet,  Unstable 
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Figure  33-16  Altitude  Effects  on  Gust  Velocity  BMS  Values 
-  High  Mountains,  Very  Stable 
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Figure  13.20  Altitude  Effects  on  Gust  Velocity  SMS  Values  -  Plains,  Very  Stable 
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Figure  13.26 


Atmospheric  Stability  Effects  on  Gust  Velocity  RMS  Values  - 
250  Feet,  Plains 
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Figure  13.28  Terrain  Effects  on  Gust  Velocity  KC  Values  -  250  Feet 
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Figure  13.32  Altitude  Effects  on  Gust  Velocity  RMS  Values  -  Desert 


HMtoliUTT  OM'JLTin  I’MMAlim 


IVE  PROBABILITY 


All  Phase  III  Data 


□  Data 

—  Distribution  Represented  by  Equation  13. ^ 


I  »  1  »  *  *  *  J —  i—l - 1 - 1 - l 

6  8  10  12  14  16  18 

GUST  VELOCITY  RMS  VALUE  -  fps 


13.3^  Lateral  Gust  Velocity  RMS  Uvniulative 
Probability  Distribution  Shape 


l4.  GUST  VELOCITY 


VERSUS  INFLIGHT 


WIND  VELOCITY 


RMS 


MEASURED 


The  relationship  between  gust  velocity  rms  values  and  recorded  inflight 
wind  velocity  was  investigated  with  respect  to  type  of  terrain  ana  alti¬ 
tude  above  the  terrain.  The  terrain  and  wind  involved  during  LO-LOCAT 
Phase  in  are  discussed  in  Section  VI. 

An  analysis  of  rms  values  as  a  function  of  wind  speed  was  accomplished  by 
two  different  methods.  The  first  method  involved  plotting  the  gust  veloc¬ 
ity  rms  values  versus  wind  speed.  The  slopes  and  intercepts  of  first-order 
least  square  lines  (Equation  l4.l)  which  best  fit  the  data  were  computed: 

at  =  y0  +  mW  (i4.i) 

where: 


m  =  slope  of  line 


y0  =  y  intercept 

The  standard  deviation  of  the  data  about  the  least- square  line  was  calcu¬ 
lated  as  follows: 


9 


D 


(14.2) 


where: 

x^  =  value  of  gust  velocity  rms  at  a  given  wind  speed 

Rt  =  value  of  least  square  line  at  the  given  wind  speed 

The  standard  deviation  of  the  data  about  the  least  square  line  indicates 
the  degree  of  data  scatter.  The  results  of  this  method  are  shown  in  Fig¬ 
ures  14.1  through  l4.8.  The  applicable  least  square  line  and  ±  one  stan¬ 
dard  deviation  lines  for  the  data  are  shown  on  each  plot. 

The  slope  of  the  least-square  line  indicates  the  effect  of  wind  speed  on 
gust  velocity  rms  values;  the  greater  the  slope,  the  greater  the  effect. 
The  wind  speed-turbulence  relationships  of  Figures  l4.1  through  14.8  are 
svsnnejrized  in  Figure  l4.<>  by  showing  the  slopes  of  the  least  square  lines 
plotted  as  a  function  of  terrain  and  altitude. 

For  the  second  method  of  analysis,  the  gust  velocity  rms  values  were  sep¬ 
arated  into  several  different  groups.  The  particular  group  within  which 
the  rms  value  for  a  given  data  sample  was  placed  was  dependent  upon  the 
corresponding  wind  velocity  for  the  sample.  Each  group  covered  a  wind 
speed  band  of  ID  fpa.  The  cumulative  probability  of  gust  velocity  rms 
was  then  calculated  for  each  group.  Probability  plots  were  constructed 
and  are  shown  in  Figures  14.10  through  14.13.  Not  all  bands  contained 
enough  data  for  statistical  reliability  (30  samples),  therefore,  all 
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combinations  of  terrain  and  altitude  could  not  be  statistically  analyzed. 

A  greater  change  In  the  vlnd  versus  gust  rms  relationship  vith  terrain 
than  with  altitude  was  Indicated  (Figure  l4.9).  Therefore,  the  data  were 
categorized  on  the  basis  of  terrain  type  only.  Figures  14.10  through 
14.13  show  the  cusulatlve  probability  distributions  of  gust  velocity  zme 
values  as  a  function  of  wind  speed  for  glvm  types  of  terrain.  The  ver¬ 
tical  and  longitudinal  gust  velocity  ras  values  exhibit  approximately  the 
same  relationship  with  respect  to  vlnd  speed  for  both  the  high  and  lov 
mountain  data.  The  lateral  coaycoent  shows  a  greater  relationship  with 
vlnd  speed  over  the  high  mountain  terrain  than  either  the  longitudinal  or 
vertical  components,  while  over  the  lov  mountain  terrain  the.  opposite  is 
true.  This  may  be  due,  In  part,  to  the  vlnd  direction  effects  discussed 
later. 

The  effect  of  altitude  on  the  vlnd  speed  versus  ras  value  relationship 
seems  to  be  small.  For  all  data,  except  those  recorded  over  lov  mountain 
terrain,  gust  velocity  ras  values  shoved  slightly  more  correlation  vith 
wind  speed  at  the  250-foot  altitude. 

Oust  velocity  rms  values  from  legs  within  the  proximity  of  predominant 
mountain  ridges  were  studied  vith  respect  to  wind  direction.  In  order 
to  study  the  effects  of  being  on  the  leeward  or  windward  side  of  a  ridge, 
only  legs  which  were  congletely  on  one  side  or  the  other  were  chosen. 

Also  considered  was  the  wind  direction  with  respect  to  the  ridge  (perpen¬ 
dicular,  diagonal,  or  parallel)  When  the  turbulence  was  being  recorded  an 
the  leeward  side.  The  data  from  legs  5,  6,  7,  and  6  of  the  Peterson  route 
are  the  most  applicable  to  this  type  of  analysis.  The  McConnell  route,  of 
course,  consisted  entirely  of  smooth  terrain.  Although  the  Edvards  and 
Orlfflss  routes  did  Include  rough  terrain  there  were  no  close  predominant 
ridges  completely  on  one  side  of  a  leg. 

The  layout  of  the  Peterson  legs  Is  shown  In  Figure  IV.  3  of  Appendix  IV. 

The  same  ridge  line  was  used  for  a  wind  direction  reference  cm  both  legs 
5  and  6  and  was  as suaed  to  run  In  a  north-south  direction.  The  predomi¬ 
nant  ridge  was  as  aimed  to  run  In  a  340*  direction  for  leg  7*  For  leg  8, 
the  vlnd  direction  was  referenced  to  the  same  ridge  line  as  for  leg  7, 
except  that  effects  of  wind  bloving  from  the  direction  of  the  mountains 
near  Colorado  Springs  were  also  taken  Into  account.  Therefore,  the  ridge 
line  was  assumed  to  curve  and  run  from  340*  to  360*  for  leg  8  data.  The 
data  were  categorized  by  wind  direction,  as  shown  in  Table  i4.1. 

The  plotted  data  are  shown  In  Figures  14.14  through  14.29.  Oust  velocity 
rms  values  were  plotted  versus  vlnd  speed  for  vlnd  direction  category.  In 
Figures  l4,l4  through  14.25,  the  data  for  different  categories  were  plotted 
separately  and  approximate  mean  lines  were  drawn  tJvrough  the  data.  These 
mean  linos  for  each  leg  are  shown  In  Figures  14.26  through  14.29.  Certain 
trends  can  be  seen  although  the  data  are  widely  scattered. 
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The  moot  significant;  differences  between  the  gust  rma  values  recorded  when 
wind  direction  was  perpendicular  to  a  mountain  ridge  and  those  recorded 
for  a  diagonal  vlnd  direction  occur  for  the  leg  ^  data.  The  reason  for 
this  may  be  the  fact  that  the  first  part  of  leg  7  does  not  run  along  any 
predominant  mountain  ridge.  The  airplane  required  only  approximately  32 
statute  miles  to  record  a  4.5  minute  turbulence  sample  although  the  legs 
were  laid  out  to  be  approximately  45  statute  miles  long.  Therefore,  ap¬ 
proximately  one-ball  of  the  data  recorded  on  leg  7  were  recorded  prior  to 
when  the  aircraft  reached  the  beginning  of  the  ridge.  Some  effects  of  the 
ridge  would  probably  be  Included  when  the  wind  direction  was  perpendicular 
to  the  ridge  or  diagonal  from  272.5°  to  317*5°  (see  Table  l4.l).  Less 
effects  should  have  been  present  when  the  wind  blew  diagonally  from  l82. 5° 
to  227*5°.  Only  approximately  8  per  cent  of  the  wind  directions  that  were 
classified  aa  diagonal  were  between  182.5°  and  227.5°.  Thus,  diagonally 
categorized  data  of  leg  7  are  probably  more  related  to  what  could  be  ex¬ 
pected  fncn  parallel  or  windward  wind  direction  data. 

TABLE  14.1 

WIND  DIRECTION  CATEGORIZATION 


■ 

Airplane  on  Leeward 

Side  of  Ridge 

mSSSSBM 

Wind  Direction 
When  Categorized 
Parallel  to  Ridge 

5  and  6 

157.5°  to  202.5° 

and 

337-5°  to  22.5° 

202.5°  to  247.5° 

and 

292.5“  to  337.7° 

247.5°  to  292.5° 

0°  to  180’ 

7 

137.5*  to  182.5° 
and 

317.5*  to  2.5° 

182.5°  to  227.5* 

and 

272.5°  to  317.5° 

227.5°  to  272.5° 

0°  to  l80° 

8 

137.5°  to  182.5° 
and 

337*5°  to  22.5° 

182. 5°  to  227.5" 

and 

272.5°  to  337-5° 

227.5°  to  272.5“ 

0°  to  160° 

Unfortunately  only  a  small  amount  of  data  were  recorded  when  the  tfind 
direction  was  either'  parallel  to  the  ridge  or  the  airplane  was  on  the 
windward  side  of  the  ridge.  The  prevailing  wind  direction  over  all  legs 
was  such  that  the  airplane  was  usually  on  the  leeward  sides  of  the  ridges. 
This  is  probably  one  reason  for  the  higher  gust  velocity  rms  values  re¬ 
corded  at  Peterson  than  at  the  other  three  locations.  Only  approximately 
17  per  cent  of  the  vlnd  direction  data  recorded  over  legs  5  >  6,  7>  and  8 
at  Peterson  were  classified  as  parallel  or  windward.  Although  these  data 
are  few,  there  are  enough  to  indicate  that  lower  gust  velocity  rms  values 
were  recorded  during  these  wind  conditions  on  legs  5  and  7  than  during 
conditions  when  the  wind  direction  was  categorized  as  perpendicular  or 
diagonal.  Also,  the  correlation  between  gust  velocity  r*ns  values  and  wind 
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speed  was  less  for  data  recorded*  on  the  windward  sides  of  the  mountain 
ridges.  This  agrees  with  the  results  shown  in  Figure  14.9  regarding  the 
effects  of  terrain  on  the  relationship  between  gust  velocity  nos  and  wind 
speed.  The  windward  and  parallel  data  recorded  on  leg  6  are  mostly  clus¬ 
tered  at  the  low  wind  speeds  and  are  Inconclusive. 

Leg  8  is  classified  as  a  plains  leg  but  data  were  evaluated  with  respect  * 
to  wind  direction  bemuse  of  its  proximity  to  a  predominant  mountain  ridge 
and  rough  terrain.  The  gust  velocity  rms  values  show  a  different  correla¬ 
tion  with  wind  speed,  depending  upon  how  the  wind  direction  was  categorized. 
The  gust  rms  values  for  the  perpendicular  sir’  diagonal  grouping  generally 
increased  slightly  with  increasing  wind  speed.  An  opposite  trend  is  noted 
f6r  the  parallel  and  windward  data,  and  these  data  contain  greater  gust 
velocity  rms  valuer  at  the  lower  wind  speeds. 

For  data  recorded  on  legs  5,  6,  and  7>  the  gust  velocity  rms  values,  and 
their  relationships  to  perpendicular  and  diagonal  wind  speeds,  are  signif¬ 
icantly  larger  for  the  lateral  component  than  for  either  the  vertical  or 
the  longitudinal  components  of  turbulence.  This  is  probably  became  the 
airplane  was  flown  parallel  to  the  ridges  on  these  legs  and,  therefore, 
the  predominant  winds  blowing  across  these  ridges  were  side  winds  with 
respect  to  the  airplane.  The  peaks  and  openings  along  these  ridges  would 
cause  an  Intermittent  side  wind,  or  lateral  wind  component,  and  increase 
the  lateral  gust  velocity  rms  values. 

Relationships  between  vertical  gust  velocity  rms  values  and  wind  speed 
were  investigated  in  Reference  14.1.  This  reference  suggests  that  for 
flight  over  smooth  terrain,  the  rms  gust-wind  speed  relationship  varies 
primarily  as  a  function  of  atmospheric  suability,  while  over  mountainous 
terrain,  the  effects  of  stability  do  not  appear  to  be  significant.  The 
LO-LOCAT  Phase  III  data  were  also  investigated  with  respect  to  terrain 
and  stability.  These  data  are  shown  in  Figure  14.30  where  the  slope  of 
the  least-square  lines  of  wind  Versus  gust  velocity  rms  data  are  plotted. 

The  data  were  plotted  for  all  three  components  of  turbulence.  The  only 
significant  variation  with  stability  occurs  for  data  obtained  over  the  low 
mountains.  Ter  these  data,  the  relationship  between  wind  and  rms  gusts 
becomes  less  ai*  the  atmosphere  becomes  more  unstable.  The  relationship 
over  high  mountains  is  slightly  leas  during  very  stable  conditions  than 
during  conditions  of  stable,  neutral,  and  unstable  atmospheric  stability. 
Stability  see.08  to  have  very  little  effect  on  the  relationship  of  data 
obtained  over  the  plains. 


152 


ITY  RMS  VALUE  -  f pa  GUST  VELOCITY  RMS  VALUE  -  /pa  GUST  VELOCITY  RMS  VALUE 


14 


•  .  *  . 

^  «•  •  • 


Longitudinal  Component 

- at  -  3.30  +  0.070  w 

- Standard  Deviation  *  1.92  fps 


0  10  20  30 

lUr 


60  70  80 


*  •  ix 

;y>^ . 
'  . ' 


c.  *-v..  . 

0  ^  «  .  *  • 


p  :  •  \  .  V.  ,  .  •  •  .  • 

"  V"  . 


0* - 1 - ^ - * — 

0  10  2  0  30  40 


Lateral  Component 

-at  -  3.26  +  0.093  w 

-  Standard  Deviation  »  2.24  fps 

_ L  _ I - J 

60  70  80 


.  *.\t> —  .  ‘  •  .- 

6  •;y/C-:  ••  " 

•  *T  '' 

1*  •  *  V ~<r, ■*" 

‘  "  Vertical  Component 

2 k-yj&h  ^  - at  -  5.30  +  O.O58  W 

.*.■  *'  *  - Standard  Deviation  *  1.62  fps 


_J - X  — 


0  10  20  30  1*0  50  5o  70  go 

WIHD  SPEED  -  fpa 

Figure  l4.1  Oust  Velocity  RME  Values  Versus  Wind  Speed  at 
250  Feet  Over  High  Mountains 


GUST  VELOCITY  RMS  VALUE  -  fpa  OUST  VELOCITY  RMS  VALUE  -  fp*  GUST  VELOCITY  RMS  VALUE 


*IKD  SPEED  .  fp* 


Fifur*  14.2  Oust  Velocity  IMS  Valu«*  V*r*u*  Wind  Sp**d  at 
750  ?**t  O r*r  Righ  Mountain* 

15* 


69  fps 


.00  fos 


1.53  fP« 


0  10  20  )0  ko  50  6o  TO  00 


Vito)  SPUD  .  f pm 

T lfur*  lV.3  (hut  Velocity  MB  Valoea  Veraua  Vlad  Speed  at 
250  Feet  Over  Low  Houattlna 


155 


GUS1  VTUXITY  rats  VALUE  -  t pa  01' ST  VELOCITY  RMS  VALUE  -  fpe  GUST  VELOCITY 


14 

12|- 


10f 


A 


Lateral  Component 
_ (jt  -  2.00  +  0.058  W 

-  Standard  Deviation  »  1.01  fpe 


6 

U- 


t- 


v 

oL 

0 

lUr 


r* 

<~t  t  ' 

..  -  t  JL  1  I 

10  20  50  i*0 


50  60  TO  80 


12 

10 

8 

6 


Vertical  Component 

- ot  -  1.5*  ♦  0.066  w 

- Standard  Deviation  *  0.98  fpe 


i  - 

kb  • 


i  V  te  * 

\  *  «  *  •  , 

^  •  •  _  «  ■  - 

;  j.  ' 

of.  •  -  ;  t\A  -  h 
■ 


St 


(>l _ 1 - -  l  - ...  I  t  •  -  i  1  *■  •  -  1 

0  10  20  50  ko  50  60  70  80 

VIJTO  SPEED  .  fpa 

figure  lfc.V  Guot  Velocity  MB  Valeee  Versua  tied  S 


peed  at 


750  Feet  Ovei  Lee  Mount alas 


GUST  VELOCITY  RMS  VALUE  -  fps  GUST  VELOCITY  RMS  VALUE  -  fos  GUST  VELOCITY  RMS  VALUE 


GUST  VELOCITY  RMS  VALUE  -  fps  GUST  VELOCITY  RMS  VALUE  -  fps  GUST  VELOCITY  RMS  VALUE 


I 


J.U| 
12 1 
10 

8^ 


Longitudinal  Component 
— ot  -  2.77  +  0.017  w 
- Standard  Deviation  ■  1,19  fps 


12 

10 

8 


Lateral  Component 
— —  ot  ■  2.53  +  0.02U  W 
- Standard  Deviation  =  1.13  fps 


12  h 

„! 

6: 


ertical  Component 
-  Gt  3  2.40  +  0.010  W 

-  Standard  Deviation  a  1.03  fps 


0- 

0  10  20  30  ko  50  60  70  80 

WIND  SPEED  -  fps 


Figure  l1*, 7  Oust  Velocity  RHS  Values  Versus  Wind  Speed  at 
250  Feet  Over  Plains 


159 


GUST  VELOCITY  RMS  VALUE  -  f ps  GUST  VELOCITY  RMS  VALUE  —  fps  GUST  VELOCITY  RMS  VAILE 


& 

Vi 


lUr 


i  12 


10- 

8- 


Longitudinal  Component 

_ at  ■  2.83  +  0.004  w 

- Standard  Deviation  -  0.93 


6 


°0  10  20  30  1*0  50  60  70  80 

l4r 


12 

10 

e 


Lateral  Component 

_ at  *»  2.8o  +  o.oio  w 

- Standard  Deviation  =  0.Q6  fps 


I 

12  r 

i 


4o  50  60  70  80 


Vertical  Component 

0 —  at  a  3.16  -  0.006  w 

- Standard  Deviation  «  1 ,03  fps 


,u  -  .  .  *  i  ‘  •  »  ' 

0  10  2  0  30  1*0  50  6o  70  8o 

WIND  SPEED  -  fps 

Figure  14.8  Gust  Velocity  RMS  Vnlues  Versus  Wind  Speed  at 
750  Feet  Over  Plains 


160 


SIjOPE  OF  GUST  VELOCITY  RMS  VALUE  VERSUS  WIND  SPEED 
LEAST  SQUARE  LINE  (^c^/AW)  -  Dimensionless 


c 


Longitudinal  Component 


0.10  - 

0.05  “ 


0 


-0.051- 


P  -  Plains 
D  -  Desert 
IM  -  Low  Mountains 
HM  -  High  Mountains 


□ 

o 


_g 


o  250  Fl  et 
0  750  Feet 


0.10 

— 

Lateral  Component 

a 

0.05 

0 

— —  _ 

1 

1 

1 

OG3| 

1 

1 

0 

0 

-0.05 

— 

1 

1 

1 

i . . 

0.10 

- 

Vertical  Component 

0.05 

_ 

□ 

0 

§ 

0 

—  — 

-  -Q- 

1 

i 

I 

1 

op 

1 

1 

—  — 

—  ~  ~ 

-0.05 

- 

1 

.  1 

1 

i 

D 

p 

LM 

HM 

TYPE  OF  TERRAIN 

Figure  lU.9  Correlation  of  Gust  Velocity  RMS  Value 
Wind  Speed  as  a  Function  of  Terrain  an 
Altitude 


l6l 


Longitudinal  Component 


■rm 

.999 

■9* 

.9T> 

■  99 


Lateral  Component 


Vertical  Component 


0  <  W  <  10 
JO  <  W  <  20 


a  g~g 

i  i  > 


.I.I.l.l.I.M.l.l.I.I 
*  «  r  I  9  X  11  U  U  1*  15 

■MT  MffTTT  M  VIM  -  ((i 
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Figure  lL.l6  Correlation  Between  Vercical  Gust  Velocity  RMS 
Values  and  Wind  Speed  with  Respect  to  Wind 
Direction  -  Leg  5  at  Peterson 
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Figdre  14.17  Correlation  Between  Longitudinal  Gust  Velocity 
RMS  Values  and  Wind  Speed  with  Respect  to  Wind 
Direction  -  Leg  6  at  Peterson 
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Figure  14. 18  Correlation  Between  Lateral  Gust  Velocity  RMS 
Values  and  Wind  Speed  with  Respect  tc  Win-* 
IHrection  -  Leg  6  at  Peterson 
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Figure  14.19  Correlation  Between  Vertical  Gust  Velocity  RMS 
Values  and  Wind  Speed  with  Respect  to  Wind 
Direction  -  Leg  6  at  Peterson 
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Figure  14.20  Correlation  Between  Longitudinal  Gust  Velocity 
RMS  Values  and  Wind  Speed  with  Respect  to  Wind 
Direction  -  Leg  7  at  Peterson 
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Figure  14.21  Correlation  Between  Lateral  Gust  Velocity  RMS 
Values  and  Wind  Speed  with  Respect  to  Wind 
Direction  -  Leg  7  at  Peterson 
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Figure  lh.22  Correlation  Between  Vertical  Gust  Velocity  RMS 
Values  and  Wind  Speed  with  Respect  to  Wind 
Direction  -  Leg  7  at  Peterson 
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Figure  14.23  Correlation  Between  Longitudinal  Gust  Velocity 
EMS  Values  and  Wind  Speed  with  Respect  to  Wind 
Direction  -  Leg  8  at  Peterson 
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Figure  14.24  Correlation  Between  Lateral  Gust  Velocity  RMS 
Values  and  Wind  Speed  with  Respect  to  Wind 
„  Direction  -  Leg  8  at  *  rnon 
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Figure  14.25  Correlation  Between  Vertical  Gust  Velocity  RMS 
Values  and  Wind  Speed  with  Respect  to  Wind 
Direction  -  Leg  8  at  Peterson 
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figure  14.26  Effects  of  Wind  Direction  on  the  Correlation 
Between  Oust  Velocity  HC3  Values  and  Wind 
Speed  -  Leg  5  at  Peterson 
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Figure  lU.27  Effects  of  Wind  Direction  on  the  Correlation 
Between  Gust  Velocity  RH5  Values  and  Wind 
Speed  -  lag  6  at  Peterson 
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Figure  lU.28  Effects  of  Wind  Direction  on  the  Correlation 
Between  Gust  Velocity  RMS  Values  and  Wind 
Speed  -  Leg  7  at  Peterson 
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Figure  lU.29  Effects  of  Wind  Direction  on  the  Corre1ation 
Between  Gust  Velocity  RMS  Values  and  Wxnd 
Speed  -  Leg  8  at  Peterson 
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Figure  lb. 30  Correlation  of  Oust  Velocity  BMS  Values  with 
Wind  Speed  as  a  Function  of  Terrain  and 
Atmospheric  Stability 
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15.  MAXIMUM  DERIVED  EQUIVALENT  GUST  VELOCITIES 


An  approximate  maximum  and/or  minimum  value  of  derived  equivalent  gust 
velocity  was  calculated  for  each  turbulence  sample.  These  values  were 
used,  to  establish  the  approximate  turbulence  level  encountered  and  as 
an  immediate  guideline  in  determining  if  gust  encounters  were  of  a  mag¬ 
nitude  to  warrant  an  airplane  structural  inspection.  Maximum  incremen¬ 
tal  acceleration  from  a  one  g  mean  was  determined  from  "quick-look" 
strip  charts  of  vertical  acceleration  at  the  airplane  center  of  gravity. 
Confutation  of  the  derived  equivalent  gust  velocity  was  accomplished 
using  the  equation: 


2  W _ 

Ude  _CL<  p0  SV,  Kg  , 


ANzcg 


where  Kg,  the  gust  alleviation  factor,  is  expressed  as: 


.86* 

5.3+* 


(15.1) 


(15.2) 


and  fi ,  the  dimensionless  airplane  mass  ratio,  is  defined  by: 


2W 

M  *  CjJ>cgs 


(15.3) 


Using  the  equation  of  state,  p  (air  density)  may  be  expressed  as: 


(15.4) 


The  values  of  Ude  equation  constants  for  the  T-33A  airplane  are  defined 
below: 

c  =  mean  aerodynamic  chord  =  6.72  feet 
g  =  gravitational  constant  =  32.2  feet/sec2 
pQ  =s  standard  day  air  density  =  .002378  slug/feet3 
S  =  wing  area  =  238  feet2 

Constant  values  were  determined  for  the  equation  variables  to  facilitate 
calculations  in  the  field.  This  was  accomplished  by  using  an  average  gross 
weight  for  the  airplane  over  the  route  ( 12, 500  pounds),  an  equivalent  air¬ 
speed  approximated  by  target  airspeed  (591.15  feet  per  second),  and  the 
slope  of  the  lift  coefficient  versus  angle  of  attack  curve  at  a  Mach  num¬ 
ber  equivalent  to  target  airspeed  and  average  altitude  (5>9l/r^di.an) .  The 
ratio  of  static  pressure  to  ambient  temperature  was  estimated  from  data 
obtained  during  Phases  I  and  II.  This  ratio  was  calculated  by  averaging 
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these  data  over  a  corresponding  seasonal  time  period  for  each  similar 
route  leg.  The  values  for  this  ratio  and  the  simplified  Ude  equation 
derived  by  substitution  of  the  assigned  constants  into  Equation  15.1 
are  given  In  Table  15.1. 

Cumulative  probability  distributions  of  Ude  values  calculated  during 
Phase  III  are  shown  in  Figures  15.1  through  15.4.  A  comparison  of  the 
Ude  cumulative  probability  distributions  to  those  of  the  maximum  recorded 
gust  velocities  is  discussed  in  the  next  section. 


The  absolute  values  of  the  maximum  derived  equivalent  gust  velocities  are 
tabulated  in  Appendix  VI  for  each  turbulence  sample. 

TABLE  15.1 


Ude  EQUATIONS  AND  EQUATION. CONSTANTS 
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Gust  Velocity  and  Ude  -  All  Categori  a 


(issijg  as  '*&$$$$  I  §  I  i  fe  §  e  B  SB  *  «  «  *  ®  ®  as  *  *  $  *$|  | 


16.  MAXIMUM  GUST  VELOCITY  DURING  EACH  SAMPLE 


The  maximum  absolute  values  of  each  gust  velocity  component  for  each  tur¬ 
bulence  sample  were  obtained.  Absolute  values  not  calculated  because  of 
low  signal-to -noise  ratios  were  represented  by  average  values  estimated 
from  typical  low  sigma  turbulence  samples.  The  cumulative  probability 
distributions  of  these  data,  as  well  as  of  the  maximum  absolute  derived 
equivalent  gust  velocities  were  calculated  and  are  shown  in  Figures  15.1 
through  15.4.  Figure  15.1  shows  the  overall  cumulative  probability  as 
obtained  from  all  categories;  Figure  15.2,  the  probability  with  respect 
to  terrain;  Figure  15. 3*  the  probability  with  respect  to  altitude;  and 
Figure  15.4,  the  probability  with  respect  to  time  of  day.  Comparison  of 
the  vertical  gust  velocity  component  with  the  maximum  derived  equivalent 
gust  velocity  indicates  fairly  good  agreement  in  practically  all  cases. 
Cong&rlson  of  the  three  calculated  gust  velocity  components  showB  that 
the  cumulative  probability  distributions  of  the  vertical  and  longitudinal 
components  are  very  similar  for  all  the  categories  presented.  The  lateral 
component,  however,  shows  probabilities  higher  than  the  other  components 
for  the  overall  category,  the  high  mountain  category,  the  time  of  day 
category,  and  the  altitude  category.  The  high  mountain  data  is  included 
in  these  categories  and  it  appears  that  the  cumulative  probability  dis¬ 
tribution  of  the  lateral  component  is  higher  than  that  of  the  longitudi¬ 
nal  or  vertical  because  of  flying  parallel  to  the  high  mountain  ridge  lines. 

Other  factors  to  be  considered  when  congaring  gust  velocity  conponents  are: 

•  High-pass  filtering  of  the  gust  velocity  components  attenuates 
the  magnitude  of  the  long  wavelength  amplitudes  of  the  longi¬ 
tudinal  component  slightly  more  than  the  vertical  or  lateral 
components  because  the  amplitudes  of  the  longitudinal  compo¬ 
nent  are  higher  at  the  long  wavelengths,  as  discussed  in  Sec¬ 
tion  V  . 

•  The  long  wavelength  amplitudes  of  the  vertical  component  are 
smaller  than  those  of  the  lateral  or  longitudinal  components 
because  of  the  airplane's  close  proximity  to  the  ground, 
allowing  measurement  of  the  compressed  wavelength  vertical 
eddys. 

•  Incomplete  removal  of  the  effects  of  airplane  motion  from  the 
gust  velocity  components  tends  to  increase  the  level  of  the 
lateral  component  during  flight  over  plains,  water,  aud  desert 
and  would  also  tend  to  increase  the  level  of  both  the  vertical 
and  lateral  components  during  contour  flight  over  the  high  and 
low  mountain  legs.  This  effect  is  believed  to  be  the  major 
source  of  error  in  true  gust  velocity  calculations.  The  mag¬ 
nitudes  of  these  error s  are  small  and  are  discussed  in  Appen¬ 
dix  III. 
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Time  histories  of  the  longitudinal ,  lateral,  and  vertical  gust  velocity 
consonants  were  plotted  for  turbulence  samples  during  which  the  maximum 
value  of  one  or  more  of  the  components  equaled  or  exceeded  ±50  feet  per 
second. 

Sandies  meeting  the  above  criteria  were  not  encountered  over  the  McConnell 
or  Griff is s  routes.  One  turbulence  sample  meeting  these  criteria  was  en¬ 
countered  over  the  Edwards  route  while  seventy  separate  occurrences  were 
encountered  during  flight  over  the  Peterson  route.  The  criteria  were  ex¬ 
ceeded  only  on  those  samples  flown  over  terrain  classified  as  high  moun¬ 
tain.  Time  histories  of  the  gust  velocity  consonants  for  the  single 
Edwards  encounter  are  presented  in  Figure  17.1  and  those  for  the  Peterson 
occurrences  are  shown  in  Figures  17.2  through  17.71*  A  data  block  Included 
on  each  figure  identifies  the  test  number  and  the  leg  over  which  these  data 
were  obtained.  Also  tabulated  ir  the  data  block  are  the  maximum  and  mini¬ 
mum  gust  velocities,  turbulence  time  series  standard  deviation,  and.  If 
available,  the  average  wind  direction  and  velocity,  average  ambient  air 
temperature,  average  true  airspeed,  average  radar  altitude,  and  category 
number. 

The  vertical  gust  velocity  conqxsient  in  Figure  17.1  should  be  particularly 
noted.  These  data  show  a  turbulence  wavelength  of  approximately  14,300 
feet  with  peak  magnitudes  of  74.5  and  -65*0  feet  per  second.  The  long 
wavelength  consonant  appears  at  a  frequency  of  approximately  0.04  cycles 
per  second  in  the  time  history.  This  same  turbulence,  if  encountered  dur¬ 
ing  Phases  I  and  H  testing,  would  have  appeared  as  0.02  cycles  per  second 
data  due  to  the  lower  speed  at  which  the  C-131  aircraft  were  flown.  Phases 
I  and  II  gust  velocity  measurements  in  this  frequency  regime  would  have  been 
attenuated  approximately  80  per  cent  due  to  the  low  frequency  data  filtering 
techniques  utilized  to  alleviate  the  effects  of  instrumentation  drift.  This 
is  a  cogent  exaiqple  of  the  fact  that  there  is  turbulence  at  low  altitude  con¬ 
taining  wavelengths  longer  than  7/000  feet  and  is  the  reason  why,  in  compar¬ 
ing  Phase  HI  results  with  Phases  I  and  II  for  equivalent  categories,  that: 

(l)  the  probability  distributions  of  gust  velocity  peaks  are  higher,  (2)  the 
probability  distributions  of  gust  velocity  rms  values  are  higher,  (3)  integral 
scale  lengths  are  longer,  (4)  Taylor’s  microscales  are  longer,  and  (5)  the 
viscous  dissipation  rates  are  less. 

The  number  of  Peterson  turbulence  samples  and  the  number  of  samples  during 
which  a  particular  gust  velocity  component  met  the  50  feet  per  second  cri¬ 
teria  are  tabulated,  with  respect  to  leg  number,  in  Table  17.1.  This  tabu¬ 
lation  shown  that  the  totaj.  number  of  severe  turbulence  encounters  over  leg 
five  exceeded  that  of  the  other  legs  by  a  factor  greater  than  two.  As 
shown  in  Appendix  IV,  leg  five  cro  cedes  south  along  the  eastern  slope  of 
the  Sangre  de  Cristo  Oulebra  Range,  The  wind  direction  at  the  time  of 
these  encounters  varied  from  254  to  319  degrees;  therefore,  the  airplane 
was  an  the  lee  side  of  the  range  where  the  most  severe  turbulence  would 
most  likely  occur.  Also,  on  leg  five,  the  lateral  component  exceeded  50 
feet  per  second  during  more  turbulence  samples  than  did  the  longitudinal 
or  vertical  components. 


TABUS  17.1 


TABULATION  OF  SAMPLES  AND  0OMPONENTS 
EQUAL  TO  OR  EXCEEDING  $0  FPS 


The  above  set  of  conditions  for  leg  5  were  analogous  to  those  encountered 
on  legs  1,  4,  6,  and  7,  viz;  for  turbulence  samples  gathered  on  the  lee 
side  of  a  mountain  range,  the  lateral  component  exceeded  50  fps  during 
more  samples  than  did  the  vertical  or  longitudinal  component.  It  should 
be  noted  that  aiiaost  always  the  lateral  component  exhibited  higher  posi¬ 
tive  than  negative  excursions  when  the  airplane  was  flying  south  (legs  4 
and  5)  and  higher  negative  than  positive  excursions  when  the  airplane  was 
flying  north  (legs  1,  6,  and  j) .  The  lateral  component  is  defined  as  pos¬ 
itive  for  gusts  from  the  left  of  the  airplane. 

The  longitudinal  component  for  leg  2,  Figures  17.10  through  17*  l8,  exceeded 
50  feet  per  second  during  more  severe  turbulence  samples  than  did  either 
the  lateral  or  vertical  conponants .  This  again  is  probably  due  to  the  south 
westerly  direction  of  the  leg  with  respect  to  the  wind  direction.  The  wind 
direction  for  the  time  histories'  shown  varied  from  24l  to  306  degrees  with 
one  sample  showing  winds  from  47  degrees. 

The  striking  similarity  of  the  gust  velocity  components  depicted  in  Figures 
i7.ll,  17.12,  and  17.13  should  also  be  noted.  These  turbulence  samples 
were  gathered  over  leg  2  of  the  Peterson  route  on  threb  separate  flights 
conducted  on  the  same  day.  The  elapsed  time  between  the  first  turbulence 
sample  gathered  during  the  dawn  flight,  and  the  last  turbulence  sample  gath¬ 
ered  during  the  mid-afternoon  flight  is  approximately  seven  hours.  Time 
histories  having  the  severe  turbulence  characteristics  similar  to  the  three 
mentioned  above  were  also  obtained  19  and  23  days  later  as  shown  in  Figures 
17.14  and  17.18. 

The  most  severe  continuous  turbulence  encountered  occurred  over  leg  5  and 
is  shown  in  Figure  17-38.  The  time  series  standard  deviations  were  15-45 
fps,  19.05  fps,  and  3.5-04  fps  for  the  longitudinal,  lateral,  and  vertical 
gust  velocity  coaponents,  respectively.  None  of  the  maximum  excursions 
occurred  during  this  particular  turbulence  sample.  These  maxlmus  values 
occurred  as  shown  in  Figures  17-1  for  the  vertical  component,  in  Figure 
17.33  longitudinal  component,  and  in  Figure  17-34  for  the  lateral  component. 
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Figure  17*1  Time  History  of  Severe  Turbulence  Encounter 


Figure  17.2  Tine  History  of  Severe  Turbulence  Encounter 
Peterson 
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Figuro  17.3  Tl»e  History  of  Severe  Turbulence  Encounter 
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Figure  17.5  Tine  History  of  Severe  Turbulence  Encounter 
Peterson 


Figure  17.8  Tine  History  of  Severs  Turbulence  fiicounter 
Peterson 


s' 

i 


l 


i 


200 


Figure  17.10  Tine  History  of  Severe  Turbulence  Encounter 
Peterson 
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Figure  17.11  Tine  History  of  Sever*  Turbulence  Encounter 
Peterson 
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Figure  17.12  Tine  History  of  Severe  Turbulence  Encounter 
Peterson 
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Figure  17.13  Tlae  History  of  Severe  Turbulence  Encounter 
Peterson 


Figure  17. lU  Time  History  of  Severe  Turbulence  Encounter 

Peterson 
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Figure  17.17  ^*l»o  History  of  Severe  Turbulence  Encounter 

Peterson 


Figure  17.18  Tine  History  of  Severe  Turbulence  Encounter 
Peterson 


Figure  17.19  Tin  History  of  Severe  Turbulence  Bn counter 
Peterson 
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Tigurm  17.21  Tine  History  of  Severe  Turbulence  Encounter 
Peterson 


Twt  196  Ambient  Air  Temp  (*F)  37. ^  u  iu  (*P*)  30.0 
D*t «  2-26-69  True  AlrapMd  (fps)  6l?.2^  u  aln  (*P*)  -52.2 
Leg  Mo.  4  Radar  Altitude  (ft)  No  Data  v  aax  (*P*)  6U.6 
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Flyura  17.22  Time  History  of  Severe  Turbulence  Encounter 
Peterson 


Figure  17*2U  TIm  History  of  Severe  Turbulence  Encounter 
Peterson 


Figure  17.25  Tine  History  of  Severe  Turbulence  Encounter 
Peterson 
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FJg-ure  17.27  Ti«e  History  of  Severe  Turbulence  Encounter 


Figure  17.28  Tlee  History  of  Severe  Turbulence  Encounter 
Peterson 
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Figure  17.29  Time  History  of  Sever*  Turbulence  Encounter 
Peterson 
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Fifurt  17 *3^  ^  Ti®®  History  of  Stvoro  Turbultnco  Encounter 


Figure  17.52  TIM  History  of  Severe  Turbulence  Encounter 
Peterson 
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Plfuro  17.36  TIm  Hlttory  of  Severe  Turbulence  Encounter 
Pet ere on 
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7.38  Time  History  of  Severe  Turbulence  Encounter 


230 


Figure  17.39  Tl»e  History  of  Severe  Turbulence  Encounter 
Peterson 


Figure  17.40  Ti»e  History  of  Severe  Turbulence  Encounter 
Peterson 
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Eilfure  17.^3  Tttn  History  of  Severe  Turbulence  Encounter 
Paterson 
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Figure  17.44  TIm  History  of  Severe  Turbulence  Encounter 
Peterson 
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Figure  17.46  Tl®e  History  of  Severe  Turbulence  Kncomter 
Peterson 
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Figure  17.50  Time  History  of  Severe  Turbulence  Encounter 
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Figure  17.51  TiiM,  History  of  Severe  Turbulence  Encounter 
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Peterson 


Figure  17 « 5^  Ti»e  History  of  Severe  Turbulence  Encounter 

Peterson 
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Figure  17.55  Tiee  History  of  Severe  Turbulence  Encounter 
Peterson 
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Fifur •  17.56  Tl»  History  of  Severe  Turbulence  Encounter 
Peterson 
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Figure  17.56  Tine  History  of  Severe  Turbulence  Encounter 
Peterson 
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Tiwe  History  of  Severe  Turbulence  Encounter 
Peterson 
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Figure  17.60  Ti*e  History  of  Severe  Turbulence  Encounter 
Peterson 
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Flgfure  17.6l  Tine  History  of  Severe  Turbulence  Encounter 
Peterson 


Figure  17. 6e  Tine  History  of  Severe  Turbulence  Encounter 
Peterson 
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Figure  17.63  Time  History  of  Severe  Turbulence  Encounter 
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Figure  17. 6U  Tine  History  of  Severe  Turbulence  Encounter 
Peterson 
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Figure  17.65  Time  History  of  Severe  Turbulence  Encounter 
Peterson 
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Figure  17.67  Tine  History  of  Severe  Turbulence  Encounter 
Peterson 


Figure  17.68  Ti»e  History  of  Severe  Turbulence  Encounter 
Peterson 
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Figure  17.69  Ti»e  History  of  Severe  Turbulence  Encounter 
Peterson 


Fljur©  17.70  TIM  History  of  Sever©  Turbulence  Encounter 
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Figure  17.71  Ti»e  History  of  Severe  Turbulence  Encounter 
Peterson 


16.  OUST  VELOCITY  MAGNITUDE  VERSUS  TEMPERATURE  CORRELATIONS 


Since  the  advent  of  high  altitude  commercial  and  military  flight  activity, 
at  least  one  program  (Reference  l8,l)  has  been  accomplished  to  relate  tur¬ 
bulence  to  air  temperatures.  The  objective  in  determining  this  relation¬ 
ship  was  to  predict  the  occurrence  of  CAT  so  that  evasive  action  could  be 
implemented.  The  feasibility  of  determining  a  correlation  of  turbulence 
with  both  air  and  ground  surface  temperatures  was  evaluated  during  the 
LO-LOCAT  Program.  Throughout  the  program,  air  and  surface  temperatures 
were  recorded. 

The  coherency  function  between  indicated  OAT  and  vertical  gust  velocity 
and  GST  and  vertical  gust  velocity  was  calculated  as  shown  by  Equation 
18.1.  This  function  indicates  the  degree  of  statistical  correlation  be¬ 
tween  two  variables.  A  value  of  1.0  would  indicate  perfect  correlation 
and  a  value  of  zero  would  Indicate  no  correlation. 

„  (18.1) 

y  (fc)  a  I  —  s 

xy  +r(k)  +y(k) 

Four  turbulence,  samples;  one  in  each  stability  class  recorded  at  250  feet 
were  chosen  for  this  analysis.  The  coherency  function  for  OAT  and  GST 
with  respect  to  the  vertical  gust  velocity  component  is  shown  in  Figures 
l8.1  through  16. h.  The  degree  of  correlation  1b  very  low  for  GST  and  only 
slightly  higher  for  OAT  at  frequencies  between  10 and  10'3  cpf.  ^ 

It  is  concluded  that  OAT  and  GST  values  recorded  during  this  program  are 
not  correlated  with  gust  velocity  rms  values.  This  nay  nave  been  due  to 
certain  limitations  in  this  data.  Poor  frequency  response  characteristics 
of  the  transducers  would  adversely  effect  any  correlation  of  temperature 
with  gust  velocity. 

The  frequency  input  to  the  radiometer  used  to  measure  ground  surface  tem¬ 
perature  could  be  altered  by  the  speed  of  the  airplane  when  flying  above 
a  terrain  whose  features  vary  significantly.  If  the  radiometer  has  rela¬ 
tively  slow  response  characteristics,  the  speed  can  have  considerable 
effect  on  surface  temperature  measurements.  Figure  18. 5  shows  the  change 
of  the  PSD  slope  of  GST  with  respect  to  ground  speed  for  several  data 
samples  recorded  over  various  "plains  like"  terrain  segments.  The  slopes 
shown  in  this  figure  were  obtained  over  that  portion  of  the  PSD  frequency 
range  where  the  power  spectra  decreased  in  an  approximate  logarithmic 
faBhion.  This  terrain  type  was  chosen  because  of  the  sharp  contrasts  in 
surface  features  which  gave  tenperature  changes.  This  figure  shows  that 
as  the  ground  speed  increases,  the  -lope  of  the  ground  surface  tempera¬ 
ture  PSD  also  increases.  This  indicates  a  frequency  response  problem 
since  increasing  speed  and  hence  increasing  frequency  increases  the  amount 
of  attenuation  of  the  higher  frequencies. 

The  PSD  of  indicated  ambient  air  temperature  out  to  0.0065  cpf  in  Figure 
18.6  shows  that  the  response  is  near  a  -5/3  logarithmic  slope,  which  is 
what  would  be  expected  /Reference  lQ.2).  The  low  frequency  true  airspeed 
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variations  are  evident  in  the  higher  power  deviations  above  the  -5/3  slope 
since  this  parameter  has  the  temperature  effects  due  to  ram  rise  included. 
Because  of  these  variations  due  to  airspeed,  it  is  probable  that  true 
ambient  temperature  could  provide  a  more  meaningful  correlation  with  gust 
velocity. 

Standard  deviations  of  ground  surface  temperature  (GST)  and  vertical  gust 
velocity  for  the  plains  leg  at  Peterson  Field  are  compared  in  Figure  18.7. 
This  figure  indicates  that  there  is  no  direct  relationship  between  the 
standard  deviations.  It  dees  show  the  increase  in  gust  velocity  rms  with 
decreasing  stability  noted  throughout  the  analysis  of  LO-LOCAT  data.  A 
similar  comparison  of  OAT  standard  deviation  versus  gust  velocity  rms 
values  showed  a  greater  degree  of  scattep. 

It  is  apparent  that  the  radiometer  and  the  OAT  data  should  be  corrected 
for  frequency  response  before  any  valid  correlation  to  gust  velocities 
can  be  determined. 
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Coherency  of  Vertical  Oust  Velocity  vith  Outside  Air 
and  Ground  Surface  Temperatures 
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Coherency  of  Vertical  Gust  Velocity  with  Outside  Air 
and  Ground  Surface  Temperatures 


Ground  Surface  Temperature  Standard  Deviation  Versus 
Gust  Velocity  Standard  Deviation 


SECTION  V 


GUST  VELOCITY  SPECTRA 


Pove  s  ectra  were  calculated  for  904  of  the  1938  total  four  and  one-half 
minute  turbulence  samples  recorded  during  this  program. 

Of  the  904  samples,  185  were  discarded  because  the  data  from  all  three 
components  exhibited  low  signal-to-noise  ratios,  instrumentation  anomalies, 
or  high  coherency  and  an  additional  98  were  discarded  because  all  three 
conponents  showed  non -homogeneity  characteristics.  This  resulted  in  a 
total  of  621  samples  remaining  for  analysis  of  the  above  characteristics 
for  each  gust  velocity  component;  4-55  of  these  samples  were  considered  to 
exhibit  valid  data  for  all  components  and  136,  69,  and  4l  additional  sam¬ 
ples  were  considered  valid  for  the  individual  longitudinal,  lateral,  and 
vertical  components,  respectively.  Thus,  broken  down  into  the  individual 
conponents,  the  total  nunfcer  of  samples  considered  valid  with  respect  to 
the  longitudinal,  lateral,  and  vertical  components  were  591>  524,  and  496, 
respectively.  These  latter  numbers  of  samples  were  used  to  obtain  normal¬ 
ized  spectra,  truncated  standard  deviations,  comparison  of  experimental- 
mathematical  Spectra,  and  average  turbulence  scale  lengths  while  the  455 
number  of  samples  were  used  in  evaluating  coherency  and  isotropy.  Power 
spectral  density  data  for  these  455  conditions  are  presented  in  Appendix  3X. 

19.  INSTRUMENTATION  ANOMALIES,  LOW  INTENSITY  TURBULENCE 

An  example  of  power  spectra,  with  suspected  instrumentation  problems,  Is 
presented  in  Figure  19.I.  It  is  obvious  that  the  spectra  are  suspect  due 
to  the  large  hump  in  the  lateral  spectrum,  which  occurs  at  the  airplane 
short  period  mode.  Experience  has  shown  that  when  spectra  such  as  that 
shown  in  Figure  19.I  are  obtained,  it  can  ultimately  be  proven  to  be 
caused  by  instrumentation  problems. 

An  example  power  spectra  plot  of  low  intensity  turbulence  is  presented  in 
Figure  19.2.  Turbulence  data  sanples  containing  low  signal-to-noise  ratios 
typically  had  values  less  than  1.4  fps,  irregularly  shaped  spectra  (i.e., 
spectra  which  do  not  have  a  definite  -5/3  logarithmic  slope  and  in  which 
the  power  turns  upward  at  the  higher  frequencies),  and  uncocmonly  large 
integral  scale  lengths. 
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20.  STATISTICAL  CONL  IEENCE  AND  IKCEPENIENCY 

■  *■  —  -  m  mm  _ _ 

Confidence 

Turbulence  sample  lengths  were  established  to  give  a  compromise  between 
statistical  confidence  and  turbulence  stationarlty  considerations.  Alias¬ 
ing  was  prevented  by  the  use  of  a  high  sampling  rate  to  give  a  large  fre¬ 
quency  span  between  the  break-point  of  the  low-pass  filter  and  the  Nyquist 
frequency.  Gust  velocities  were  calculated  using  a  sampling  rate  of  100 
sps  which  results  in  a  Nyquist  frequency  of  50  cps.  PSD's  were  calculated 
by  averaging  adjacent  points  which  resulted  in  a  sampling  rate  of  50  sps 
for  the  PSD  calculations.  The  break-point  of  low-pass  filters  used  was 
22  cps  the  Nyquist  frequency  for  the  PSD  calculations  was  25  cps. 

With  negligible  aliasing,  the  stability  of  uniformly  spaced  power  spectral 
density  estimates  may  be  Judged  by  analogy  with  a  chi-square  variate  with 
d  degrees  of  freedom  where: 


d*  (20. l) 


With  the  number  of  data  points  n  =  13,500,  A  f  =  0.04l£  cps  and  f„  =  25  cps, 
approximately  45  degrees  of  freedom  were  obtained  for  individual  samples. 

Statistical  reliability  was  also  evaluated  by  defining  an  interval  of  con¬ 
fidence  on  either  side  of  the  measured  spectrum.  The  interval  indicates 
the  amount  of  variation  that  may  be  expected  in  the  spectrum. 

Based  on  90  per  cent  confidence  the  equivalent  degrees  of  freedom  are  de¬ 
fined  In  Reference  20.1  as: 


(20.2) 

(20.3) 

(20.4) 
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and: 


If  it  is  assumed  that  the  pooling  of  turbulence  samples  in  a  category  is 
equivalent  to  placing  them  end  to  end,  then  the  statistical  confidence  for 
averaged  spectra  ie  quite  high.  Degrees  of  freedom  are  increased  by  a 
factor  equal  to  the  number  of  sauries  pooled.  Confidence  bands  are  reduced 
as  shown  in  Figure  20.1.  This  occurs  since  m  remains  the  same  and  n  is  in¬ 
creased  by  a  factor  equal  to  the  number  of  samples  pooled.  r 

Independency 

The  statistical  independence  of  the  gust  components  u,  v,  and  w  was  evalu¬ 
ated  \asing  the  coherency  function.  The  coherency  function  between  u  and  v, 
end  v  and  w,  was  calculated,  e.g.: 


Ok) 


[vwf 

*,(*) 


(20.13) 


Two  parameters  are  said  to  be  completely  dependent  (coherent)  if  their  co¬ 
herency  function  is  equal  to  1.0  and  completely  independent  (not  coherent) 
if  the  coherency  is  equal  to  zero.  Due  to  statistical  variations  in  the 
recorded  data,  electronic  noise,  etc.,  a  coherency  of  zero  is  impossible  to 
achieve.  The  gust  velocity  components  are  therefore  considered  to  be  inde¬ 
pendent  if  the  coherency  function  approaches  zero. 

The  coherency  function  was  evaluated  by  category.  This  was  accomplished 
by  comparing  average  data  from  different  categories.  The  average  coher¬ 
ency  function  was  calculated  by  determining  the  value  of  the  coherency 
function  of  each  sample  at  spatial  frequencies  of  5  x  10  "5 ,  7  x  10 "5, 

1  x  10  *4,  1.5  x  10'4,  2  x  10-4,  3  x  10'4,  4  x  10*4,  6  x  10 ’4,  1  x  10'3 , 

1.5  x  10‘3 ,  2  x  10'3 ,  4  x  10'3,  7  x  10*3  ,  lx  10-2 ,  and  1.5  x  10‘2  cpf. 

The  values  obtained  at  each  of  these  frequencies  were  then  averaged  and 
the  standard  deviation  calculated. 

No  significant  or  consistent  difference  in  coherency  was  associated  with 
different  categories.  The  magnitude  of  the  coherency  function  was  rea¬ 
sonably  small  indicating  statistical  independency,  especially  at  the 
shorter  wavelengths.  The  maximum  value  of  the  average  was  less  than  0.15 
as  can  be  seen  in  Figure  20.2.  These  data  represent  the  average  of  coher¬ 
ency  values  for  all  categories.  Statistical  independency  for  individual 
turbulence  samples  is  shown  in  Appendix  DC. 
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Figure  20.1  Confidence  Bands  Associated  vith  Pooled 
Sandies 


21.  HOMOGENEITY 


Only  turbulence  samples  considered  to  be  homogeneous  can  be  used  to  accu¬ 
rately  define  valid  scale  lengths  and  spectra  shape.  Those  samples  con¬ 
sidered  to  be  non-hbmogeneous  were  eliminated  from  the  data  being  used  to 
define  scale  lengths  and  spectra  shape. 

Homogeneity  Is  the  property  which  in  space  coordinates  corresponds  to  the 
statlonarlty  property  In  time  coordinates.  The  entire  analysis  of  gust 
velocity  time  histories  as  recorded  by  an  airplane  Is  based  upon  the  va¬ 
lidity  of  Taylor’s  hypothesis.  Assigning  Taylor's  hypothesis  to  be  valid, 
then  spatial  homogeneity  of  the  turbulence  sample  Is  Indicated  by  the  sta- 
tionarity  of  the  time  history.  Detailed  discussions  of  Taylor's  hypothesis 
may  be  found  In  References  21.1  end  18 .2.  An  indication  of  turbulence  homo¬ 
geneity  was  obtained  by  dividing  the  un -normal! zed  spectrum  of  the  middle 
third  of  a  data  sample  by  the  un-normalized  spectrum  of  the  entire  sample 
for  each  gust  velocity  component.  Homogeneity  is  indicated  when  the  spec¬ 
trum  ratios  are  1.0.  Each  turbulence  sample  for  which  power  spectra  were 
calculated  was  Individually  evaluated  for  homogeneous  characteristics. 

The  method  used  for  evaluating  these  characteristics  Involved  plotting 
the  spectrum  ratios  versus  spatial  frequency  (see  Appendix  IX  for  example 
plots)  and  visually  examining  each  plot.  Those  spectrun  ratio  values  which 
did  not  approximate  one  throughout  the  frequency  range  were  considered  to 
be  Indicative  of  nonhomogeneous  turbulence . 

Homogeneity  characteristics  of  the  turbulence  were  studied  by  calculating 
the  mean  value  of  each  spectrum  ratio.  Mean  values  which  were  less  than 
one  were  replaced  by  their  reciprocals.  Cumulative  probabilities  of  these 
mean  values  were  then  computed,  and  plotted  for  evaluation. 

Note  that,  if  the  shapes  of  the  spectra  are  the  same,  the  mean  of  the  ratio 
of  two  un-normalized  spectra  is  equal  to  the  ratio  of  the  variances  from 
each  spectrum.  That  is: 


♦<»>.  < 


(21.1) 


where:  subscripts  1  and  2  refer  to  the  middle  third  of  the 
data  sample  and  tl  "  entire  data  sample,  respectively. 

Since  essentially  all  the  homogeneity  characteristic  plots  (Appendix  IX) 
have  a  zero  slope,  the  two  spectra  are  either  coincident  or  parallel. 
Therefore,  it  is  concluded  that  the  ratio  of  the  variance  of  the  middle 
third  of  the  data  sample  to  the  variance  of  the  entire  sample  could  be 
employed  as  an  indicator  of  homogeneity  and  would  give  results  similar 
to  the  results  from  the  spectrum  ratio  method  which  was  used  throughout 
the  L0-L0CAT  program. 
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The  cumulative  probabilities  of  the  gpectrun  ratio  mean  values  are  shown 
for  all  three  components  In  Figure  21.1.  These  data  Include  the  turbulence 
samples  which  were  considered  to  be  both  homogeneous  and  nonhomo geneous . 

It  was  determined  from  the  static narity  run  test,  Section  5,  that  an  average 
of  77  per  cent  of  the  data  may  be  considered  to  be  stationary.  Based  on 
this  percentage,  it  can  be  seen  from  Figure  21.1  that  the  data  may  be  con¬ 
sidered  to  be  homogeneous  when  the  spectra  ratio  mean  is  between  1.0  and 
approximately  1.4.  When  the  ratio  is  greater  than  1.4,  the  data  are  non- 

homogeneous  according  to  the  run  test  criteria. 

\ 

As  discussed  above,  the  homogeneity  characteristics  of  the  L0-L0CAT  data 
were  evaluated  visually  from  the  spectrum  ratio  plots.  Only  the  data  for 
the  samples  considered  to  be  homogeneous  sure  plotted  in  Figure  21,2.  Ap¬ 
proximately  93  per  cent  of  these  data  are  between  the  values  of  1.0  and 
1.4.  Thus,  the  visual  method  of  homogeneity  determination  and  the  sta- 
tionarity  run  test  give  approximately  the  same  homogeneity  characteristic 
results. 

Figures  21.3  through  21.6  show  variations  of  homogeneity  as  functions  of 
the  categories  of  terrain,  absolute  altitude,  atmospheric  stability,  and 
tire  of  day.  The  data  contained  in  these  figures  are  those  from  botn  the 
homogeneous  and  nonhomogeneous  turbulence  samples.  All  of  these  figures 
are  comprised  of  the  same  data  and  vary  only  by  which  a  particulAr  cate¬ 
gory  component  is  held  constant.  The  data  in  these  figures  show  that  the 
turbulence  became  more  homogeneous  as:  (l)  terrain  roughness  decreased; 

(2)  atmospheric  stability  decreased;  and  (3)  time  of  day  progressed  from 
dawn  to  mid-afternoon.  Altitude  appears  to  have  very  little  effect  on 
the  homogeneity  characteristics  of  the  data. 
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22.  ISOTROPY 


A  turbulent  field  which  is  isotropic  is  one  in  which  there  is  no  preferred 
direction  and  therefore  all  average  functions  describing  this  field  must 
remain  unchanged  regardless  of  a  rotation  or  reflection  of  the  coordinate 
system.  Experimental  evidence  has  shown  that  most  localized  turbulence  is 
nearly  Isotropic  (local  isotropy). 

According  to  References  21.1,  18.2,  and  22.1,  the  isotropickcharacteristics 
will  follow  <frw(k)/|>v(k)  =  1  and  4>u(k)/<f>v(k)  will  vary  from^j.O  to  0.75  with 
frequency  and  scale  length.  The  characteristics  of  the  experimental  data 
were  established  by  dividing  the  ^normalized  power  spectra  of  u  and  w  by 
the  unnormalized  power  spectrum  of  v.  In  this  manner,  if  the  turbulence  is 
isotropic,  the  w/v  spectra  ratio  will  be  equal  to  1.0  arid  the  u/v  ratio  will 
approach  a  form  similar  to  that  in  Figure  22.1.  Figure  22.1  illustrates  how 
the  u/v  ratio  varies  with  the  scale  length  and  was  calculated  using  the  von 
Karntua  equations  for  turbulence  spectra. 

Unno  realized  spectra  ratios  were  determined  at  selected  spatial  frequency 
values  for  each  homogeneous  turbulence  sample.  The  selected  frequencies 
were  the  same  as  those  used  in  the  coherency  unalysis  (Section  20).  The 
mean  value  of  the  spectra  ratios  at  each  frequency  was  calculated  and  plot¬ 
ted  versus  frequency  for  selected  categories .  Isotropy  characteristics  for 
the  combined  categories  and  with  respect  to  terrain,  altitude,  stability, 
and  time  of  day  are  shown  in  Figures  22.2  through  22.6.  Isotropy  charac¬ 
teristics  for  the  individual  turbulence  samples  are  presented  in  Appendix 
IX. 


Theoretical  longitudinal  to  lateral  spectra  ratio  values  indicative  of  Iso¬ 
tropic  turbulence  were  determined  from  the  von  Kansan  expressions  (see  Sec¬ 
tion  28)  and  are  shown  in  Figures  22.2  through  22.6  for  comparison  purposes. 
The  von  Kanuan  spectra  ratios  shown  in  Figures  22.2  through  22.6  were  deter¬ 
mined  from  Figure  22.1  by  using  the  average  of  the  longitudinal  and  lateral 
scale  lengths  applicable  to  the  set  of  categorized  data  analyzed.  As  stated 
previously,  if  the  turbulence  is  isotropic,  then  the  vertical  to  lateral 
spectra  ratios  should  be  equal  to  1.0  throughout  the  frequency  range  regard¬ 
less  of  scale  length. 

Standard  deviations  of  the  ratios  from  the  mean  at  each  frequency  were  cal¬ 
culated.  These  values  based  on  all  of  the  vertical  to  lateral  ratios  are 
shown  in  Figure  22.2.  Standard  deviations  of  the  longitudinal  to  lateral 
ratios  are  not  shown.  This  ratio  varies  as  a  function  of  scale  length. 
Therefore,  standard  deviations  in  the  lower  frequency  range  are  increased 
due  to  the  scale  length  variations  (see  Figure  22. l)  and  do  not  give  a 
true  indication  of  variations  in  isotropic  characteristics. 

There  Is  a  trend  for  the  turbulence  to  be  slightly  more  isotropic  over 
rough  terrain  than  over  smooth  terrain  and  at  750  feet  compared  to  250 
feet.  There  is  very  little  difference  in  the  isotropy  characteristics 
of  the  data  recorded  during  the  different  a  ionospheric  stability  condi¬ 
tions.  The  turbulence  data  recorded  at  dawn  show  much  less  isotropy 
than  the  data  recorded  at  either  mid-morning  or  mid-afternoon.  It  should 
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be  noted,  however,  that  85  per  cent  of  the  dawn  data  analyzed  for  iso¬ 
tropic  characteristics  were  obtained  over  the  Peterson  route.  This  is 
because  the  majority  of  turbulence  samples  gathered  at  dawn  over  the 
McConnell  route  had  low  signal-to-noise  ratios  and  were,  therefore, 
not  included  In  the  analysis.  Also,  the  number  of  dawn  flights  over  the 
Edwards  route  was  limited  due  to  inclement  weather.  The  discussion  in 
Section  17  points  out  that  lateral  gust  velocities  at  Peterson  were 
larger  than  the  other  two  components  of  turbulence.  Thus,  the  vertical 
to  lateral  and  longitudinal  to  lateral  unnormalized  spectra  ratios 
from  the  Peterson  data  and,  consequently,  the  dawn  data  are  less  than 
required  for  isotropic  conditions. 
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23.  NORMALIZED  SPECTRA 

The  pover  spectra  were  normalized  by  dividing  the  pover  magnitude  by  L 
and  plotting  versus  kL.  Spectra  shift  19  and  down  in  power  with  a  re  spec¬ 
tive  increase  or  decrease  in  *tJ ,  and  change  shape  with  a  change  in  L. 
Therefore,  normalizing  in  this  manner  permitted  canparison  of  the  spectra 
for  different  categories  without  having  to  account  for  the  Influence  of 
and  L. 

Change  In  the  shape  of  the  normalized  spectra  with  category  variation 
was  evaluated  by  conqparlng  the  average  normalized  spectra  for  the  differ¬ 
ent  category  classifications.  The  average  spectrua  for  a  category  was  de¬ 
termined  in  the  same  way  that  average  coherency  functions  were  determined 
in  Section  20. 

2 

Spectra,  normalized  by  »t  L,  obtained  over  different  types  of  terrain  were 
compared  for  high  and  low  mountains,  and  rough  and  smooth  terrain.  Bough 
terrain  refers  to  the  combination  of  high  and  low  mountain  data  and  smooth 
terrain  refers  to  the  combination  of  plains,  desert,  and  water  data.  No 
consistent  or  significant  difference  in  spectra  shape  resulted  because  of 
terrain.  Figures  23.1  through  23.3  Bhov  comparisons  of  those  categories 
which  exhibited  the  maximum  differences. __ 

The  lateral  gust  velocity  spectra  in  Figure  23/6  exhibit  a  peculiarity 
that  is  also  present  to  a  leaser  degree  in  the  vertical  gust  velocity 
spectra.  This  appears  as  a  dip  at  a  kL  of  0.2/  and  is  believed  to  be  caused 
in  part  by  inccnplete  removal  of  probe  motiptlwhich  occurs  at  approxi¬ 
mately  0.3  cps.  This  0.3  cps  probe  motion,  accounting  for  variations 
in  ground  speed  and  scale  length,  falls  in  the  kL  range  from  approxi-  1 
raately  0.l4  to  0.^*6  kL.  An  example  plot  of  sideslip  differential  pres-  V 
sure  illustrating  the  probe  motion  is  presented  in  Figure  23. *+. 

The  removal  of  probe  motion  is  more  difficult  In  ^lateral  gust  velocity 
than  in  vertical  or  longitudinal  because  of  the  greater  amount  of  motion 
in  the  lateral  direction  in  this  particular  frequency  range. 

It  has  been  found  in  the  past,  Reference  1.2,  that  certain  errors  generated 
in  calculated  gust  velocities  are  a  function  of  the  unount  of  airplane 
motion.  Since  airplane  motion  and  hence,  probe  motion,  is  inherently 
larger  laterally  than  vertically  or  longitudinally,  lateral  gust  velocity 
is  inherently  more  inaccurate. 

Further  Indication  that  airplane  motion  contributes  to  the  dip  in  the 
-pectrc  yvr  obtained  by  comparing  the  normalized  spectra  of  low  intensity 
turbulence  samples  to  those  of  high  Intensity  turbulence  level.  The  aa- 
s’jsption  here  is  that  the  Lower  the  intensity  of  turbulence  the  less 
error  is  induced  by  airplane  motion. 

Figure  23.5  illustrates  this  cccBparisac.  The  data  represented  by  squares 
were  obtained^  from  10  samples  exhibiting  a  relatively  low  intensity  of 
turbulence  ( *t  *1.76  fpe)  but  still  sufficiently  high  so  as  not  to  be  ad¬ 
versely  affected  by  noise.  Tbe  solid  line  connecting  the  circles  repre¬ 
sents  the  average  of  over  100  high  intensity  turbulence  samples,  *  5.13 
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This  comparison  Indicates  the  samples  representing  the  low  turbulence 
intensities  are  in  better  agreement  with  the  results  from  Phases  I  and 
II  and  with  the  von  Kaxnan  spectra.  In  particular,  the  lover  intensity 
turbulence  data  are  in  better  agreement  with  the  von  Karmen  mathematical 
expressions  at  kL  values  below  0.1  than  with  the  average  of  the  high  tur¬ 
bulence  data. 

An  example  plot  showing  the  dip  in  the  spectra  for  an  individual  sample 
is  presented  in  Figure  23.6.  Numerous  samples  containing  this  type  spec¬ 
tra  shape  produced  the  effect  seen  in  Figure  23.2.  The  location  of  the 
dip  in  each  sample  shifts  due  to  different  ground  speeds  and  scale  lengths 
and  the  averaging  effect  which  produces  a  smoother  function  over  a  broad 
range  of  kL  values. 

2 

Conparison  of  average  spectra,  normalized  by  L,  to  determine  effects 
caused  by  altitude  are  presented  in  Figures  23.7  through  23.9.  Figure 
23.7  contains  a  comparison  for  longitudinal  gust  velocity.  No  signifi¬ 
cant  or  consistent  change  in  longitudinal  gust  velocity  can  be  attrib¬ 
uted  to  the  effect  of  altitude.  The  maximum  change  due  to  altitude  was 
noted  for  the  lateral  gust  velocity  component  as  shown  in  Figure  23.8. 

Here,  however,  the  difference  is  again  Influenced,  by  increased  airplane 
motion  effects  while  contour  flying  at  2 50  feet  as  compared  to  750  feet. 

The  difference  obtained  in  the  vertical  spectra  for  the  two  different  alti¬ 
tudes  is. presented  in  Figure  23.9*  The  vertical  gust  spectra  for  the  250- 
foot  data  are  consistently  lower  in  power  in  the  kL  range  of  0.1  to  0.4 
and  are  consistently  higher  at  kL  values  less  than  0.06.  The  difference 
in  the  spectra  in  the  kL  range  <*0001  0.1  to  0.4  is  similar  to  the  charac¬ 
teristic  previously  discussed  for  lateral  spectra.  The  existence  of  air¬ 
plane  motion  in  the  250-foot  altitude  data  is  the  result  of  a  generally 
higher  level  of  turbulence  at  250  feet  than  at  750  feet.  A  comparison  of 
low  Intensity  turbulence  to  high  intensity  turbulence  data,  both  obtained 
at  250  feet,  did  produce  indications  of  increased  probe  motion  in  the  kL 
range  of  0.06  and  less.  This  is  shown  in  Figure  23  10.  The  conclusion 
is  that  there  is  no  effect  duo  to  altitude  and  the  difference  Been  are  due 
mostly  to  probe  motion.  This  seems  valid  since  lew  frequency  power  in  vertical 
gust  spectra  would  not  be  expected  to  increase  with  decreasing  frequency. 

/ 

As  in  the  case  of  altitude,  no  consistent  or  significant  changes  in  the 
shape  of  the  spectra  were  found  to  be  caused  by  atmospheric  stability 
variations.  The  spectra  comparison  for  the  different  stabilities  is 
presented  in  Figures  23.11  through  23.13.  The  maximum  difference  ob¬ 
tained  in  the  comparison  for  time  of  day  effects  on  the  spectra  shape 
are  presented  in  Figures  23. l4  through  23.16. 

Plots  showing  the  scatter  bands  for  a  selected  category  combination  are 
presented  in  Figures  23.17  through  23.19*  These  data  show  the  mean, 
and  (m  ±  •)  obtained  by  the  averaging  process  previous]  y  discussed.  The 
scatter  bands  for  all  other  categories  analyzed  were  in  good  agreement. 

The  Irregularity  in  lateral  gust  spectra  is  shown  in  Figure  23.18. 

Plots  of  normalized  PSD's  for  the  turbulence  samples  exhibiting  valid 
data  for  all  three  gust  velocity  cceqponents  are  presented  In  Appendix  IX. 
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Figure  2j.9  Vertical  Power  Spectra  Variation  eltft  Altitude 
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24.  CROSS  SPECTRA 


The  cross  spectra  of  longitudinal  with  vertical  gust  velocity  and  longitudinal 
with  lateral  gust  velocity  were  calculated  primarily  for  use  in  the  computation 
of  coherency  functions  associated  with  the  turbulence  samples  for  which  power 
spectra  were  to  be  calculated.  A  limited  number  of  the  turbulence  samples  were 
plotted  in  the  forms  of  cross-spectra,  cospectra  and  quadrature  spectra  in  an  at¬ 
tempt  to  investigate  the  contribution  of  energy  from  each  component  tc  the  corre¬ 
lation  between  the  components  being  analyzed  and  to  possibly  define  an  "isotropic 
limit"  representing  the  wave  number  beyond  which  there  is  no  significant  contri¬ 
bution  to  the  Reynolds  stress.  It  is  suggested  in  Reference  18,2  that  both  the 


energy  contribution  and  the  isotropic  limit  may  be  obtained  from  cross-spectrum 
information.  /A 


According  to  Reference  13.2,  the  longitudinal -vertical  quadrature  spectrum  ' 
(imaginary  part  cf  the  cross  spectrum)  will  give  an  indication  as  to  what  part 
of  the  turbulence  eddies  are  being  passed  through  when  gust  velocities  are 


measured.  This  is  visualized  by  assuming  eddies  of  equal  wavelength  arranged 
as  shown  in  Figure  24.1.  If  these  eddies  are  penetrated  below  center,  the 
gust  velocity  time  histories  of  Figure  24.2  will  result.  If  they  are  pene¬ 
trated  above  center,  the  time  histories  of  Figure  24.3  result. 


In  direct  frequency  terms,  the  quad-spectrum  can  be  thought  of  as  the  average 
product  of  u  and  w1  within  a  narrow  frequency  interval  divided  by  the  interval. 
w1  is  w  displaced  to  the  left  90  degrees.  Note  that  when  w  is  displaced  90  de¬ 
grees  to  the  left,  passing  through  the  bottom  of  the  eddy  gives  a  negative  quad¬ 
rature  power  value  (u  and  w*  are  out  of  phase),  and  passing  through  the  top  gives 
a  positive  value  (u  and  w‘  are  in  phase). 

Representative  plots  of  the  cross-spectra,  cospectra,  and  quad-spectra  ob¬ 
tained  during  -.0-L0CAT  Riase  III  are  presented  in  Figures  24.4  through  24.15. 
Trends  of  the  quad-spectra  are  shewn  in  Figure  24.16.  The  trends  were  ob¬ 
tained  by  smoothing  each  longitudinal-vertical  quadrature  spectra  (Figures 
24.6,  24.9,  24.12,  and  24.15).  These  curves  were  then  averaged  for  the 
categories  noted.  All  curves  show  the  bottom  of  the  largest  eddies  being 
sensed,  on  the  average,  as  would  be  expected  since  the  eddy  radius  is  much 
larger  than  the  measuring  height. 

Mere  isotropy  is  shown  at  750  feet  than  at  250  feet,  as  was  found  by  the  analysis 
of  isotropy  in  Section  22.  Reference  18.2  suggests  an  isotropic  limit  of  1.7 
times  the  height.  Three  values,  1270  and  425  feet,  are  shown  in  Figure  24.16. 

The  curves  indicate  that  eddies  smaller  than  1270  and  425  'feet  were  sensed  both 
above  ana  below  their  centers,  giving  center  indications  on  the  average.  The 
fact  that  the  quadrature  power  of  750-foot  data  reaches  the  center  indication 
at  the  longer  wavelength  and  stays  there  with  decreasing  wavelength  indicates 
that  these  data  are  more  isotropic  than  the  250-foot  data. 


More  isotropy  iB  indicated  over  high  mountains  than  over  the  plains.  The  quad¬ 
rature  spectra  for  plains  deviates  from  a  center  indication  at  a  shorter  wave¬ 
length  than  the  quadrature  spectra  for  high  mountains.  The  plains  curve  has 
a  top  indication,  on  the  average,  for  wavelengths  between  1700  and  500  feet. 

This  indicates  that  wavelengths  from  500  feet  to  at  least  1700  feet  were  com¬ 
pressed  in  the  vertical  direction,  VteveiengthB  larger  than  1700  feet  for  all 
the  curves  were  probably  compressed;  however,  this  can  not  be  concluded  from  the 
•lata  since  750  feet  was  the  maximum  measuring  height.  These  results  agree  with 
those  >.btairiv.d  during  the  general  analysis  of  isotropy  in  Section  22. 
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Figure  24.1  Schematic  of  Eddies  Arranged  Contiguously 
Next  to  the  Ground 


Figure  2^. 2  Gust  Velocities  Encountered  vhen  Passing 
Tnrough  Contiguous  Eddies  Below  Center 


Figure  24.3  Gust  Velocities  Encountered  when  Passing 
Through  Contiguous  Eddies  Above  Center 


Figure  2U,i*  Cross,  Co,  sad  Quadrature  Power  Relationship 
for  Lateral  -  Longitudinal  Gust  Velocity 
Test  28,  Leg  6,  Category  a 1222^ 
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Figure  2k. 5  Cross,  Co,  and  Quadrature  °ower  Relationship 
for  Lateral  -  Vertical  Gust  Velocity  Test  28, 
Leg  6,  Category  412224 


Figure  24.6  Cross t  Co,  and  Quadrature  Power  I".  lationship 
for  Longitudinal  -  Vertical  Gust  Velocity 
Test  28,  Leg  6,  Category  412224 
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Figure  24.7  Cross,  Co,  and  Quadrature  Power  Relationship 
for  Lateral  -  Longitudinal  Oust  Velocity 
Test  51,  Leg  3,  Category  52?231 


10^ 


W-5 


nr* 


Figure  24,8  Cross ,  Co,  and  Quadrature  Power  Relationship 
for  Lateral  -  Vertical  Oust  Velocity  Teet  SI, 
Leg  3,  Category  32323* 


323 


Figure  24.9  Cross,  Co,  and  Quadrature  Peers r  Relationship 
for  Longitudinal  -  Vertical  Oust  Velocity 
Test  51 ,  !<•*  3.  Category  32323* 
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Figure  PU.li  Cross,  Co,  and  Quadrature  Power  Relationship 
for  Lateral  -  Vertical  Oust  Velocity  Teat  66, 
La«  2 ,  Category  312551 
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Fleur*  24.12  Cro«s,  Co,  and  Quadrature  Fowar  Relationship 
for  longitudinal  -  Vortical  O^at  Velocity 
Teat  66,  L*e  2,  Cat**ory  312331 
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Figure  24.1J  Cross,  Co,  snd  Quadrature  Power  Relationship 
for  Lateral  »  Long ltudinai  Oust  Veloc.  tjr 
Test  66,  Leg  6,  Category  U3331 
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Figure  2U,lk  Cross,  Co,  and  Quadrature  Power  Relationship 
for  Lateral  -  Vertical  Oust  Velocity  Tost  66 
Leg  6,  Category  113/51 
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Wave  Length 


2$.  TRUNCATED  STANDARD  DEVIATIONS 


The  truncated  standard  deviation  is  the  square  root  of  the  area  under  the 
truncated  spectrum.  The  lower  truncation  limit  of  the  Phase  III  spefctra 
in  cpa  was  increased  by  a  factor  of  two  over  that  oi'  Phases  I  and  11. 

This  was  done  to  make  the  lover  limit  in  terms  of  spatial  frequency  essen¬ 
tially  the  same  despite  the  approximate  factor  of  two  difference  in  speed 
between  the  Phases  I  and  II  aircraft  and 'the  Phase  in  aircraft. 

The  truncated  standard  deviation  was  calculated  as  follows: 

Phases  I  and  II 


r* 

I  <*>(k)  dk 


(25.1) 


Phase  III 


%I 

*(k)  dk 


(25.2) 


where:  Vj  and  Vjjj  are  the  ground  speeds  flown  during 

Phases  I  and  II  and  during  Phase  III,  respectively. 

It  should  be  noted  that  the  upper  limit  of  integration  was  held  constant. 
This  was  not  changed  because  the  area  under  the  truncated  spectrum  at  the 
higher  spatial  frequencies  contributes  only  a  small  amount  to  the  total 
area  under  the  spectrum.  However,  because  of  this,  the  truncated  stan¬ 
dard  deviations  obtained  from  samples  of  like  categories  would  be  expected 
to  be  slightly  higher  for  Phases  I  and  II  than  for  Phase  III  If  the  speed 
for  the  Phase  III  category  is  exactly  double  that  for  the  Phases  I  and  II 
c  ategory. 

Since  the  truncated  standard  deviation  is  calculated  from  that  portion  of 
the  spectrum  in  the  inert ial  subrange  defined  by  ♦(k)  -  Ak‘5/3  and  since 
Vjii  *»  2Vj,  the  relationship  between  equation  25.1  and  25.2  can  be  ex¬ 
pressed  as: 


0.965  9i„ 


(25.3) 
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Cumulative  probabilities  of  the  truncated  standard  deviations  from  Phases  I 
and  II  spectra  are  compared  in  Figures  25. 1  and  25.2  to  those  obtained  from 
the  Phase  III  spectra.  The  low  mountain  and  plains  category  data  were  uti¬ 
lized  for  this  comparison  because  these  data  were  considered  to  be  compara¬ 
ble  with  respect  to  the  geophysical  categories  for  the  different  phases  and 
consisted  of  enough  sables  to  be  considered  statistically  reliable.  The 
comparisons  show  that,  for  like  geophysical  categories,  the  truncated  stan¬ 
dard  deviation  values  obtained  from  Phases  I  and  II  spectra  are,  in  most 
instances,  equal  to  or  slightly  greater  than  those  obtained  from  Phase  III 
spectra. 


Figure  25.1  Truncated  Standard  Deviation  Variation  with  Phase  -  Low  Mountains 


26.  GUST  VELOCITY  DISTANCE  HISTORY  POWER  SPECTRA 


The  power  spectral  density  of  clear  air  turbulence  defines  the  energy  con¬ 
tent  of  the  turbulence  at  given  wavelengths  (l/k) .  The  wavelength  being 
in  terms  of  distance  traversed  over  the  ground  per  cycle  of  turbulence. 

In  relating  the  wavelength  of  the  turbulence  to  distance  traversed  over 
the  ground  rather  than  through  the  air  mass,  a  more  accurate  definition 
of  the  turbulence  is  obtained.  The  assumption  is  made,  of  course,  that 
the  turbulence  is  relatively  stationary  with  respect  to  the  ground  and 
does  not  move  with  the  air  mass.  This  assumption  seems  valid  since  the 
type  of  turbulence  being  meaeured  is  either  of  a  mechanical  or  convective 
nature  and  the  turbulence  generated  in  either  case  would  be  more  stationary 
in  relation  to  the  ground  than  to  the  air  mass. 

If  the  turbulence  is  generated  by  wind  flow  over  a  mountain  or  some  simi¬ 
lar  object  (mechanical),  the  turbulence  will  persist  until  the  protuberance 
in  the  wind  flow  corridor  is  eliminated  or  a  predominant  change  in  the  pre¬ 
vailing  wind  flow  occurs.  In  this  case,  the  turbulence  in  the  general  vi¬ 
cinity  of  the  mountain  will  not  move  with  the  air  mass.  If  the  turbulence 
is  generated  by  convection,  the  turbulence  will  remain  as  long  as  a  tem¬ 
perature  gradient  remains.  Since  the  gradient  exists  due  to  the  difference 
in  ground  and  air  temperature,  then  the  turbulence  must  remain  relatively 
stationary  with  respect  to  the  ground. 

This  leads  to  the  conclusion  that  turbulence  properties  may  be  defined 
more  accurately  if  studied  in  terms  of  distance  rather  than  time.  Some 
investigators  contend  that  if  the  airplane  speed  over  the  ground  does  not 
vary  significantly  throughout  the  data  collection  interval  that  one  may 
work  with  the  time  function  and  subsequently  convert  the  spectra  to  a  func¬ 
tion  of  spatial  (distance)  frequency.  This  hypothesis  was  proven  to  be 
true.  Power  spectra  of  each  of  the  three  components  of  gust  velocity  were 
calculated  for  four  samples  with  the  gust  velocity  being  in  terms  of  both 
time  and  distance.  These  samples  were  selected  on  the  basis  of  having 
large  variations  in  ground  speed,  being  homogeneous,  and  having  an  ade¬ 
quate  signal- to-noise  ratio. 

The  power  spectrum  as  a  function  of  spatial  frequency,  k,  was  obtained  by 
using  the  time  function  and  calculating  each  spectrum  in  terms  of  fre¬ 
quency,  f,  in  cps  and  then  converting  to  k  by  dividing  f  by  the  average 
ground  speed  for  the  L-l/2-minute  sample.  The  power  spectrum  was  also 
calculated  on  gust  velocity  distance  history  data.  The  gust  velocity  time 
history  was  converted  to  distance  history  using  instantaneous  ground  speed 
values  prior  to  spectrum  calculations.  The  resulting  spectra  were  in  terms 
of  spatial  frequency  k. 

These  normalized  spectra  for  gust  velocity  distance  history  and  time  his¬ 
tory  data  are  compared  in  Figures  26.1  through  2b. 4.  The  turbulence  sam¬ 
ples  cover  a  gust  velocity  ras  range  from  2,64  to  3*5$  f£s  with  scale 
lengths  varying  from  199  to  659  feet.  The  pertinent  data,  along  with  the 
maximum  speed  change  during  each  sample,  (AV),  is  presented  in  Table  2b.  1. 
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As  can  be  seen  from  the  spectra  comparisons  in  Figir  38  26.1  through  26.4, 
there  is  essentially  no  difference  in  the  spectra  even  for  a  sample  hav¬ 
ing  a  AV  =  106  fps.  It  should  be  noted,  however,  that  this  is  the  maximum 
change  and  for  this  particular  case  was  encountered  over  a  90-second  time 
interval  while  contour  flying  over  the  mountains. 

Although  scale  lengths  were  not  calculated  using  the  distance  histories, 
the  normalized  spectra  practically  overlay,  and  therefore,  the  scale 
lengths  would  not  differ.  This  good  agreement  indicates  that  bhe  use 
of  average  speeds  for  the  calculations  of  k  is  Justified. 


TABLE  26.1 


COMPARISON  OF  DISTANCE  HISTORY  AND 
TIME  HISTORY  STANDARD  EEVIATIONS 


Vertical 


Standard  Deviation  - 


Time 

distance 

Var 

History 

History 

u 

2.72 

2.72 

V 

2.78 

2.78 

V 

3.44 

3.44 

fps 


c 

6 

6 

6 


-  Distance  History 

- Time  History 


Figure  26. k  Time  History-Distance  History  Pover  Spectra  Comparison 
Test  87  Leg  7 


27.  SPECTRA  DURING  HIGH  SPEED  FLIGHT 


One  objective  of  the  Phase  III  Program  was  to  define  low  altitude  turbulence  to 
wavelengths  in  excess  of  21,500  feet  for  a  limited  number  of  turbulence  samples 
recorded  while  flying  at  the  maximum  speed.  The  airplane  maximum  speed,  however, 
.over  than  originally  expected  and  the  maximum  wavelength  measured  was  there¬ 
fore  only  16,000  feet.  The  average  wavelength  for  Phase  III  Is  approxi¬ 

mately  15,000  feet,  whereas  it  is  16,700  feet  for  the  turbulence  samples  recorded 
at  high  speed.  The  lover  maximum  speed  is  attributed  to  an  Increase  in  the  air¬ 
plane  drag  characteristics  caused  by  the  Installation  of  externally  mounted  in¬ 
strumentation. 

A  total  of  5^  turbulence  samples,  during  which  the  ground  speed  was  equal  to  or 
greater  than  675  fps,  were  selected  for  analysis.  These  samples  were  obtained 
over  Leg  6  at  both  the  McConnell  and  Peterson  routes  using  maximum  engine  thrust. 
These  legs  were  selected  to  give  maximum  flight  safety. 

Twenty-nine  of  the  54  selected  samples  exhibited  the  necessary  homogeneity  char¬ 
acteristics  required  for  a  valid  analysis.  The  time  series  rma  values  and  scale 
lengths  for  these  data  are  presented  in  Table  27.1.  The  average  rms  values  are 
3.69,  3.75*  and  3.30  fpa  for  u,  v,  and  v,  respectively.  The  corresponding  aver¬ 
age  scale  lengths  are  782,  502,  and  362  feet.  These  values  are  in  general  agree¬ 
ment  with  the  rest  of  the  Phase  in  data  recorded  over  the  plains. 

The  resulting  normalized  spectra  and  isotropy  and  coherency  characteristics  are 
presented  in  Figures  27.1  through  27.3.  These  daca  agree  in  general  with  the  rest 
of  the  Phase  III  data.  Lateral  gust  velocity  spectrum  values  obtained  at  high 
speed  are  compared  in  Figure  27.1  to  the  normalized  spectrum  obtained  at  250  feet 
above  the  plains  and  desert  terrain.  The  spectrum  for  250  feet  was  used  since  a 
high  percentage  (7 6ft)  of  the  high  speed  samples  were  obtained  at  250  feet. 

These  high  speed  samples  were  included  with  the  lower  speed  samples  for  statisti¬ 
cal  analyses  discussed  elsewhere  in  this  report. 


IABIE  27.1 

TURBULENCE  SAMPLES  HECOREED  AT  HIGH  SPEED 


Teat 

No. 

Leg 

No. 

Category 

No. 

teudmum 

Wavelength 

X  (feet) 

■ 

L 

k (feet) 

L_«I 

V 

w 

u 

V 

V 

24 

8 

414324 

16200 

4.22 

3.19 

695 

636 

239 

32 

8 

424324 

16400 

3.08 

3.38 

Hl8 

531 

718 

4o 

8 

412234 

17100 

2.07 

2.31 

2.50 

308 

335 

353 

46 

8 

412234 

16800 

2.46 

3.14 

2.68 

480 

400 

294 

161 

8 

414343 

16900 

3.00 

2.76 

2.59 

640 

435 

285 

162 

8 

422143 

16600 

2.98 

3.18 

... 

767 

750 

« 

163 

8 

424243 

16500 

4.07 

5.53 

4.02 

546 

895 

375 

i  1 66 

8 

411243 

16300 

2.94 

2.04 

-- 

1191 

321 

170 

8 

411143 

16900 

3.15 

-- 

267 

— 

171 

8 

414243 

17000 

5.35 

5.37 

4.33 

850 

467 

336 

172 

8 

413343 

17300 

4.55 

4.80 

— 

722 

701 

— 

185 

8 

414343 

16200 

3.72 

3.51 

3.36 

796 

419 

269 

187 

8 

423243 

16500 

2.47 

2.85 

2.20 

777 

673 

452 

188 

8 

423343 

16300 

3.49 

3.65 

3.61 

964 

687 

547 

190 

8 

413343 

16600 

4.51 

5.37 

— 

649 

768 

— 

209 

8 

424243 

16500 

2.87 

2.58 

2.54 

1078 

483 

469 

211 

8 

412143 

16900 

3.02 

3.37 

2.67 

251 

216 

231 

212 

8 

4lte43 

17100 

5.80 

5.71 

4.64 

838 

470 

318 

213 

8 

413343 

16500 

5.51 

4.99 

3.85 

861 

338 

244 

220 

8 

413313 

16500 

4.09 

3.59 

3.4l 

1180 

382 

381 

224 

8 

413213 

16600 

3.23 

3.34 

3.29 

489 

356 

300 

225 

8 

414313 

16700 

5-53 

5.58 

3.93 

1218 

766 

331 

230 

8 

412213 

16600 

2.85 

2.69 

2.55 

592 

328 

238 

231 

8 

414313 

16500 

4.65 

3.92 

2.88 

15^7 

612 

327 

233 

8 

423213 

16800 

4.26 

4.47 

4.57 

968 

639 

664 

236 

8 

413213 

16300 

3.15 

2.86 

3.31 

664 

257 

346 

24l 

8 

413213 

17500 

3-35 

3.31 

3.36 

693 

374 

300 

246 

8 

414213 

17400 

3.35 

3.57 

3.03 

797 

4?3 

316 

255 

8 

413213 

16500 

3.48 

3.38 

3.35 

729 

351 

362 

<*00/o 


Figure  27.1  High  Speed  Flight  Guet  Velocity  Power  Spectre 
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Average  Coherency  for  High  Speed  Flight  Sauries 


28.  spectra  comparison 


Use  of  the  power  spectral  density  approach  in  aircraft  design  predicates 
the  need  for  a  mathematical  model  of  the  PSD  curve.  Many  mathematical 
expressions  have  been  proposed  to  represent  the  atmospheric  turbulence 
spectra.  The  expressions  most  widely  vised  in  the  aircraft  industry  at 
the  present  time  include  expressions  suggested  by  Theodore  von  Karman 
(Reference  28. l)  and  H.  L.  Dryden  (Reference  28.2).  These  expressions, 
as  shown  below,  vary  according  to  the  gust  velocity  component  being  rep¬ 
resented. 

von  Karman  Vertical  and  Lateral  - 


[2  +  377.5(lk»  fc)2] 
[l  +  70.78  (LXvk)2]u/s 


(28.1) 


von  Karman  Longitudinal 


KW  I  i+LKu 

L  '*!  J,’  [l-70.78(LKuk)2]5/‘ 

Dryden  Vertical  and  Lateral  - 

r*y00l  _  LPv[g  +  6(2*LPvk)23 

1  Jd  [l+(2»LDvk)2]2 

Dryden  Longitudinal  - 

[*M]  4LPu 

't*  JD  * 


(28.2) 


(28.3) 


(28.4) 


As  suggested  In  Reference  28.3,  the  equations  may  be  represented  by  the 
following  general  forms: 

Vertical  and  Lateral  - 


fry  (k)  ^  2L[l  +  8*2a2  k2(n  +  l)} 

9t,  (1  +  4r2&2  k2)  n  *3/2 


(28.5) 
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(28.6) 


Theoretical  spectra  defined  by  Equations  28.1  through  28.8  were  compared 
to  the  experimentally  determined  spectra  from  the  Phase  III  program. 
These  comparisons  were  made  by  dividing  each  experimental  spectrum  by 
the  mathematical  expression  and  plotting  this  ratio  vermis  spatial  fre¬ 
quency.  A  ratio  of  1.0  indicates  perfect  agreement  between  the  experi¬ 
mental  and  the  mathematically  defined  spectra.  These  plots  are  shown 
in  Appendix  IX. 

The  von  Karman  expressions  were  analyzed  throughout  the  entire  Phase  III 
program.  The  Dryden  expressions  were  analyzed  for  data  obtained  at  the 
McConnell,  Edwards,  and  Peterson  routes.  Ihe  Lixaley-Panofsky  and  Buah- 
Pacofsky  expressions  were  evaluated  for  data  obtained  at  the  Griffiss 
location. 

3W 


Average  values  of  these  experimental  to  mathematical  ratios  were  computed 
and  plotted.  This  was  done  by  determining  the  value  of  the  ratios  for 
each  turbulence  sample  at  certain  values  of  spatial  frequency.  These 
ratios  vere  then  averaged  at  each  of  the  chosen  frequencies.  The  results 
are  shown  in  Figures  28.1  through  28.3*  These  figures  show  the  mean  values 
and  the  standard  deviation  of  the  ratios  about  these  mean  values.  All  of 
the  mean  values  shown  in  Figures  2b. 1  through  28.3  are  plotted  in  Figure 
28.4  for  purposes  of  comparison. 


L4I 


- ±  One  Standard  Deviation 


Longitudinal  Component 


Lateral  Component 


Vertical  Component 


J _ _ L _ I _ L..J _ L  Li.  1 _ l  1  11111 


10-3 

SPATIAL  imuiENCTY  -  cpf 


Figure  28.2  Comparison  of  Experimental  and  Dryden 
Power  Spectra 


f 


Figure  28.3  Comparison  of  Experimental  Pow«.r  Spectra  with 

Lumley-Panofsky  and  Busch-Panofsky  Power  Spectra 
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29.  REPRESENTATIVE  SPECTRA  SHAPE  1 

4 

The  power  spectral  density  of  turbulence  has  been  found  to  be  a  function  t 
of  spatial  frequency  to  the  -5/3  power  in  the  inertial  subrange.  This  ,'i 
is  Illustrated  In  the  Phases  1  and  II  data  and  by  the  Phase  III  data.  ;  J 

Because  of  this,  the  von  Karman  mathematical  expressions  give  a  good  I 

representation  of  all  components  of  the  turbulence  spectra  within  this  j 

frequency  range.  In  the  low  frequency,  long  wavelength  portion  of  the  3 

spectra,  the  vertical  and  lateral  expressions  show  sqoa  deviation  from  j 
the  measured  data.  It  appears  this  is  due  in  part  to  difficulties  encoun-t? 
tered  in  accounting  for  all  the  vertical  end  lateral  nation  of  the  air-  (!' 
plane  at  low  frequency.  Longitudinal  spectra  agree  well  with  the  von  / 
Karman  mathematical  expression. 

The  correlation  between  von  Karman  and  experimental  spectra  Is  shown  in 
Figure  28.1.  The  disagreement  in  thfl  low  frequency  portion  of  the  lat-  ;; 
eral  and  vertical  spectra  where  the  von  Karman  expression  underestimates  , 
the  spectra  power  is  discussed  further  in  Sections  46  ana  47. 

The  Dryden  expressions  show  a  greater  amount  of  disagreement  with  the 
experimentally  determined  data,  especially  at  the  higher  frequencies 
(Figure  26.2)  where  the  Dryden  expressions  have  a  -2  logarithmic  slope  / 
rather  than  the  -5/3  slope  of  the  data.  At  the,  lower  frequencies,  the  j 
Dryden  expressions  show  the  same  type  of  disagreement  is  those  suggested 
by  von  Kansan  in  that  they  underestimate  the  power  of  the  spectra.  This 
underestimation  is  slightly  more  for  the  Dryden  expressions,  however. 

The  Lu&ley-Panofaky  expression  for  the  longitudinal  coaponent  of  turbu¬ 
lence  shows  good  agreement  with  the  experimentally  determined  data  at 
the  higher  frequencies  (Figure  28.3).  However,  at  the  lower  frequency 
values  this  expression  deviates  from  the  date  significantly.  In  this 
area  the  mathematical  expr'.ssi'>n  underest imtes  the  experimentally  de¬ 
termined  spectra  by  a  factor  of  approximately  2.5. 

The  Busch-Panofsky  mathematical  expressions  for  vertical,  and  lateral 
spectra  are  in  good  agreement  with  the  experimental  data  (Figure  28,3). 

This  agreement  extends  throughout  the  frequency  range  investigated  ex¬ 
cept  for  the  vertical  component  at  the  highest  end  of  the  frequency  range. 

At  the  Mgfr  frequencies  the  experimentally  determined  vertical  spectra 
showed  some  deviation  from  the  theoretical  -5/3  slope  (see  Sections  46 
and  47) .  At  the  low  fi\  quency  end  of  the  spectra  these  expressions  show 
the  best  agreement  of  t„ose  studied. 

i 

The  results  discussed  k**re  for  the  von  Kansan  and  Dryden  expressions  are 
similar  to  those  obtained  during  LD-LOCAT  Phases  I  and  T1  (Reference  1.2). 
The  Ltmley-Panofsiy  xpression,  which  was  also  evaluated  during  Phases  I 
and  II,  shows  more  disagreement  with  the  experimentally  determined  spectra 
at  the  lower  frequencies  for  the  Phase  XU  daua.  The  Phase  III  data  ex¬ 
tend  to  lower  frequency  values  than  the  Phases  I  and  U  data  because  of 
the  higher  speed  aircraft  used  as  the  instrumentation  platform.  The 
Phases  I  and  II  data  extended  to  a  low  spatirl  frequency  of  approximately 
1,5  x  10  **  cpf,  whereas  tne  Phase  IJJ  lata  in  Figures  28.1  through  28.4 
are  plotted  to  a  low  frequency  value  at  J  x  ID  cpf.  Thu  Busch -Panof sky 
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expressions  were  not  analyzed  during  Phases  I  and  II.  Instead,  an  expres¬ 
sion  suggested  by  U.  0.  Lappe  (Reference  l4.l)  was  Investigated  and  the 
results  showed  the  sane  type  of  disagreement  as  that  discussed  above  for 
the  Dryden  expressions.  v 

It  is  felt  that  the  low  frequency  portion  of  the  spectrum  should  be  stud¬ 
ied  further  through  the  use  of  higher  speod  aircraft.  This  would  pro¬ 
vide  a  more  accurate  record  of  the  long  wavelength  data.  Spectra  models 
based  on  the  data  from  all  three  phases  are  shown  and  compared  to  the 
von  Karmen  spectra  in  Section  46. 
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30.  TURBUUEMCE  SCALE  LENGTH 

In  order  to  use  either  the  von  Ka/.aan  or  Dryden  mathematical  expressions 
(Equations  28.1  through  2.8.4),  a  value  for  the  integral  scale  of  turbu¬ 
lence,  L,  must  be  known.  L  is  defined  ae  the  area  under  the  autocorre¬ 
lation  curve: 


(30.1) 


and, 


(30.2) 


where: 

f(,  r)  =  longitudinal  autocorrelation  function 

g(r)  =  lateral  or  vortical  eutocorrelation  function 

The  magnitude  of  the  integral  3cale,  or  ncale  length,  is  an  indication  of 
the  average  eddy  dice.  It  can' be  sho'n  that  these  scale  lengths  in  the 
three  orthogonal  directions  are  related  by  LX«2L/»2LZ  fpr  isotropic 
turbulence  and  that  Lu  -  ix,  Lv  *=  2Ly,  and  L„  =  2L7>.  It  should  be  pointed 
out  that  Equations  28.1  through  28.4  sit  written  in  terms  of  longitudinal 
scale  length.  Thus  Lu,  L..,  and  Lw  are  longitudinal  scale  lengths  calcu¬ 
lated  from  the  longitudinal,  lateral,  and  vertical  gust  velocity  compo¬ 
nents,  respectively. • 

Scale  lengths  may  1  omputed  from  the  autocorrelation  function  if  Taylor’s 
hypothesis  is  valid  (Reference  22. l) .  They  may  also  be  confuted  from  the 
spectrum  in  the  following  manner;  Standard  deviation  has  been  defined  as 
the  square  root  of  the  area  under  the  spectrum.  For  the  computation  of 
scale  length,  a  truncated  standard  deviation,  <rT,  covering  only  the  inertia], 
subrange  of  the  spectrin  is  defined  in  Section  25  as: 


(30.3) 


k3  *  10  /  V 

For  calculations  of  L  using  the  von  K&rmon  express! one.  Equations  28.1  and 
28.2  were  simplified  by  noting  that,  in  the  Inertial  subrange,  377*5  (hK)2 
and  70. 7®  (Lk)2  »  1.  Thus,  the  onl'1'-  terms  remining  in  the  numerators 
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P 

Jkg 


*(k)  dk 


1/2 


where: 


k2  *  0.667/  V 


L>  - 1  -  - j 


g(r) dr 


/ 


f  (r)  dr 


i 

i 

i 

\  and  denominators  of  these  expressions  were  those  terms  containing  L  and/or 

<  k.  The  resulting  expressions  for  <t>  (k),  from  either  Equation  28.1  or  28.2 

'  were  substituted  in  Equation  30.3*  The  resulting  expression  was  integrated 

;  and  solved  for  scale  length,  yielding  Equations  30.4  and  30.5,  respectively. 

Longitudinal  Scale  Length  -  Vertical  or  Lateral  Component 


(30.4) 


Longitudinal  Scale  Length  -  Longitudinal  Conponent 


I<Ku  *  0.0717 


(30.5) 


Scale  lengths  used  in  .the  Bryden  expressions  were  determined  in  the  same 
manner.  Equations  28.3  and  28.4  were  simplified  by  noting  that  (2  w LK) 2 
»  1  in  the  inertial  subrange.  The  resulting  equations  obtained  by  in¬ 
serting  these  simplified  equations  into  Equation  30.3  and  integrating  were: 

Longitudinal  Scale  Length  -  Vertical  or  Lateral  Conponent 


t”'0,lsa  (£).("£■ - 1)  (306) 


Longitudinal  Scale  Length  -  Longitudinal  Conponent 


(30.7) 


Equations  30.4.  and  30.6  are  expressions  for  the  longitudinal  scale  length 
as  computed  from  the  vertical  or  lateral  gust  velocity  components  while 
Equations  30,5  and  30.7  representive  longitudinal  scale  length  as  computed 
from  the  longitudinal  component  of  gust  velocity. 

As  mentioned  previously,  scale  lengths  could  also  be  calculated  from  the 
area  under  the  gust  velocity  autocorrelation  function.  Basically,  the 
autocovariance  function  R(r)  was  normalized  by  dividing  it  by  R(0)  whore: 
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*(k)  dk  -  <rt2 


(30.8) 


B(0) 


•/; 


thus  yielding  the  autocorrelation  function  R(t)/b(0).  The  longitudinal 
scale  length  could  then  be  computed  from  the  longitudinal  component: 


L 


[R(r)/R(0)]dr 


(30.9) 


The  area  was  multiplied  by  ground  speed  to  convert  from  units  of  time  to 
units  of  distance.  Equation  30.9  1 s  for  the  calculation  cf  lon¬ 
gitudinal  scale  length  from  tne  longitudinal  component.  If  calculated 
from  either  of  the  transverse  components,  the  ucale  length  obtained  from 
Equation  30.9  oust  be  multiplied  by  two  to  get  longitudinal  scale  length. 


During  the  L0-L0CAT  Program,  scale  lengths  weru  computed  from  the  power 
spectrum,  as  shown  by  Equations  30.4  through  30.7'  However,  autocorrela¬ 
tion  functions  were  calculated  for  a  few  randomly  selected  turbulence 
samples,  and  scale  lengths  were  compiled  using  Equation  30.9*  The  scale 
lengths  calculated  by  the  two  different  methods  were  compared. 

A 

The  autocorrelation  functions  from  one  of  the  sangples  used  for  these 
scale  length  calculations  are  shown  in  Figures  30.1  through  30.3.  The 
areas  under  the  curves  were  manually  computed  to  the  point  where  the 
function  crosses  zero.  This  assures  that  the  function  oscillates  about 
zero  from  that  point  to  infinity  with  positive  and  negative  areas  tend¬ 
ing  to  cancel  each  other.  The  functions  obviously  cannot  be  calculated 
out  to  infinity.  Likewise,  calculating  them  out  to  greater  and  greater 
lag  times  means  decreasing  the  value  of  &t,  the  distance  between  fre¬ 
quency  samples,  and  thus,  decreasing  the  confidence  level  of  the  data. 
Therefore,  it  is  not  known  whether  these  particular  functions  do,  in 
fact,  oscillate  about  zero. 

Comparisons  of  scale  lengths  calculated  by  the  two  methods  are  shown  in 
Figures  30.4  through  30.6.  Scale  length  values  agree  well  when  calculated 
from 2 the  longitudinal  oon?>onent  of  turbulence.  However,  scale  lengths 
computed  from  the  lateral  and  vertical  components  differ  accordlr  g  to  the 
method  used.  As  can  be  seen  in  Figures  30.5  and  30.6,  the  values  obtained 
from  the  autocorrelation  functions  are  somewhat  larger  In  magnitude  than 
those  obtained  from  the  power  spectra. 

Both  of  these  methods  are  based  on  the  assumption  that  the  von  Karmen  ex¬ 
pressions  describe  gust  velocity  spectra  shapes.  Since  the  autocorrela¬ 
tion  function  is  the  Fourier  transform  of  the  power  spectrum,  it  would 
seem  that. the  differences  in  scale  lengths  from  the  two  methods  of 
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calculation  occur  because  the  area  under  the  autocorrelation  function 
beyond  where  the  function  crosses  zero  could  not  be  considered. 

Reference  22.1  suggests  another  method  for  calculating  scale  length  from 
the  autocorrelation  function.  This  method  considers  only  the  point  at 
which  the  autocorrelation  function  crosses  zero,  i.e.: 

L  =  0.to2  rc  (30.10) 

where: 

rc  =  t  c  V 

and  rc  =  value  of  t  at  which  autocorrelation  function  crosses  zero 

Comparisons  of  scale  lengths  calculated  using  this  method  and  using  the 
power  spectra  method  are  shown  in  Figures  30.7  through  30.9.  It  should 
be  noted  that  when  Equation  30.10  was  used  for  calculation  of  scale  length 
from  the  vertical  and  lateral  autocorrelation  functions,  the  re suit b  were 
multiplied  by  two  to  correspond  to  the  longitudinal  scale  length. 

Compared  to  scale  lengths  calculated  from  the  area  under  the  autocorrela¬ 
tion  function.  Equation  30.10  showed  less  agreement  with  the  power  spectra 
method  for  the  scale  length  from  the  vertical  and  lateral  components  and 
about  the  same  agreement  for  scale  length  from  the  longitudinal  conqponent. 

Still  another  method  for  calculating  scale  length  was  suggested  in  Refer¬ 
ence  28.4.  This  method  also  involved  use  of  the  power  spectra.  Scale 
length  was  calculated  by  multiplying  the  power  density  of  the  spectra  by 
spatial  frequency,  k,  and  plotting  this  product  versus  l/k.  The  scale 
length  was  defined  as  being  proportional  to  the  value  of  l/k  at  which 
k<f>(k)  is  a  maximum. 

L  =  A  L  (30.11) 

K 


where: 


A  =  constant  of  proportionality  K 

2 

kB  =s  value  of  k  where  k+(k)/*t  is  a  maximum 

The  von  Karman  spectra  expressions  were  used  to  determine  A.  Equations 
28.1  and  28.2  were  multiplied  by  k,  differentiated,  set  equal  to  zero, 
and  solved,  for  L.  The  resulting  express  ions  for  scale  length  were: 


LKu  =  0.1456  ~~  (from  longitudinal  component)  (30.12) 

*#> 

LKv=  0.2119  i-  (from  lateral  and  vertical  components)  (30.13) 


Scale  lengths  calculated  by  this  method  are  compared  to  those  calculated 
from  Equations  30.4  and  30.5  in  Figures  30.10  through  30.12.  The  agree¬ 
ment  for  scale  length  from  the  longitudinal  component  is  good.  Values  do 
not  agree  well  for  scale  lengths  computed  from  the  vertical  and  lateral 
components,  however.  This  may  be  due  to  excessive  airplane  motion  included 
in  the  vertical  and  lateral  gust  velocities  (See  Section  23).  The  method 
of  Equations  30.12  and  30.13  involved  use  of  spectra  values  at  the  rela¬ 
tively  low  frequencies  while  the  scale  lengths  calculated  by  the  method  of 
Equations  30.4  and  30.5  were  a  function  of  the  spectra  in  the  inertial  sub¬ 
range  where  the  data  were  considered  to  be  more  accurate. 

o 

Another  problem  that  arose  concerning  the  vise  of  Equations  30.12  and  30.13 
was  the  difficulty  in  determining  a.  value  for  l/kB.  The  turbulence  samples 
chosen  for  making  these  comparisons  *were  the  same  as  those  used  for  the  com¬ 
parisons  in  Figures  30.4  through  30. 9»  However,  not  all  of  these  samples 
could  be  used.  For  seme  of  these  samples,  there  was  no  distinct  peak  in 
the  value  of  k*(k)/#j;  ,  so  that  an  accurate  value  of  l/kB  could  not  be 
determined.  The  values  of  the  spectra  were  such  that  the  maximum  value 
of  k*(k)/#*  in  this  frequency  range  extended  over  a  large  range  of  l/km 
values.  This  occurred  in  over  half  of  the  turbulence  samples  selected 
and  these  samples  were  not  Included  in  the  analysis. 

Since  the  von  Karmen  expressions  gave  a  good  representation  of  the  experi¬ 
mentally  determined  spectra  for  all  three  components,  the  scale  lengths 
from  these  expressions  were  analyzed  statistically  by  category.  The  cumu¬ 
lative  probability  distributions  of  these  scale  lengths  for  all  the  Phase 
IU  data,  regardless  of  category,  are  shown  in  Figure  30.13* 

The  von  Karman  scale  lengths  are  ihown  as  functions  of  terrain,  absolute 
altitude,  atmospheric  stability,  and  time  of  day  in  Figures  30. l4  through 
30.17.  The  general  trend  of  the  data  in  these  figures  indicates  that  scale 
length  tends  to  increase  with  increasing  terrain  roughness,  increasing  abso¬ 
lute  altitude,  and  decreasing  atmospheric  stability. 

Scale  lengths  recorded  at  dawn  are  greater  than  those  recorded  during  the 
other  two  times  of  day  when  calculated  from  the  lateral  component.  When 
calculated  from  the  longitudinal,  and  vertical  components,  the  probabilities 
at  dawn  are  less  for  the  smaller  values  of  scale  length,  but  tend  to  become 
greater  than  at  mid-morning  or  mid-afternoon  at  the  larger  scale  length 
values.  This  seems  somewhat  inconsistent  with  *t  statistics.  For  a  given 
airspeed  (therefore,  given  values  of  k2  and  k3),  scale  length  is  a  function 
of  (  "J  #T)3  (see  Equations  30.4  and  30.5).  As  shown  in  Section  13,  *t 
increases  from  dawn  to  mid-morning  to  mid-afternoon.  Thus,  it  would  seem 
that  the  smallest  values  of  Beale  length  would  occur  at  dawn.  The  reason 
why  this  is  not  so  requires  consideration  of  how  the  turbulence  samples 
were  chosen  for  the  scale  length  calculations.  As  previously  discussed 
in  Section  19,  the  power  spectra  were  edited  and  analyzed  according  to 
the  attributes  of  the  data.  Low  intensity  turbulence  samples  were  gen¬ 
erally  honhcxaogeneous  and  exhibited  high  coherence  characteristics  which 
resulted  in  turbulence  spectra  of  generally  poor  quality.  Therefore, 
these  data  were  not  considered  valid  for  spectral  analysis.  Many  of 
these  low  intensity  turbulence  samples  were  obtained  at  dawn  since,  as 
previously  mentioned,  the  intensity  of  turbulence  at  this  time  of  day  was 
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generally  less  than  at  the  other  times  of  day.  However,  the  turbulence 
samples  obtained  at  Peterson  over  the  high  mountains  have  been  shown  to 
consist  of  gust  velocities  greater  than  at  the  other  locations.  Many  of 
the  dawn  turbulence  samples  obtained  at  Peterson  were  of  great  enough 
magnitude  to  be  used  for  spectral  analysis.  Thus,  approximately  85  per 
cent  of  the  dawn  turbulence  sauries  chosen  for  spectral  analysis,  and 
consequently,  scale  length  calculation,  were  obtained  at  Peterson.  The 
data  obtained  at  mid-morning  and  mid-afternoon  were  more  evenly  distrib¬ 
uted  among  the  four  locations.  This  caused  the  gust  velocities  in  the 
dawn  category  to  be  higher  than  at  the  other  times  of  day  for  the  power 
spectra  turbulence  samples.  The  values  of  »t  used  for  scale  length  cal¬ 
culations  are  shown  as  functions  of  time  of  day  in  Figure  30. 18. 

The  conclusions  concerning  variation  of  scale  length  with  time  of  day 
should  be  made  with  respect  to  the  considerations  discussed  above.  The 
distributions  of  and  <rT  for  all  PSD  samples  (not  only  those  chosen  for 
analysis)  tend  to  indicate  that  scale  length  values  increase  from  dawn  to 
mid-morning  to  mid-afternoon. 

.  * 
During  Phases  I  and  II  of  the  L0-L0CAT  Program,  (Reference  1.2),  scale 
lengths  at  the  250-foot  altitude  were  found  to  be  primarily  a  function 
of  terrain,  increasing  with  greater  terrain  roughness.  At  the  750-foot 
altitude,  atmospheric  stability  appeared  to  have  the  greatest  effect  on 
scale  lengths.  As  the  atmosphere  became  more  unstable  the  scale  lengths 
showed  a  tendency  to  increase  in  magnitude. 

The  Phase  III  data  were  investigated  for  the  same  trends.  Figures  30. 19 
and  30.20  show  the  variations  of  scale  length  with  terrain  for  data  ob¬ 
tained  at  each  altitude.  Scale  lengths  at  the  250-foot  altitude  tend  to 
increase  with  increasing  terrain  roughness.  Those  recorded  at  750  feet 
show  the  seme  tendency,  but  to  a  somewhat  lesser  degree.  Scale  lengths 
are  shown  as  a  function  of  absolute  altitude  and  atmospheric  stability 
in  Figures  30.21  and  30.22.  The  only  distinct  variation  that  occurs  is 
for  scale  lengths  calculated  from  the  longitudinal  component  at  the  250- 
foot  altitude.  The  reason  for  the  very  stable  data  showing  higher  prob¬ 
abilities  at  the  larger  values  of  scale  length  is  analogous  to  the  situa¬ 
tion  existing  in  the  scale  length  analysis  with  respect  to  time  of  day. 
Otherwise,  no  consistent  or  predominant  trends  are  noted  in  these  data. 

Thus  it  would  seem  that,  of  the  categories  investigated,  the  primary  vari¬ 
ations  ir>  scale  length  are  a  result  of  changes  in  type  'of  terrain  and  the 
altitude  above  the  terrain.  The  effects  of  terrain  at  each  altitude,  dis¬ 
cussed  above,  are  shown  in  Figures  30.19  and  30.20.  The  effect  of  altitude 
over  each  type  of  terrain  is  shown  in  Figures  30.23  through  30.25.  Scale 
lengths  from  the  vertical  component  seem  to  be  affected  by  altitude  slightly 
more  than  those  calculated  from  the  other  two  components.  Larger  scale 
length  values  were  recorded  at  the  higher  altitude.  The  results  are  sum¬ 
marized  in  Figure  30.26  where  mean  scale  length  values  are  shown.  A  sum¬ 
mary  of  the  altitude-stability  analysis  is  shown  in  Figure  30. 27.  There 
were  not  enough  scale  length  data  obtained  over  the  desert  or  water  to 
permit  a  statistical  analysis. 
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The  following  conclusions  are  included  in  Reference  30.1  for  the  behavior 
of  longitudinal  scale  length  from  the  vertical  component  as  a  function  of 
altitude  and  stability: 

•  Scale  length  increases  with  altitude  for  atmospheric  stability 
conditions  other  than  very  stable. 

•  Under  very  stable  conditions,  scale  length  is  nearly  constant 
with  altitude. 

With  respect  to  these  conclusions,  the  scale  length  from  the  vertical  com¬ 
ponent  does  show  a  substantial  increase  with  increase  in  altitude  for  the 
stability  conditions  other  than  very ^jjtable  (Figure  30.27).  However,  a 
similar  relationship  is  seen  for  $Ke  very  stable  conditions,  although  the 

effects  of  altitude  appear  to  be<  iess. 

/ 

The  scale  lengths  from  the  l^gitudinal  and  lateral  components  do  not  show 
these  same  variations  with  Alt it udn  for  the  various  stability  conditions. 
These  scale  lengths  show  the  greatest  increases  with  increased  altitude  for 
a  stable  atmosphere.  Altitude  change  has  no  effect  on  these  values  for 
neutral  stability.  Scale  lengths  from  the  longitudinal  component  had 
greater  mean  values  for  the  data  at  250  feet  than  for  the  data  at  750 
feet  for  the  very  stable  and  unstable  atmospheric  conditions. 

According  to  Reference  30.2,  the  scale  lengths  for  clear  air  turbulence 
below  an  altitude  of  2,500  feet  for  use  in  von  Kannan  equations  are  to 
be  calculated  as: 

Vertical  Component 

L)Cy  *  H  (30.1k) 

Longitudinal  and  Lateral  Components 

Lku*  LKv  ■  18k  H^3  (30.15) 

It  should  be  noted  taat  different  equations  are  used  for  the  w  and  v  com¬ 
ponents.  Therefore,  in  the  following  discussion,  the  subscript  v  is  used 
to  signify  the  lateral  component  and  w  the  vertical  component. 

A  comparison  between  the  scale  lengths  reconmended  by  these  equations  and 
those  obtained  during  the  LO-LOCAT  Phase  III  Program  are  shown  in  Table 
30.1.  V 
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TABLE  30.1 


L0-L0CAI  PHASE  ITL  SCALE  IEROTHS  FOR  VON  KAEMAN  EXPRESSIONS 
COMPARED  TO  THOSE  RECOMMENDED  BY  MIL-F-OO8785A  (USAF) 


Absolute 

Altitude(H) ~  Ft. 

Reference  30.2 
Recommended 

Scale  length  ~  Ft . 

Average  LO-L0CAT 
Phase  m 

Scale  length  ~  Ft. 

iai 

1 

mgm 

m 

n 

L„ 

K  w 

250  • 

1159 

1159 

250 

715 

569 

394 

750 

1670 

750 

•714 

632 

508 

With  the  exception  of  scale  length  computed  from  the  vertical  component, 
scaie  lengths  calculated  from  Equations  30.14  and  30.15  over  estimate 
the  LO-LOCAT  Phase  III  experimentally  determined  data. 

Although  the  Dryden  expressions  do  not  fit  the  LO-LOCAT  data  as  well  as 
the  von  Kansan  express!  jns,  the  Dryden  expressions  are  used  in  aircraft 
design.  Reference  30.2  states  that  the  use  of  the  Dryden  expressions  is 
permissible  when  it  is  not  feasible  to  use  the  von  Karmen  forms.  Accord¬ 
ing  to  this  reference,  the  Beales  of  turbulence  to  be  used  in  the  Dryden 
equations  for  flight  below  1,000  feet  are: 

lDu=h  (30.16) 

lDu  =  LDv  =  100H1/3  (30.17) 


Table  30.2  showB  how  the  Phase  III  data  compare  with  values  calculated 
from  Equations  30.16  and  30.17* 


TABLE  30.2 


IO-DOCAT  PHASE  III  SCALE  LENGTHS  FOR  DRYTEN  EXPRESSIONS 
COMPARED  TO  m>SE  RECOMMENDED  BY  MIL-F-OO8785A  (USAF) 


Absolute 

Altitude (H)  ~  Ft. 

Reference  30.2  - 
Recommended 

Scale  length  ~  Ft. 

BB— 

^Du 

131 

iJD  v 

sm 

Sv 

LDv 

250 

630 

630 

250 

509 

502 

389 

750 

'909 

909 

750 

512 

538 

_ 

465 
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The  Dryden  scale  lengths,  with  the  exception  of  those  calculated  from  the 
vertical  component,  are  over  estimated  by  the  Reference  30.2  equations,  as 
was  the  case  with  the  von  Karman  scale  lengths.  Equations  30.  l4  through 
30.17  do  indicate  an  Increase  in  scale  length  for  a  corresponding  increase 
in  altitude.  This  trend  was  also  noticed  in  the  L0-L0CAT  data.  Likewise, 
the  Reference  30.2  equations  agree  with  the  LO-IOCAT  data  in  that  the  scale 
length  from  the  vertical  conqponent  is  more  dependent  on  altitude  than  are 
the  scale  lengths  determined  from  either  the  longitudinal  or  lateral  com¬ 
ponents. 

The  cumulative  probability  distributions  of  the  Dryden  scale  lengths  are 
shown  Jn  Figcne  30.28. 

Under  isotropic  conditions: 

Lu  =>  Lv  =*  Lw  (30. l9) 

However,  Reference  13. 1  also  states  that  es  the  ground  level  is  approached, 
the  scale  length,  as  computed  from  the  vertical  component,  decreases.  Ref¬ 
erence  13.1  suggests  the  following  expressions  as  -the  approximate  relation¬ 
ships  between  scale  lengths  calculated  from  the  different  components  of 
turbulence: 

For  50  s  H  <  500  feet: 

Lu  Lv 

=  17  15  1*30  -  0.0006H  (30.19) 


For  E  >  500  feet: 


(30.20) 


For  absolute  altitude  levels  of  250  feet  and  750  feet,  these  expressions 
yield  ratios  of  1.15  and  1.00,  respectively.  The  values  are  compared  to 
values  from  the  DO -DOC AT  Phase  III  experimentally  determined  data  in  Table 

30.3. 


TABLE  30.3 

D0-D0CAI  PHASE  III  TCALE  LENGTH  RATIOS  COMPARED  TO 
THOSE  RECOMMENCED  BY  REFERENCE  13.1 _ 


DD-TOCAT  Phase  III  Data' 


Absolute 

Altitude  (g)  ~?c. 


250 


Values  Predicted  By 
Equations  30»19  and  30*20 


1.15 


'Hu 


'K« 


1.82 


h(Cr 

L~ 

K  v 


1.44 


'Du 


Dv 


1.31 


Lev 

Ldv 
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750 


1. 00 


1.40  1.25 


1.10 


1.29 

1.16 


Although  "Che  ratios  computed  from  the  LO-LOCAT  Phase  III  data  are  some¬ 
what  larger  t'-'an  those  predicted  by  Equations  30.19  and  30.20,  the  fact 
that  the  ratios  e  *e  reduced  in  magnitude  with  increasing  altitude  above 
the  terrain  is  apparent.  The  ratios  computed  for  the  Dryden  scale  lengths 
shot  better  agreement  with  the  values  predicted  by  Equations  30.19  and 
30.20  than  do  the  von  Karman  ratios.  If  the  same  form  as  Equations  30.19 
and  30.20  is  used,  approximations  for  the  Phase  III  relationships  would  be 
as  follows: 


t-Xu  ^Kv 

- —  =  —  =  1.73  -  0.0005H 


1  Kv 


'  Kv 


(30.21) 


^Du 

Ldv 


Ld, 

I'D  v 


=  1.39  -  .00034H 


(30.22) 


These  are  only  approximations.  Equations  30.21  and  30.22  assume  a  linear 
relationship  between  the  scale  lengths  at  250  and  750  feet  and  were  de¬ 
rived  by  assuming  tt  following  average  values: 


JKu 


LKv 
LKv 


7—  =  -  =  1.63  at  250  feet 

L  Ky 


(30.23) 


L  Ku 

L~ 


Kv 


1  Kv 


=  1.32  at  750  feet 


(30.24) 


JDu 


■Dv 


J->Dv  / 

-  =  1.30  at  250  feet 

L  Dv 


(30.25) 


L  Du 

4  Dw 


'Dv 


'Dw 


=  1.13  at  750  feet 


(30.26) 


The  values  were  obtained  by  averaging  the  appropriate  values  in  Table  30.3. 


According  to  Luppe  i^i  Reference  30.3  the  type  of  terrain  nay  be  taken  into 
consideration  by  'using  the  following  formula  for  calculating  L: 

L  =  h0  +  H  Lh  •  (30.27) 

where  Lh  and  hc  are  defined  a a  follows: 
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Terrain  Claes 


Smooth 

~/3 

Low  Mountains 

1/2 

High  Mountains 

1/8 

Values  from  Equation  30,27  are  compared 
lengths  in  Table  30.4. 

135 

300 

675 

to  Phase  IH  von  Kansan  scale 


t 

f 


Scale  lengths  obtained  over  the  plains  legB  vere  used  for  the  smooth  ter¬ 
rain  classification  in  Table  30.4. 

TABLE  30.4 

IO-IOCAT  PHASE  III  3CAIE  LENGTHS  COMPARED  TO  THOSE 


RECOMMENDED  BY  AFFDL-TR -67-122 


Absolute 

Altitude(H)  -  Ft. 

Type  of 
Terrain 

L  From 
Equation 
30.27 
~  Ft. 

■ 

■E53S5!s±EE9B 

mm 

— T  ^ 

Lr 

Kv 

250 

Smooth 

402 

655 

456 

291 

250 

Lo.  Mount. 

425 

671 

518 

4ao 

250  j 

Hi.  Mount, 

706 

775 

636 

452 

750  1 

Smooth 

635 

64o 

572 

503 

750 

Lo.  Mount. 

675 

697 

608 

476 

750 

Hi.  Mount. 

769 

_ 

765 

663 

520 

It  should  be  pointed  out  that  L,  as  calculated  in  Reference  30.3  vas  de¬ 
signed  for  use  in  the  mathematical  spectra  expression  suggested  therein, 
and  not  for  use  in  the  voni  Kansan  equations.  The  Reference  30.3  mathe¬ 
matical  expression  vas  evaluated  during  Phases  I  and  II  and  it  vas  found 
that  experimental  data  did  not  correlate  veil  with  that  expression.  How¬ 
ever,  the  fact  that  terrain  is  considered  in  the  calculation  of  L  is  in 
harmony  with  the  results  obtained  during  the  IC-LOCAT  Progrea.  The  re¬ 
sults  from  Equation  30.27  agree  better  with  the  experimentally  determined 
data  than  did  results  from  Equations  30. 14  and  30.15,  where  only  altitude 
is  considered. 

The  average  values  of  the  von  Karmen  and  Lryden  scale  lengths,  regardless 
of  geophysical  category,  are  shown  in  Table  30.5*  It  should  be  kept  in  mind 
that  these  mean  values  are  biased  due  to  the  uneven  distribution  of  data  in 
the  different  categories. 
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TABUS  30.5 


LO-LOGAT  PEASE  III  AVERAGE  SCALE  LENOTHS 


— —  11  j4  1  1'  i'll 

t. 

m 

LKv 

LKw 

J'Du 

T 

Dv 

m 

598 

450 

5H 

519 

427 

In  Section  28,  mathematical  expressions  (Equations  28. 7  and  28.8)  for  the 
PSD  suggested  by  Busch  and  Panofsky  and  by  Lumley  and  Panofsky  were  analyzed 
Scale  length,  as  such,  is  not  contained  in  these  expressions.  Instead,  Cv 
is  used  in  the  Busch-Panofsky  expression  and  C  is  used  in  the  Lumley- 
Panofsky  expression.  These  parameters  were  calculated  in  a  manner  similar 
to  the  von  Karnan  and  Dryden  scale  lengths.  Equations  28.7  and  28.8  were 
simplified  by  noting  that  in  the  inertial  subrange  1.5(Cvk)5/3  »  1,  and 
(2950  k)5/3  »  1,  respectively.  If  these  expressions  Eire  substituted  into 

the  truncated  standard  deviation  expression  (Equation  30.3)  and  integrated, 
the  following  result  is  obtained:  * 


(30.28) 


(30.29) 


Thus,  C  and  Cv  are  related  to  scale  length  and,  from  the  preceding  equa¬ 
tions  (30.28  and  30.29)  and  Equations  30.4  and  30. 5,  the  following  rela¬ 
tionships  were  derived; 

Cv  »  4.70  LKv  (30.30) 

C  .  0.415  (L  )1/3 

f.  u 


The  cumulative  probability  distributions  of  nhese  parameters  are  presented 
in  Figures  30.29  through  30. 3i.  In  Reference  18.2,  a  value  of  2.4s  is  sug¬ 
gested  for  C.  For  the  L0-LDCAT  Phase  III  Program,  a  mean  value  of  3.54  was 
determined.  The  mean  values  cf  C,  calculated  from  the  vertical  and  lateral 
components  are  1805  and  2243,  respectively. 

It  should  be  pointed  out  that  all  data  analyzed  within  this  discursion  were 
derived  from  only  homogeneous  turbulence  samples  (see  Section  19).  Scale 
lengths  fer  lew  frequency  data*are  4A«cuas«i  in  Section  52t> 
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T  -  Seconds 


Figure  30.3  Example  Vertical  Gust  Velocity 
Autocorrelation  Function  Plot 
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Figure  30.5  Comparison  of  von  Kansan  and  Autocorrelation 
Function  Derived  Longitudinal  Scale  Lengths 
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Figure  30.6  Comparison  of  von  Karmen  and  Autocorrelation 
Function  Derived  Longitudinal  Scale  Lengths 
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Figure  30.8  Coap&rison  of  -on  Karaan  and  Autocorrelatioi 
Function  Derived  Longitudinal  Sc  ale  Lengths 
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Figure  30.11  Comparison  of  von  Kantian  ana  l/k  Derived 
longitudinal  Sc^le  Lengths 
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Figure  jO.13  /on  Karajan  Longitudinal  Scale  Length 
Cuarulative  Probability 
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Figure  30.15  von  Karaan  longitudinal  Scale  Length 

Cumulative  Probability  Associated  with 
Absolute  Altitude 
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Figure  30.18  Effects  of  Tie*  of  Day  on  Cumulative 

Probability  of  Gust  Velocity  EHS  Values 
M  Used  in  Scale  Length  Calculations 
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Figure  30.21  Effect#  of  Atmospheric  Stability  on 
CuRuLativa  Probability  of  von  JCanaan 
Longitudinal  Scale  Length#  from  Eat  a 
Recorded  at  250  feet 


Calculated  from 
Longitudinal  Component 


Calculated  from 
Lateral  Component 


kso 


-Li-l-U  j,  t  1  i  1  i  1  .  t 
fay  flop  woo  mo  i*oq  l6oo  iaoo  aoou 

nu  mn  -  r*at 


.005 

.oos 


>0001  L  j-  I  l  1  l  1  1  1  l  1  I  1  »  I 

0  200  too  600  doo  WOO  1200  ltoo 


nil  UBRUH  -  Feet 


Calculated  from 
Vertical  Component 


A 

.  i .  i .  i ...  i .  i .  i  .  j  .  > 

«c  tto  to  un  ubo  itoo  itoe  1M6  moo  •> 


figure  30.2*  Effect*  of  Absolute  Altitude  on  Cumulative 
Probability  of  von  Karmn  Long?"  udJLnal 
-*  Scale  Iangth*  from  Lata  Recorded  Over 
Lev  Mountain* 
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Figure  30*25  Effects  of  Absolute  Altitude  on  Cumulative 
Probability  of  too  Karma  Longitudinal 
Scale  Lengths  frco  Data  Recorded  Over 
Plains 
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Figure  30.26  Variations  of  Mean  von  Karman  Longitudinal 
Scale  Length  with  Absolute  Altitude  and 
Terrain  Type 
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Figure  30,27  Variations  of  Mean  von  Karman  Longitudinal 
Scale  Length  vlth  Absolute  Altitude  and 
Atmospheric  Stability 
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Figure  30.28  Dryden  Longitudinal  Scale  Length 
Cumulative  Probability 
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Figure  30.29  Cumulative  Probability  of  L>»l*y-Panofaky 
Scale  Parameter,  C,  a a  Calculated  froa 
tbe  Loogituiiaal  Cocqpo&wt 


31.  TUKBUIZNCE  MICROSCALES 


The  integral  scale  length  is  defined  as  "being  equal  to  the  area  under  the 
gust  velocity  autocorrelation  curve.  The  Taylor  microscale  can  also  "be 
defined  from  the  autocorrelation  curve  as  the  intercept  on  the  y  cods  of 
a  parabola  fitted  to  the  vertex  of  the  curve.  The  Taylor  microscale  rep¬ 
resents  the  average  eddy  size  where  viscous  forces  dissipate  turbulent 
energy. 

<4 

Taylor  microscales  were  calculated  using  the  following  equation  (Refer¬ 
ence  18.2): 


X  =  15  v 


(31.1) 


where  v  Is  the  kinematic  viscosity  and  <  is  the  viscous  dissipation  rate. 
Viscous  dissipation  rates  pertaining  to  Phase  III  data  are  discussed  in 
Section  3^» 

The  cumulative  probability  distributions  of  Taylor's  microscale  for  the 
time  of  day  categories  are  shown  in  Figure  31.1.  As  has  been  the  case 
in  the  past  (Reference  31.l)>  the  distribution  is  shown  to  be  very  close 
to  a  log-normal.  There  is  a  slight  difference  in  the  distribution  for 
dawn  than  for  other  times  of  day.  The  distributions  by  altitude  and 
stability  in  Figure  31. 2  show  a  definite  difference  between  very  stable 
and  the  other  stabilities.  It  appears  that  the  microscale  in  the  dissi¬ 
pation  range  increases  as  altitude  increases  and  ns  the  atmosphere  becomes 
more  stable.  Ugh  mountains  and  plains  distributions  shown  in  Figure  31.3 
indicate  that  the  microscale  does  not  vary  significantly  between  smooth 
and  rough  terrain. 

Reference  31*1  suggests  that  the  turbulence  spectra  should  be  a  function 
of  the  kinematic  viscosity  and  the  rate  at  which  the  energy  is  dissipated. 
A  scale  length  formulated  in  the  form  of: 


is  referred  to  as  the  Kolmogorov  microscale .  This  length  is  much  smaller 
than  Taylor’s  microscale  and  represents  an  eddy  size  where  viscous  forces 
predominate*  dissipating  turbulence  fluctuations  into  heat.  Figure  31.4 
Shows  the  cumulative  probability  of  the  Xoloogcrov  microscale  for  the  time 
of  day  categories  indicating  slight  changes  with  time  of  day.  Figure  31.5 
gives  the  probability  as  a  function  of  altitude  and  stability*  shoving  that 
the  values  increase  slightly  as  altitude  increases  and  stability  decreases. 
Terrain  effects  shown  in  Figure  31.6  indicate  a  trend  toward  longer  micro- 
scales  far  plains  type  terrain. 


Reference  18.2  Indicates  that  values  on  the  order  of  10  cm  (0.328  feet)  and 
approximately  1  m  (0.00328  feet)  for  Taylor  and  Kolmogorov  microscales,  re¬ 
spectively,  are  to  be  expected  in  the  atmosphere  "near  the  surface."  Table 
31.1  summarizes  the  average  values  as  determined  during  Phase  III. 

TABLE  31.1 


AVERAGE  TAXLOR  AMD  KDIHX30R0V  MICROSCALES 


TQtitu3e 

Taylor 

iColmogorov 

(Feet) 

(Feet) 

(Feet) 

250 

0.84 

0.0033 

750 

0.90 

0.0036 

Comparison  of  ftiaae  III  data  with  Biases  I  and  II  data  reveals  that-  there 
is  a  20  per  cent  increase  in  the  average  Taylor  microscale  for  data  re¬ 
corded  atTu^Q  feet  during  Phase  HI.  This  can  be  primarily  attributed  to 
the  hi^h^jpust  velocity  standard  deviations  recorded  during  Hiaae  III. 

The  average  Taylor  microscales  obtained  during  Phase  HI  are  compared  in 
Figure  31*7  with  those  presented  by  Taylor  in  Reference  31-2  and  by 
MacCready  in  Reference  31.3  as  well  as  with  those  obtained  during  Phases 
I  and  II  (Reference  1-2),  the  BSffll  Twrer  fly-bys  (Reference  Section  4l), 
and  the  flights  supporting  the  Operation  Rough  Rider  Project  (NSSL). 

the  average  Kolmogorov  microscales  obtained  during  Hiase  III  are  compared 
in  Figure  31.8  with  these  presented  in  Reference  18.2  by  Luraley-Panofsky 
as  well  as  with  those  obtained  during  Biases  I  and  II,  the  BREN  Tcwer  fly¬ 
bys,  and  the  flights  supporting  the  Operation  Rough  Rider  Projyct  (NSSL). 
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Fig-ore  31.1  Variation  of  Taylor  Microscale  for  Time  of  Day 
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Figure  31.2  Variation  of  Taylor  Microscale  for  Altitude 
^  and  Atmospheric  Stability 
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Figure  31. ^  Variation  of  Kolmogorov  Microscale  for  Time  of  Day 
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Figure  31.7  Variation  of  Taylor  Midroecale  with  Altitude 
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Figure  31,8  Variation  of  Kolmogorov  Microscale  with  Altitude 
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